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On a large class of lattices (such as the sawtooth chain, the kagome and the pyrochlore 
lattices) the quantum Heisenberg and the repulsive Hubbard models may host a com¬ 
pletely dispersionless (flat) energy band in the single-particle spectrum. The flat-band 
states can be viewed as completely localized within a finite volume (trap) of the lattice 
and allow for construction of many-particle states, roughly speaking, by occupying the 
traps with particles. If the flat band happens to be the lowest-energy one the manifold 
of such many-body states will often determine the ground-state and low-temperature 
physics of the models at hand even in the presence of strong interactions. The localized 
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nature of these many-body states makes possible the mapping of this subset of eigenstates 
onto a corresponding classical hard-core system. As a result, the ground-state and low- 
temperature properties of the strongly correlated flat-band systems can be analyzed in 
detail using concepts and tools of classical statistical mechanics (e.g., classical lattice-gas 
approach or percolation approach), in contrast to more challenging quantum many-body 
techniques usually necessary to examine strongly correlated quantum systems. 

In this review we recapitulate the basic features of the flat-band spin systems and 
briefly summarize earlier studies in the field. Main emphasis is made on recent develop¬ 
ments which include results for both spin and electron flat-band models. In particular, 
for flat-band spin systems we highlight field-driven phase transitions for frustrated quan¬ 
tum Heisenberg antiferromagnets at low temperatures, chiral flat-band states, as well as 
the effect of a slight dispersion of a previously strictly flat band due to nonideal lattice 
geometry. For electronic systems, we discuss the universal low-temperature behavior of 
several flat-band Hubbard models, the emergence of ground-state ferromagnetism in the 
square-lattice Tasaki-Hubbard model and the related Pauli-correlated percolation prob¬ 
lem, as well as the dispersion-driven ground-state ferromagnetism in flat-band Hubbard 
systems. Closely related studies and possible experimental realizations of the flat-band 
physics are also described briefly. 

Keywords : Heisenberg model; Hubbard model; flat band; Mielke-Tasaki ferromagnetism. 

PACS numbers: 71.27.+a, 75.10.-b, 75.10.Jm, 75.10.Lp, 71.10.Fd 


1. Introduction and Outline 

Flat-band physics attracts increasing attention nowadays. A good evidence for this 
is a number of reviews which have been p ublished recently, mainly with a special fo¬ 
cus on topological flat-band models However, the interest in strongly correlated 
flat-band systems came up long ago. For example, it was known since early 1990s 
that linear spin-wave theory for the antiferromagnetic Heisenberg model on the 
kagome lattice leads to flat modes!®* This flat mode of the kagome antiferromagnet 
corresponds to a coherent motion of spins on a hexagon while the spins surround¬ 
ing the hexagon are ordered ferromagnetically. 0 However, these calculations were 
approximate (linear spin-wave theory, i.e., large-s limit) and refer to the zero-field 
case. Several years later new prominent features of quantum Heisenberg antiferro¬ 
magnets on a class of highly frustrated latticed were observed around the satura¬ 
tion field. These are the ground-state magnetization jump at the saturation fieldj^a 
magnetic -field driv en spin-Peierls instability,® a finite residual entropy at the satura¬ 
tion fieldjnHHHH an{ j an unconventional low-temperature thermodynamics ! 11112| 13 ] 
A crucial feature leading to these effects is the fact that the lowest band of one- 
magnon excitations above the fully polarized ferromagnetic state (which becomes 
the ground state for magnetic fields above the saturation field) is completely dis¬ 
persionless (flat). In many cases the flat-band states can be visualized as the states 
located within a small part of the lattice (trap) and this localization in turn is a 
consequence of destructive quantum interference which occurs due to a special ge- 


a Magnetic interactions are frustrated, if a spin cannot arrange its orientation such that it profits 
from the interaction with its neighbors as, for instance, in the case of antiferromagnetic interactions 
on the triangular lattice. For a more in depth discussion we refer to Refs. ElO 
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ometry of the lattice. The very existence of localized magnons as the lowest-energy 
one-particle states opens an interesting perspective to construct and fully character¬ 
ize many-magnon ground states of the considered frustrated quantum Heisenberg 
antiferromagnets. Moreover, the set of relevant low-energy many-magnon states of 
the quantum spin systems can be mapped onto corresponding (highly non-trivial) 
classical systems and as a result one may use completely different methods from 
the toolbox of classical statistical mechanics to study the spin systems at hand. 
The striking difference t o the previous linear spin-wave studies of flat-band anti¬ 
ferromagnets at zero fielJEl 

is that there are flat-band many-body states at high 
magnetic fields which are exact ground states for arbitrary spin quantum number 
s and that the abovementioned features are true quantum effects which disappear 
in the classical limit s —> oo. 

Later on it was recognized! ^ that the flat-band quantum spin systems are 
covertly related to another class of strongly correlated electron models, the so- 
called flat-band Hubbard ferromagnets, which were discussed even earlier j 15| 16| 17 
While the seminal papers by A. Mielke and H. Tasaki are focused mainly on rigorous 
proofs of the ferro mag netic ground states for the flat-band Hubbard modefGIISI-E 
and their stability)!^ similarities to localized-magnon systems lead to the discus¬ 
sion of the ground-state degenera cy and universal low-temperature behavior of flat- 
band Hubbard systemsJJHEniH Furthermore, a percolation representation for the 
ground-state ferromagnetism suggested by A. Mielke and H. Tasaki^ has been 
recognized to present a ne w class of cla ssical percolation problems - the so-called 
Pauli-correlated percolation.^2 23 24 2o 

In the theory of frustrated magnetisrrP*? 27 28 ^he localized-magnon systems 
represent a special class of models which in a special regime (high fields, low tem¬ 
peratures) admit a rather detailed study of their properties which is based on the 
rigorous knowledge of the many-body localized-magnon states. In the broad field 
of the theory of the Hubbard modelpOMQI anc g j n particular, of ferromagnetism of 
the Hubbard modelpl 32 33 ^he Mielke-Tasaki flat-band ferromagnetism provides 
one specific root to ferromagnetism which also allows rather rigorous analysis. Both 
fields appear to be related: From the mathematical point of view the description 
of spin and electron systems has many similarities, since a key property of both 
flat-band systems is that relevant quantum many-body degrees of freedom can be 
mapped on classical ones. 

Earlier studies on localized magnons were reviewed some years ago, see Refs. [34], 
f35l and [36j In the present paper we intend to provide a snapshot of the state of 
art in the field presenting new directions with larger focus on electron systems. 
Naturally, the main attention is paid to summarizing studies in which the authors 
are involved or closely related ones. We also provide a detailed introduction to the 
basic concepts. 

The outline of this review is as follows. We begin with a short introduction 
to localized-magnon and localized-electron systems focusing on lattices, Hamilto¬ 
nians, concepts etc. and recall the main results which were reviewed earlier (mag- 
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netization jump, la ttice instability, residual entropy, low-temperature peak of spe- 
cific heat) l ’*4 | ao | ob | gee g ec . [ 2 ] Concerning new results for spin systems, we illus¬ 
trate the emergence of Isin g d egrees of freedom in frustrated ladder and bilayer 
Heisenbe rg ant iferromagnets^ in Sec. 0 chiral localized magnons (and localized 
electronsj2S139 j n Se c. [4j and discuss the effect of small deviations from the flat- 
band geometrj® II j n Sec. [ 5 ] Next we discuss some studies stimulated by the 
localized-magnon picture (Sec. HJ) discussing in particular flat-band systems at den¬ 
sities slightly above the localized-magnon-crystal density42i and the magn etization 
plateaus in the spin-i Heisenberg antiferromagnet on the kagome lattice.SSEIESl 
Very recently it has been found that localized-magnon states can be the relevant 
low-energy states even in zero ma gnetic field for a sawtooth chain with ferro- and 
antiferromagnetic interactions We will briefly discuss this very new example for 
localized-magnon states in Sec. [ 6 ] Next, in Sec. 0 we deal with solid-state sys¬ 
tems which should exhibit localized-magnon physics. The natural mineral azurite 
CU 3 ( 003 ) 2 (OH ) 2 is the m ost promising exa mple for experimental solid-state obser¬ 
vation of flat-band physics P ^^^ I SO I Sl^ l-yyjg illustrate the theoretical predictions 
based on a diamond-chain Heisenberg model. Another localized-magnon compound 
has been synthesized recently, Ba 2 CoSi 2 C> 6 Cl 2 .E^ We also briefly discuss this system 
in Sec. 0 Note, however, that these compounds do not fulfill the strict flat-band 
criteria, but show some characteristic features such as the magnetization jump. 

In the second part of this review we move to electron systems. In Sec. 0 we 
consider two different classes of the flat-band Hubbard model an d discuss their uni¬ 
versal behavior at low electron densities and low temperatures 11SI2S.2 X Section 0 
reviews recent work on the existence of the ground-state ferromagnetism in a par¬ 
ticular two-dimensional flat-band Hubbard model which can be described by a new 
type of percolation called Pauli-correlated percolation.^] j n Sec. 03 we consider 
Hubbard systems for which a flat band does not lead to ground-state ferromag¬ 
netism, but the ground state becomes ferromagnetic, if the flat band acquires small 
dispersion and the on-site Coulomb repulsion is sufficiently large (dispersion-driven 
ferromagnetism)EH Furthermore, we discuss several related topics (Sec. [IT]). First 
we give in Sec. 0L] a brief account on a general construction sche me of a class of 
exact ground states of specific strongly correlated electron systems^^^Then we 
discuss some recent results on flat-band systems with randomness®® 5fi|59] ag we n 
as the transport properties of flat-band two-terminal nanodevices.®^ In Sec. [T 2 l we 
discuss possible experimental realization of theoretical scenarios for flat-band Hub¬ 
bard models l 61 [ 62|MlM165|6616 7 |68|69|. ZQI FinaUy, in Sec. 02 we close by a summary 
and sketch a perspective for future work. 
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2. Flat Band and Localized Eigenstates 

2.1. One- and many-particle localized eigenstates in frustrated 
Heisenberg antiferromagnets 

We begin with a brief account of the concepts used to discuss the localized-magnon 
states. We consider the Heisenberg model in an external magnetic field 

N 

H = ^2 J ijSi - Sj - hSZ , S z = J2s Z i- ( 1 ) 

«) i=l 


> 0 denotes 


Here s,; = (sf,sf,sf) is the spin-i operator attached to the site i, Ji 3 
the antiferromagnetic exchange interaction between nearest-neighboring sites, h is 
the dimensionless external magnetic field, the first sum in Eq. © runs over all 
neighboring sites of the N-site lattice. S z commutes with the Hamiltonian H and 
hence the z-component of the total spin is a good quantum number. In the subspace 
with S z = ^ there is only one state, that is the fully polarized state |FM) = | f 
... t) which plays the role of the vacuum state for magnons. In the subspace with 
/S ' 2 = — 1 we have N one-magnon states, and their explicit form and their energy 

can be found exactly. The one-magnon energies crucially depend on the lattice 
geometry. 

Let us mention here that the concept of localized-magnon states is not restricted 
to spin quantum number s = \ and isotropic Heisenberg systems, see, e.g., Refs. [81 
mi However, the related physical effects (see below) are most pronounced in the 
extreme quantum case s = i and vanish in the classical limit s —> oo. 

We consider lattices which support a completely dispersionless lowest-energy 
magnon band. There is quite a lot of lattices having flat bands. In Figs. [T] and [2] 
we show some of them: The sawtooth lattice)"^ the kagome chain}©] the frustrated 
diamond chainpH the frustrated two-leg ladder,^- the double-tetrahedra chain}©! 
and the frustrated (cylindrical) three-let { ladder"' 7 "'" (one-dimensional lat tices) and 
the Tasaki lattice}©! the kagome lattice}©!© the square-kagome lattice©HMJ anc j 
the frustrated bilayeJ^ (two-dimensional lattices). Other lattices which were also 
discussed in the localized-magnon co ntex t are: The dimer-plaquette chain;®®] the 
star lattice,^! the checkerboard lattice}^ the sorrel netJ^Hthe pyrochlore lattice!©! 
There are also lattices which where not considered in the localized-magnon cont ext 
so far, although, they belong to this class, e.g., the triangulated kagome lattice .®"' 7 
Traces of a flat band can be observed also for magnetic molecules, see, e.g., Refs. [581 
[851 We mention, that most of the flat-band lattices were considered also in the gen¬ 
eral context of frustrated quantum magnetism. For example, the spin-i Heisenberg 
antiferromagnets on the kagomej©!©! the starj©land the square-kagome®^ lattices 
are known as examples of two-dimensional quantum spin systems with a magneti¬ 
cally disordered ground state. 

To construct more flat-band lattices one may use Mielke’s line-graph 
prescription®^ or Tasaki’s decoration procedure!®® In this sense we can define two 
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m, 3 






Fig. 1. (Color online) One-dimensional lattices which have a flat lowest-energy one-particle band: 
The sawtooth lattice, the kagome chain, the frustrated diamond chain, the frustrated two-leg 
ladder, the double-tetrahedra chain, and the frustrated (cylindrical) three-leg ladder. The sites are 
enumerated by two indices, where the first one enumerates the cells and the second one enumerates 
the sites within a cell. Alternatively, a one-index numbering of sites is also possible; e.g., for the 
sawtooth chain by m, 1 —»• 2 j, m, 2 —> 2j + 1 etc. and j = 0 ,..., ^ — 1 or for the kagome chain 
by m, 1 — > 3 j, ra, 2 —>■ 3 j + 1, m, 3 —> 3j + 2 etc. and j = 0,..., — 1. While for the kagome 
chain all exchange (hopping) parameters J ( t ) are identical, there are two different bond strengths, 
Ji (£i) and J 2 (£ 2 ), in the other lattices. For the sawtooth chain a special flat-band condition is 
required, namely the exchange (hopping) parameter along the zigzag path J 2 (£ 2 ) has to be 2 
( y/2 ) times larger than the exchange (hopping) parameter Ji (£ 1 ) along the basal straight line. 
For the frustrated diamond chain, the frustrated two-leg ladder, the double-tetrahedra chain, and 
the frustrated three-leg ladder the exchange (hopping) parameter J 2 (£ 2 ) marked by the bold red 
lines must exceed a threshold value to have the flat band as the lowest one (note that for these 
lattices the bold red lines at the same time denote the trapping cells). For the diamond chain we 
also distinguish blue and green lines. For the distorted diamond chain also considered here both 
bond strengths are different, where blue lines indicate exchange (hopping) parameters Ji (t 1 ) and 
green lines exchange (hopping) parameters J 3 (£ 3 ). In case of J\ 7 ^ J 3 (£1 7 ^ £ 3 ) the lowest band 
becomes dispersive. The lowest band is flat (flat-band geometry) if blue and green bonds mark 
identical bond strengths, i.e., J\ = J 3 (£1 = £ 3 ). 


classes of flat-band lattices: One class is given by Mielke’s line graphs (e.g., the 
kagome chain, the kagome lattice or the checkerboard lattice) whereas another one 
can be constructed by Tasaki’s decoration procedure (e.g., the sawtooth chain or the 
two-dimensional Tasaki lattice). Note further, that some of the flat-band lattices be¬ 
long to the class of models with local conservation laws (e.g., the diamond chain, the 
dimer-plaquette chain, the frustrated ladder, and the frustrated bilayer). In Refs.JE: 
and [20] the construction starts from isolated cells which are further arranged in a 
infinite lattice. An important feature of flat-band lattices considered here is that 
the trapping cell is surrounded by triangles, see Figs. [Hand[2] Just this triangular 
geometry of antiferromagnetic bonds prohibits the magnon to escape from the trap 
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Fig. 2 . (Color online) Two-dimensional lattices which have a flat lowest-energy one-particle band: 
The Tasaki lattice, the kagome lattice, the square-kagome lattice, and the frustrated bilayer. The 
sites are enumerated by two indices, where the first vector index enumerates the cells and the 
second one enumerates the sites within a cell. While for the kagome lattice all exchange (hopping) 
parameters J ( t ) are identical, there are two different bond strengths in the other lattices. For the 
Tasaki lattice (two-dimensional counterpart of the sawtooth chain) a special flat-band condition is 
required, namely, the exchange (hopping) parameter along the zigzag path J2 (£2) has to be 4 ( 2 ) 
times larger than the exchange (hopping) parameter J\ (£1) along the basal lines. For the square- 
kagome lattice and the frustrated bilayer the exchange (hopping) parameter J2 (£2) marked by the 
bold red lines must exceed a threshold value to have the flat band as the lowest one (note that for 
these lattices the bold red lines at the same time denote the trapping cells). For the square-kagome 
lattice we also distinguish blue and green lines. For the flat-band geometry blue and green bonds 
mark identical bond strengths. For the distorted square-kagome lattice also considered here both 
bond strengths are different, where blue lines indicate exchange (hopping) parameters Ji (£1) and 
green lines exchange (hopping) parameters J3 (£3). In case of J\ 7^ J3 (£1 / £3) the lowest band 
becomes dispersive. 

due to destructive quantum interference. Obviously, antiferromagnetic couplings on 
a triangle are frustrated, i.e., frustration is a common feature of all lattices discussed 
in this paper. Interestingly, the same mechanism can lead t o a com plete localization 
of electrons in a magnetic field (“Aharonov-Bohm cage”) J201SH One more scheme 
for generating flat-band lattices has been suggested very recently.^ it starts from 
a lattice with dispersive bands, then a set of Fano defects is added, and finally by 
an appropriate transformations of the Fano lattice a corresponding flat-band lattice 
can be obtained, see Sec. mzB It is also interesting to note that flat bands are 


b The presence of two types of states, the localized (dispersionless) ones and the extended (disper¬ 
sive) ones, allo ws to use the universal language of the Fano-resonance phenomenon, known from 
atomic physics see Ref. 1931 


































April 14, 2015 0:38 WSPC/INSTRUCTION FILE de'ri'marevised 


8 O. Derzhko, J. Richter, M. Maksymenko 


often associated with symmetries, see Sec. VII in Ref. [94) 
The localized-magnon state can be written as follows 


|Lm.) =,^a i s i | t •••!)• 


( 2 ) 


Here the sum runs over the sites belonging to one trapping cell C (e.g., the ver¬ 
tical bond for the diamond chain, the frustrated two-leg ladder or bilayer or the 


hexagon/square for the kagome/square-kagome lattice). We denote the energy of 


this state by £i (the flat-band energy). To cage the magnon in the trap C we have 
to fulfill two conditions, namely, XXeC ^n a i ~ 0 f° r VZ ^ C and = constant 

for Mi £ C, see Ref. [8. 

Localized magnons are building blocks for constructing exact many-magnon 
states. Clearly, a state consisting of many localized magnons is an exact many- 
magnon state, if the localized magnons are sufficiently separated of each other, such 
that they do not interact. In case that neighboring trapping cells share a common 
site (e.g., the sawtooth chain or the kagome lattice) the occupation of neighboring 
cells is forbidden (hard-core rule), whereas for lattices with isolated trapping cells 
(e.g., the diamond chain but not the frustrated two-leg ladder or the frustrated 
bilayer) neighboring cells can be occupied (no hard-core rule or more precisely only 
a hard-monomer rule). According to the number of trapping cells Af = aN, where 
the proportionality factor a depends on the lattice, there is a maximum number of 
localized magnons n max oc Af that corresponds to the complete occupation of the 
trapping cells (according to possible hard-core rules). 

It can be shown that the multi-magnon eigenstates built by n localized magnons 
under some conditio ns a re the lowest-energy states in the subspaces with corre¬ 
sponding S z = y — Moreover, it can be shown almost in all cases ri 
that the localized-magnon states are linearly independent (but not necessarily or¬ 
thogonal). 

The multi-magnon state with maximum n = n max is an exceptional state called 
magnon crystal (complete ordered occupation of traps). It resembles a Wigner crys¬ 
tal and it is typically protected by an excitation gap. As a result it corresponds to a 
magnetization plateau state preceding the magnetization jump associated with the 
localized-magnon states, cf. also Sec. 16.21 

An important step to analyze the physics related to localized-magnon states 
consists in the possibility to map the co mplete m anifold of the localized-magnon 
states on classical hard-core lattice gases [10112113 ] p or i ns tance, for the sawtooth 
chain with Af = ^ trapping cells the mapping is onto the hard-dimer problem on a 
simple linear chain of Af sites, and for the frustrated bilayer the corresponding hard¬ 
core problem is that of hard squares on a square lattice (see also Sec. [3j). Using then 
well-elaborated statistical mechanics of hard-core lattice gases one easily calculates 
the low-temperature thermodynamics of the initial quantum spin system, see, e.g., 

Refs, nm Ha m 

In the case of electron systems, the one-particle problem is quite similar. More 
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precisely, the one-electron problem corresponds to one-magnon problem with XY 
interaction in the Hamiltonian ©, i.e., with sfsj + sfsj = ^(s /"sj +si" s p instead 
of Sj ■ Sj in Eq. ©■ However, the construction of many-particle states (i.e., the 
occupation rule of the trapping cells) is different. While for the Heisenberg spin 
system we have to avoid the interactions between trapping cells to remain in the 
ground-state manifold, in the case of the Hubbard system we have to avoid the 
on-site Hubbard repulsion between electrons with different spins in order not to 
run out of the set of ground states. We discuss the one-electron and many-electron 
states of flat-band Hubbard systems in Secs. [8] and El 


2.2. Thermodynamics of flat-band antiferromagnets 

Localized-magnon states constitute a substantial part of all 2 N states of iV-site spin- 
5 flat-band spin systems. Moreover, in the presence of a high magnetic field around 
the saturation value h sa t they are the low-energy ones and in many cases separated 
from the other states by an energy gap. Therefore they dominate the properties of 
the flat-band spin system in the low-tenrperature high-field regime. 

The main observable effects at T = 0 are: A magnetization jump at the satura¬ 
tion field j^l a magnetic-field driven lattice instability,& and a residual entropyP^The 
ground-state magnetization jump between the values M = ^ and M = — n max 

(here n max is the maximal possible number of localized magnons) is in an obvi¬ 
ous manner related to the structure of the energy spectrum of the spin system at 
hand. Namely, the energy of the n-localized-magnon state in the presence of a mag¬ 
netic field, Efm — n£\ — h (^ — n), becomes n-independent at the saturation field 
h sa .t = £i and hence the magnetization S z = — n acquires any value between 

■y and -y — fimaj- A state with n localized magnons has a huge degeneracy, since 
the n localized magnons can be placed over the lattice in many ways. At the sat¬ 
uration field the ground-state degeneracy dramatically increases and corresponds 
to the grand-canonical partition function of an auxiliary classical gas of hard-core 
objects at zero value of the chemical potential. This quantity grows exponentially 
with lattice size N that leads to a residual entropy. 

Moreover, an appropriate lattice distortion fitting to the structure of the local¬ 
ized magnons and thus preserving the localization condition leads to a lowering of 
magnetic energy that is proportional to the displacement of sites. By calculating the 
total energy which consists of the magnetic and elastic parts, one can check whether 
the lattice distortion is favorable or not. It appears that, e.g., the kagome lattice 
exhibits a field-driven lattice deformation (spin-Peierls transition) in the vicinity of 
the saturation field which is accompanied by a hysteresis phenomenon.® 

At nonzero but low temperatures in the vicinity of the saturation field an extra 
low-temperature maximum of the specific heat occurs, indicating that the manifold 
of localized eigenstates sets an additional low-energy scale. Importantly, if the lat¬ 
tice dimension is two (or three) an order-di sorder phase transition occurs at low 
temperatures just below the saturation field. — — The phase transition is related 
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to the ordering of the hard-core objects on an auxiliary lattice, which mimics the 
ordering of the localized magnons. 

Further detai ls can be found in the original papers or review articles on localized- 
magnon systems.®^ As will be shown below some of these features appear also 
in flat-band Hubbard systems demonstrating that the construction and characteri¬ 
zation of many-particle localized states is precisely at the heart of the fact that, to 
some extent, a unified consideration of spin and electron models is possible. How¬ 
ever, there are also important differences owing to the spin degree of freedom of 
electrons, the Pauli principle, and the on-site repulsion which may lead to a co¬ 
herent spin state at T = 0, i.e., to the flat-band ferromagnetism of Mielke and 

TasakiE™ 


3. Emergent Ising Degrees of Freedom 

As it has been discussed already, the manifold of localized-magnon states can be 
completely characterized, i.e., their energies are obtained trivially from the flat- 
band energy eq and their degeneracies can be obtained as a canonical partition 
function of a certain hard-core lattice-gas model. In many cases the energy levels 
above the localized-magnon manifold are unknown. However, in some cases one can 
construct the lowest-energy excited states, too. By including these lowest-energy 
excited states into a lattice-gas approach one can substantially extend the localized- 
magnon description. Instead of a lattice gas with infinite repulsion (hard-core rule) 
we then arrive at a lattice gas with finite repulsion.— 

To be specific, we consider the spiir-^ Heisenberg antiferromagnet m on the 
frustrated two-leg ladder (see Fig. [l]) and the frustrated bilayer (see Fig. [2]). For 
these models one can introduce the total spin on a vertical J 2 -bond t m = s mi i+s mi 2 , 
where m = 1 ,..., AT = ^ runs either over the sites of a simple chain (in case of 
the ladder) or the sites of square lattice (bilayer case). Then the Hamiltonian of the 
models becomes 

+ Ji t m • tz, (3) 

(ml) 

where the second sum runs over the nearest-neighbor bonds of the linear chain or 
the square lattice, respectively. As it is evident from Eq. ©, = t m (t m + 1), 

m = are good quantum numbers and the Hamiltonian of the models 

depends on the set of quantum numbers {t m }, t m = 0,1. As a result, a subset 
of 2^ low-lying eigenstates can be calculated exactly ^ 0 -I and their contribution to 
thermodynamics can be estimated using a classical lattice-gas model. 

Let us consider the relevant eigenstates in more detail. First we consider states 
which consist of n singlets on the vertical bonds mi,... ,m n (i.e., tj n = 0 ) and 
N — n fully polarized triplets (i.e., tjq = 2, t* n = 1) on the remaining vertical bonds. 
If we forbid the occupation of neighboring vertical bonds by singlets (hard-core 
rule) we are faced with independent localized-magnon states and the singlet on a 
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vertical bond is the localized magnon. The energy of the independent localized- 
magnon states is E l ™ = £fm — ne i, s\ = hi = /i sat with £fm = ATJi + Af^-, 
£i = J 2 + 2Ji (ladder) and Efm = 2A/”Ji + £\ = J 2 + 4Ji (bilayer). The de¬ 

generacy of the independent localized-magnon states gj\f(n) is given by the canonical 
partition function of n hard-core objects (hard dimers or hard squares) 

on the lattice (simple chain or square lattice) of A f sites. The independent-localized 
magnons are linearly independent^® and constitute the ground-state manifold for 
S z = N — 1,..., y if J 2 > 2Ji (ladder) or J 2 > 4Ji (bilayer). 

Allowing the localized magnons to be nearest neighbors, we arrive at another 
set of eigenstates, the interacting localized-magnon states. For magnetization values 
S z = M - n, n = 2,..., 4r, the energies of the interacting localized-magnon excited 
states are: E” = E l ™ + vJ\, where v is the number of pairs of neighboring localized 
magnons. The interacting localized-magnon states are the low-lying excited states 
for these values of S z if J 2 is large enough, i.e., J 2 > J 2 (strong-coupling regime). 
From exact-diagonalization data for finite systems we have found for J 2 above which 
the strong-coupling regime holds the values J 2 « 3.00Ji (ladder) and J 2 « 4.65Ji 
(bilayer). Obviously, for the models at hand the localized-magnon ground states are 
separated by a gap from excitations in the strong-coupling regime, i.e., for J 2 > J 2 . 
Note, however, that the existence of an excitation gap is not necessarily a common 
feature of localized-magnon systems. 

For lower values of the magnetization, S z = 4(- — r, r = 1,..., 4r, where no 
independent localized-magnon states exist, the class of interacting localized-magnon 
states contains the ground-state manifold as well as low-lying excited states in the 
strong-coupling regime. The ground-state manifold is built by n = + r (r = 

1 ,..., 4r) localized magnons, where, e.g., magnons occupy one sublattice of the 
underlying lattice (simple chain or the square lattice) completely, and the remaining 
r localized magnons sit on the other sublattice. The energy of this state is Eu +r = 
—— r J 2 , J 2 = h 2 . The low-lying excited states are constructed from the ground 
state by rearranging the localized magnons to increase the number of neighboring 
magnons. Then each new pair of neighboring localized magnons increases the energy 
by Ji. 

Note that the interacting localized-magnon states can be visualized as par¬ 
tially overlapping hard-core objects in contrast to the independent localized-magnon 
states which can be visualized as non-overlapping hard-core objects. We are inter¬ 
ested in the partition function Z(T,h,N) of the frustrated quantum spin system 
which, after taking into account that in the strong-coupling regime the independent 
and interacting localized-magnon states dominate at low temperatures and high 
fields, can be written within the frames of a lattice-gas model of classical particles 
with finite nearest-neighbor repulsion V = J\ 


Z{T, h, N) ~ £ ■■■ E e ' 

ni—0,1 n at—0,1 


%M —hAS +(h — hi ) E m n m+Jl 12(ml) n m n l 


= e- —Z lg (T,g,M), (4) 
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Fig. 3. (Color online) Phase diagrams as they follow from the lattice-gas description of the s = ^ 
Heisenberg antiferromagnet in a magnetic field on the frustrated two-leg ladder with J\ = 1, 
J 2 = 4 (left) and on the frustrated bilayer with J\ = 1, J 2 = 5 (right). A staggered occupation 
of vertical bonds by localized magnons (a kind of “antiferromagnetic” long-range order) occurs 
along the red line between fi 2 — J 2 — 4 and h± = J 2 + 2 J\ = 6 at T = 0 (ladder) or below the red 
critical line T c (h) with starting point /i 2 = J 2 = 5, T = 0 and endpoint hi = J 2 + 4Ji = 9, T = 0 
(bilayer), whereas a uniform occupation of vertical bonds (a kind of “ferromagnetic” long-range 
order) occurs along the green lines h < /i 2 at T = 0 and h\ < h at T = 0. The remaining part of 
the phase diagrams corresponds to a disordered phase. Moreover, we show the lines below which 
the exact-diagonalization data and the lattice-gas predictions for the specific heat of the finite 
system of N = 16 sites coincide with the accuracy up to 5%, 2%, and 1% (blue dashed lines from 
top to bottom). 


where Si g (T, /i, V) = X] ni=0 i • ■ • S nj v=o 1 e ^ ~ is the grand-canonical parti¬ 
tion function of the classical lattice-gas model with the Hamiltonian 'H({n m }) = 
n m + j\ Y^(mi) n m,ni, fJ, = hi — h, and hi = £1 is the saturation field h sat . 
Introducing the Ising variables a m = 2n m —l (i.e., n m = l+ ° rn ) one arrives at a clas¬ 
sical antiferromagnetic Ising model in a uniform magnetic field. Extensive studies of 
the relation between the frustrated quantum spin models and the classical lattice- 
gas models for finite systems (up to N = 32) show that the effective description 
[i.e., the (non-frustrated) Ising model] can reproduce the thermodynamic properties 
of the initial frustrated quantum spin model in the strong-coupling regime!^ 

Since the emergent classical models may be easily studied either analytically 
(ladder) or using classical Monte Carlo simulations (bilayer), a lot of results for 
these models are available (see, e.g., Refs. f97lf98lb Hence we can discuss the proper¬ 
ties of the frustrated quantum antiferromagnets at hand on the basis of this knowl¬ 
edge. In Figs. [3] and 0] we report the phase diagram of both models in the half 
plane “field-temperature”, the low-temperature magnetization curves and the tem¬ 
perature dependence of the specific heat below the saturation field for the bilayer. 
An interesting low-temperature behavior in the magnetic-field region /12 < h < hi 
is related to an ordering of localized magnons. Indeed, both auxiliary lattices are 
bipartite ones and in the two-dimensional case (i.e., for the square lattice) the oc¬ 
cupation of the two sublattices exhibits a finite-temperature order-disorder phase 
transition corresponding to that one which is present in the square-lattice Ising 
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Fig. 4. (Color online) Zero-temperature magnetization curves for the ladder (Ji = 1, J 2 = 4) 
and bilayer (Ji = 1, J 2 = 5) models (left) and specific heat as a function of temperature for the 
bilayer model for various magnetic fields below the saturation field (right). The plateau at 41 = 1 
is between /12 = J 2 and the saturation field h\. The specific-heat data were obtained by Monte 
Carlo simulations of the classical lattice-gas model with finite repulsion V = J\ for systems of up 
to 800 X 800 sites exploiting usual Metropolis algorithm with up to 3 000 000 Monte Carlo steps. 
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antiferromagnet in a uniform magnetic field!^ Crossing the critical line T c (h) the 
thermodynamic quantities for the frustrated bilayer exhibit singularities: The spe¬ 
cific heat shows a logarithmic singularity, the staggered magnetization (order pa¬ 
rameter) within the (“antiferromagnetically”) ordered phase decays to zero with the 
critical exponent j3 = |, and the staggered susceptibility diverges with the exponent 
7 =7. These field and temperature dependences of various thermodynamic quanti¬ 
ties clearly indicate the ordering of the localized magnons on vertical bonds of the 
frustrated bilayer which may occur for T < Xbnsager = 2111 (^ 2 + 1 ) ^ 1 ~ 0.567296Ji. 
The investigations presented in Ref. [3T| can be obviously extended to other two- 
dimensional lattices, for example, frustrated bilayers consisting of two triangular, 
honeycomb, kagome etc. layers. Furthermore, these findings may be of relevance to 
understand the phy sical properties of the recently synthesized spin dimer magnet 
E^CoSbOgC^!^ Although, this compound has a pronounced XY anisotropy of 
the antiferromagnetic Heisenberg interactions and certainly also a deviation from 
ideal flat-band geometry remnant features of the discussed bilayer system could be 
present, cf. also Sec. 17.21 


4. Chiral Localized Magnons 

Most of the investigations of localized-magnon systems consider the case of a bi¬ 
partite trapping cell, e.g., a single bond, a square or a hexagon. Such cells have 
a non-degenerate one-magnon ground state which remains the ground state of the 
whole lattice after suitable connection of the cells into the infinite lattice. However, 
one can also consider a non-bipartite trapping cell, e.g., a triangular trap. The trian- 
gular trap has a twofold degenerate ground state that leads to new effectsjSBEHESET] 
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This degeneracy is associated to the chirality. The chirality operator for the triangle 
may be introduced as follows: 


X = ^=(si • [s 2 x s 3 ]) 


(5) 


(see the corresponding triangular cells in Fig. [[]). As a result, the lowest-energy 
magnon band having the energy £i i2 (k) = — §J 2 — J\ = —£ 1,2 (double-tetrahedra 
chain, see Fig. © or £i, 2 (k) = — |J 2 — 3Ji = —£i )2 (frustrated triangular tube, 
see Fig. |T|) is twofold degenerate, i.e., the one-magnon flat-band states are 2A/ r -fold 
degenerate. The flat-band states are located within triangles and can be written 
explicitly as follows: 


|T)m — (s m ,l + UJS m ,2 + Uj2s m, 3 ) I t - - ■ t)> 

I ) m = ( S m,l + U}2s m ,2 + UJS m, 3 ) I t • • ■ t), (6) 

where w = e^~ l . Using Eqs. (0 and 0 one can check that Xm|±)m = ±|±) m . 

Moreover, the flat-band states become the lowest-energy ones if J 2 > 2 J 3 . 

We pass to the low-temperature high-field thermodynamics which is conditioned 
by the set of independent localized-magnon states. The formula for the (dominating) 
contribution of the localized-magnon states to the partition function reads 


Z(T,h,N) 


n max 

E 

n —0 


m{n)e~ 


s-h-n(hi-h) 


(7) 


where hi = h SSb t = £i j2 is the saturation field and gj\[(n) is the number of localized 
states of n magnons on a lattice of N cells. For gjg(n) now we have to take into 
account also the chiral degrees of freedom. That yields gjg{n) = 2 n C) J- for the double- 
tetrahedra chain and gjg(n) = 2 n Z\ 1( i(n,Af) for the frustrated triangular tube, where 
C^f = ('^') is the binomial coefficient and Zhd(n.,A0 denotes the canonical partition 
function of n hard dimers on a A/”-site linear chain. Here the factor 2 n stems from the 
chirality. The sum in Eq. © can be easily calculatecPSland, as a result, one gets the 
thermodynamics of the models at hand in the low-temperature high-field regime. 
The considered spin systems exhibit the typical features of localized-magnon sys¬ 
tems: The magnetization jump, the residual entropy, and the extra low-temperature 
peak in the specific heat. Furthermore, similar to the case of the frustrated two-leg 
ladder (see the previous section), one can elaborate for the frustrated three-leg lad¬ 
der a soft lattice-gas model instead of the hard-dimer model to take into account 
the low-energy excitations!® Based on finite-size calculations the critical value of 
J 2 above which the strong-coupling regime holds was estimated as ~ 2.68Ji. 
Note that in contrast to the two-leg ladder, the lattice-gas model for the three-leg 
ladder does not possess the symmetry around h\ and h 2 , because of the chirality. 

The degeneracy owing to the chirality can be lifted, e.g., by a Dzyaloshinskii- 
Moriya interaction. If we add a small perturbation in form of a Dzyaloshinskii- 
Moriya interaction D = De z = (0,0,14) between neighboring spins of each 
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Fig. 5. (Color online) Low-temperature dependence of the specific heat for the double-tetrahedra 
chain (left) and frustrated triangular tube (right) slightly below the saturation field h = 0.98/ii 
(J\f = 4, J\ = 1, J 2 = 5, D = 0, 0.025, 0.05, 0.1). Exact-diagonalization data (symbols) are 
compared to analytical predictions (lines). 

triangular trap, the partition function is again given by Eq. ©, however, 
with gjg(n) = (2 cosh ^p-) n C^ for the double-tetrahedra chain and gjg(n) = 
(2 coshA/”) for the frustrated triangular tube. The introduced per¬ 
turbation leads to a splitting of the ground-state levels and, therefore, one more 
low-energy scale arises. As a result, the low-temperature features close to the satu¬ 
ration field /i sa t = h\ are more subtle. For example, the temperature profiles of the 
specific heat show more tiny features which can be seen in Fig. [5] For a special set 
of parameters the temperature dependence of the specific heat may exhibit three 
(four) maxima for the double-tetrahedra chain (frustrated triangular tube). The 
low-temperature maxima are excellently described within the presented effective 
low-energy theory. 

The considerations of this section refer to the antiferromagnetic Heisenberg 
model. Later on in Sec. l8.llwe will discuss how they can be straightforwardly applied 
to the Hubbard model.^SSS ^6 


5. Beyond the Ideal Flat-Band Geometry 

The localized-magnon picture emerges only in the case of strict localization of 
magnons or, in other words, when the lowest magnon band is strictly flat. Of course, 
this requires certain relations between exchange couplings which hardly can be sat¬ 
isfied in a real-life case. Among the reasons for that one may mention nonideal 
geometry or some extra exchange couplings which lead to a violation of the local¬ 
ization condition. However, one may expect to find a magnetic compound which is 
quite close to the localized-magnon prototype. Therefore it is important to consider 
spin models with a small deviation from the strictly dispersionless case, i.e., the 
former flat band acquires a slight dispersion, and to figure out how the localized- 
magnon features are modified by (weak) dispersion. 

In what follows we discuss the effect of small deviations from the condi¬ 
tions for localization con side ring three models which belong to the so-called 
monomer universality classnamely, the diamond chai n (see Fig. [I]), the dimer- 
plaquette chain, and the square-kagome lattice (see Fig. [2]). ^ HU For large exchange 
bonds on the trapping cell, J 2 , an appropriate approach is the strong-coupling 
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approximatioii' 7 ° 99. which, in fact, was elaborated for a distorted frustrated dia¬ 
mond spin chain in Refs. [3T] and [1001 Let us discuss this approach in some detail 
considering as an example a distorted frustrated diamond chain, see Fig. [Q Inspired 
by the azurite geometry, we assume that J\ ^ J 3 and set Ji + J 3 = 2 J < J 2 (for 
a more general nonideal geometry see Ref. HU). The starting point of the strong¬ 
coupling approximation is the simple problem of two spins interacting via a J 2 bond. 
In a high magnetic held the state |u) = | tit 2 ) has the lowest energy ^ — h whereas 
just below the “bare” saturation Held ho = J 2 the state | d) = ^=(| tils) — | Iit 2 )) 
(localized magnon) becomes the lowest-energy one with the energy — The col¬ 
lection of spin dimers on ,/ 2 -bonds and the isolated spins at the site m, 3, see Fig. [TJ 
at h = ho constitutes the “main” Hamiltonian F ma ; n whereas the remaining J\ and 
J 3 interactions represent the perturbation V = H — H ma j n . We introduce the model 
space which is spanned by the ground states of the main Hamiltonian |y>o) and the 
projector onto this space P = |v ? o){v ? o| 


P — ^m^m] 


Pm=Vm ® (| l3><t3 I ) m , ^ = (l«)H + \d)(d\) r 


( 8 ) 


We are interested in an effective Hamiltonian H e g- which acts in the model space only 
but gives the ground-state energy of the Hamiltonian H. The effective Hamiltonian 
can be found within perturbation theorjff2mn2II03| 


F eff 


PHP + PV^2 


W a ){y a \ vp 

£■0 £a 


(9) 


where \<p a ), o. ^ 0, denotes the excited states of F ma ; n . Moreover, it is convenient 
to introduce (pseudo)spin-f operators 

T z = ±(\u)(u\-\d)(d\), T+ = \u)(d\, T~=\d)(u\ (10) 

for each cell m. As a result, one arrives at the spin-i isotropic XY model in a 
^-aligned Held with the Hamiltonian 

H eS = MC - h Y, r* + j E ( T - T - + T - T ") > (H) 

m (ran) 


where the last sum runs over all nearest-neighbor bonds on a simple chain. The pa¬ 
rameters of the effective Hamiltonian are given by h = h—hi— ^^j^ , J = , 

and hi = J 2 + J, J = \ — , see Refs. 11001 [5l l :4TJ For the dimer-plaquette chain 
and the square-kagome lattice the effective model is again given by the Hamiltonian 
in Eq. CD, however, with specific values of parameters, see Ref. 1411 

In Fig.[G]we illustrate the quality of the strong-coupling predictions (long-dashed 
green lines). For this purpose we compare strong-coupling-approach predictions with 
exact-diagonalization data for the initial model. The comparison of data for the 
magnetization and the specific heat reported in Fig.[ 6 ]shows that the strong-coupling 
theory gives correct qualitative predictions but it does not provide quantitative 
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Fig. 6 . (Color online) Magnetization curves (left) and low-temperature dependence of the specific 
heat (right) for a distorted frustrated diamond chain (J\f = 6 , Ji = 0.85, J 2 = 3, J 3 = 1.15). 
Exact-diagonalization data (black bold solid curves) versus effective-model predictions, namely, 
strong-coupling approach (1111) (green long-dashed curves) and localized-magnon approach J12D 
(blue short-dashed curves) and m (red dotted curves). 


details. That is not astonishing after all, since we start from the limit of small Ji 
and J 3 in comparison to J 2 and consider only the second-order term in Eq. ©• 
The strong-coupling calculation does not use the special features of localized- 
magnon picture, where we assume Ji = </ 3 = J < 4r. Starting from the localized- 
magnon description we need only | J 3 — Ji| to be small. Hence, one may expect that 
such an approach would contain more relevant physics than the strong-coupling one 
already in the “main” (i.e., unperturbed) problem. To follow this idea we may apply 
Eq. © for the following splitting of the Hamiltonian H: H ma ; n corresponds to the 
ideal geometry (i.e., with J = ' 7l j ,' 73 instead of J, and J 3 at h = hi with hi denoting 
the saturation field for the ideal-geometry case) and V = H — H ma i n . However, we 
are immediately faced with a difficulty in calculating the second term in the r.h.s. 
of Eq. ©, since only the ground states of H ma \ n (i.e., the set of localized-magnon 
states) are known, whereas the excited states | ip a ), a ^ 0, remain unknown. To 
overcome this obstacle we start from the Hamiltonian 


U = VHV, V = ® m.Vm- 


( 12 ) 


Here the projector V m [see Eq. ©] restricts the states of the cell m to only two 
states | u) and |d). After this approximation we can split the Hamiltonian U into 
the main part K ma i n which corresponds to the ideal geometry and the perturbation 
V = U — "Hmain- Now all eigenstates of "Hmam are known and we can easily calculate 
the second term in the r.h.s. of Eq. ©• As a result, for the considered frustrated 
quantum antiferromagnets we arrive at the spin-i XXZ Heisenberg model 


u cS = ATC - h ]T r* + )T [J (T*T* + TO) + J Z TO] ■ 

m ( mn ) 


(13) 
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The parameters of Heft for the distorted diamond chain are given by 

C = - fr -^- ( J 3- Jl > 2 f * + 1 ^ 
4 2 16j2 / 

4 2 16J 2 \ J 2 ) 

17 7 (J3 — J1) 2 , T , 7 7 "^3 + ^1 

h = h — hi -—-, hi = J 2 + J, J = 


4J 2 


2 ’ 


, (J 3 - Jr) 2 1 (J 3 - J 1) 2 (, , J , 

4J 2 l--£ 4J 2 V J 2 

j- = A_A! 1 _ d = (L + ...,. 


4 J 2 V J 2 


(14) 


To compare with the strong-coupling approach we give in Eq. (1141) except the full 
expressions also the leading terms of an expansion in . Obviously, the lowest-order 
terms agree with the strong-coupling approximation m Note that the next-order 
terms in Eq. (1141) correspond to the next-order terms within the strong-coupling 
scheme obtained in Ref. llOOt The effective Hamiltonian parameters for the dimer- 
plaquette chain and the square-kagome lattice (even for a more general nonideal 
geometry) can be found in Ref. UK From Fig. [G] one can clearly see that the 
perturbation-theory approach based on the localized-magnon picture yields much 
better results than the strong-coupling approach, and it can provide even a quan¬ 
titative description of the initial model. We mention, that the ideal-geometry limit 
of the effective model gives an alternative (but identical) description of the thermo¬ 
dynamics previously discussed, e.g., in the reviews [5511351 . 

We conclude that the nonideal geometry leads to new effective (unfrustrated) 
quantum spin models which are well known in statistical mechanics. Hence, we may 
use the broad knowledge on their properties to understand the low-temperature 
high-field behavior of the frustrated quantum antiferromagnets at hand. The most 
interesting model is the two-dimensional square-kagome lattice, since the ef fectiv e 
models exhibit the famous Berezinskii-Kosterlitz-Thouless (BKT) transition!^ 4 

Let us briefly outline the relevant results concerning the BKT transition. Since 
the early 1970ies it is known that the classical two-dimensional isotropic XY model 
undergoes a transition from bound vortex-antivortex pairs at low temp eratu res to 
unpaired vortices and antivortices at some critical temperature Tbkt-^ 4 Below 
Tbkt (superfluid phase) the system is characterized by quasi-long-range order, i.e., 
correlations decay algebraically at large distances without emerging of a nonvan¬ 
ishing order parameter. Above Tbkt (normal phase) the system is disordered with 

b 

an exponential increase of the correlation length £ as T —> Tbkt + 0, £ oc 
with t = T x^kt T ' 4416 cr itical temperature for the classical square-lattice isotropic 
XY model (without magnetic field) is Tbkt ~ 0.893|J|. 10'^ In the quantum spin- 
^ case the BKT critical behavior also occurs, however, the critical temperature 
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is lower, Tbkt ~ 0.34|J|1^® The quantum model is gapless with an excitation 
spectrum which is linear in the momentum. The specific heat c(T) shows T 2 be¬ 
havior as T —> 0, increases very rapidly around Tbkt, and exhibits a finite peak at 
temperature T* somewhat above Tbkt- This kind of the low-temperature thermo¬ 
dynamics survives for not too large ^-aligned magnetic field |h| < 2|J|1^S1 Also for 
spin-i square-lattice XXZ model with dominating isotropic XY interaction in a 
^-aligned magnetic field the BKT transition appears, 

Now we transfer this knowledge to the distorted square-kagome Heisenberg an- 
tiferromagnet. We use the position of the maximum T* in the specific heat as an 
indicator of the BKT transition to construct a sketch of the phase diagram of the 
model, see Fig. 0 From the temperature dependence of the specific heat (see left 
panel of Fig. 0 ) calculated by exact diagonalizations and quantum Monte Carlo sim¬ 
ulations for the effective model (fl3l) one observes an agreement with the expected 
T 2 behavior for T —> 0 in a small region below the saturation field (and an expo¬ 
nential decay outside this region). For gapped systems there is a broad maximum in 
the specific heat. In the gapless field region the maximum occurs at lower temper¬ 
atures and it is more pronounced becoming peak-like. Furthermore, in the gapless 
field region the temperature profiles show noticeable finite-size effects. A tentative 
BKT-transition line TbktW (green short-dashed curve in the right panel of Fig.0 
drawn on the basis of the specific heat data provides a generic phase diagram of the 
distorted square-kagome Heisenberg antiferromagnet at low temperatures and high 
fields. The largest transition temperature appears for zero effective field h, which 
corresponds to h = hi + 1 _ 1 j (hi = 2 J 2 + J, J = 73 (/ 7l ) of the initial 

frustrated model. A deviation from this field yields a decrease of Tbkt, and finally 
Tbkt becomes zero. Thus, in the highly frustrated quantum Heisenberg antiferro¬ 
magnet on the square-kagome lattice with deviations from ideal flat-band geometry 
the BKT transition may appear in a certain region of finite magnetic-field values, 
i.e., the model exhibits a magnetic-field-driven BKT transition. 

Interestingly, recent experimental and theoretical studies of the complex spin- 
dimer material C 36 H 48 C 112 F 6 N 8 O 12 S 2 (having a spin interaction network corre¬ 
sponding to a two-dimensional distorted ho neyc omb lattice) have given evidence of 
a field-induced BKT physics in this systemff® 


6. Closely Related Studies on Spin Systems 
6.1. Bose condensation in flat bands 

One interesting study related to localized-magnon physics that has been recently 
presented by S. D. Huber and E. Altman in Ref. l42l deals with Bose condensation 
in flat bands. The authors of Ref. [42] consider a Hamiltonian of the form 

H = E IM (&& + H e.) + | 

(u) * 


(15) 
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Fig. 7. (Color online) Distorted square-kagome lattice with J\ — 0.8, J 2 = 2, J 3 = 1.2. Left 
panel: Low-temperature specific heat per site c(T, h ) at high fields obtained by quantum Monte 
Carlo simulations for the effective Hamiltonian d. The line thickness increases with increase of 
A f ■ For comparison, we also show exact-diagonalization data for V = 20 by very thin curves with 
circles. Right panel: Sketch of the phase diagram in the “field-temperature” half plane. The sketch 
of the phase boundary (thick dashed green line) is obtained from the position of the maximum T* 
in c(T, h) that is shown by the thick solid blue curve (A f = 20) and the blue symbols (triangles: 
J\[ = 100, squares: J\[ = 256, pentagons: J\[ = 400, circles: A f = 576). The endpoints of phase 
boundary at h m 4.996 and h ps 5.027 are indicated by filled violet circles. 
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where b\ (pi) is a bosonic creation (annihilation) operator on lattice site i. Usu¬ 
ally weakly interacting bosons form a condensate in the lowest eigenstate of the 
kinetic-energy operator. However, if the hopping parameters |1 provide a flat low¬ 
est band, the nature of the ground state will be determined by the interaction U 
acting within the hugely degenerate flat-band-states manifold. S. D. Huber and 
E. Altman co nsid er the one-dimensional sawtooth lattice and the two-dimensional 
kagome lattice.S^Both lattices have a lowest flat band in the one-particle spectrum. 
Such bosoni c mo dels with flat bands are of direct relevance for ultracold atoms in_op- 
tical latticeJl®] or diluted frustrated spin -1 magnets with fully polarized spins LliOl 
Exact many-body ground states can be constructed by occupying some subsets 
of spatially non-overlapping trapping cells, i.e., at sufficiently low lattice fillings 
(densities of bosons) v. Again, the multiple possibilities to put the bosons into the 
non-overlapping localized states lead to a massive degeneracy of the ground state. 
There is a critical density, ry^x = n ™ x = \ {N is the number of lattice sites) for the 
sawtooth chain and iz max = 5 for the kagome lattice at which a close-packed crystal 
state is formed. So far, there is a close correspondence to the spin systems discussed 
in the preceding sections. S. D. Huber and E. Altman now consider lattice fillings 
slightly above close-packed filling iz max . The wave function of any additional boson 
must overlap spatially with those of other bosons and then an exact construction of 
eigenstates is not possible. The question addressed in Ref. [42] is: What is the struc¬ 
ture of the ground states at lattice fillings slightly above v max l In particular, do 
the added bosons condense to form a superfluid on top of the close-packed crystal 
state? We will briefly discuss this question now. The strate gy to attack the problem 
follows a previous study of the quantum Hall effect The authors project the 
Hamiltonian m onto the flat band and obtain this way an effective low-energy 
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Hamiltonian which depends only on the (weak) on-site interaction U. The resulting 
low-energy model is defined on a new lattice, it is free of frustration and can be 
analyzed further using standard methods. 

We start with the sawtooth-chain model. On the one-particle level the Hamilto¬ 
nian m corresponds to the spin-i isotropic XY Hamiltonian, since the commu¬ 
tation relations are not important. One easily finds that the localized states exist 
if t 2 = y/21\. They have the form \Vj) = — ^(b] 4 J _ 1 — \f2\2g j + b\ j)|vac) and their 
energy is —2t\ = — 2t. To achieve a better comparison with the original paper of 
S. D. Huber and E. Altman we use here the indices of Ref. 02J These indices A, j 
and B,j correspond to to, 2 and to, 1, respectively, in Fig. [T] and the lattice has 
27V" sites. On the many-particle level the particle statistics becomes important. One 
can construct many-boson ground states of the interacting Hamiltonian by occupy¬ 
ing only non-overlapping | Vj) states. In particular, the crystal state [charge-density 
wave (CDW) state, all non-overlapping trapping cells are occupied] is given by 
— t 

icdwi) = n/=i 14 b|vac). The states | Vj) are linearly independent and complete 
in the flat-band subspace, however, they are not orthogonal. It is convenient to 
introduce an alternative set of localized states within the flat-band subspace which 
form an orthonormal basis 


w} |vac) = [ W A (. r j - r i) b A,j + W %( r j - r i) b B,j l vaC )> 


w A (rj ) = ^~ T d K cos%e^-^, Wg ( r .) = -±- f dn sin%e^' ( 16 ) 


27T 


27T 


[for explicit expression for d K see Eq. (4) in Ref. 02]. Although these Wannier states 
(or W-states) locally have the structure of the E-states, the W-states decay ex¬ 
ponentially with a localization length of ^ « log(2.15)ap^ where a is the lattice 
constant. 

Now additional particles are introduced into the system and a small interaction 
U < A is assumed (for the opposite limit U > A, see Ref. I112p . where A = 2f is 
the energy gap from the flat band to the next band. To implement the projection 
onto the flat band, we express the bosonic operators in terms of the Wannier ba¬ 
sis, b \(B) j = 12m w A(B)( r J — r m)Wm + higher bands, and neglect the contribution 
from the higher bands. The resulting Hamiltonian contains a renormalized on-site 
interaction U « 0.40C/, further-range interactions as well as assisted-hopping terms, 
where all effective parameters are proportional to U. Since the renormalized on-site 
interaction U is the largest parameter, S. D. Huber and E. Altman replace the 
on-site repulsion by a hard-core constraint and arrive at a spin-^ chain model 


w- 

II 

to 

rdw 

J z rr z (j z + - _ 

A + T\ 

K4l^±2+ H - C -) 

1 a i a i+1 ' 2 

V • d 


+/ sp 1 


( a t-l a i+l + H - C 0 


(17) 
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with the set of parameters 

U « 0.40f7(—>■ oo), I z « 0.11217, J dw « -0.02517, 7 sp « -0.01117. (18) 


We turn to the kagome-lattice model. Diagonalizing the tight-binding 
part of Hamiltonian we get three bands, £o(k) = —2 1, e±(k) = f{l ± 

a /3 + 2[cos k\ + cos(fci — fo) + cosfo]}. Localized states are located within single 
hexagons. Importantly, there is a band touching at the T point, i.e., the flat band 
is not gapped. As a result, the Wannier states are not exponentially localized but 
have power-law tails (the slowest decay is ~ jiy). Projecting onto the flat band we 
are faced with a further complication, because of the vanishing gap in the one-boson 
spectrum between the flat band and the second band: In the case of the kagome 
lattice we do not have the small parameter -k-. However, the authors of Ref. 35] ar¬ 
gue that for fillings v > | the interactions open an effective band gap which finally 
controls the projection onto the flat band. Again a hard-core constraint is supposed, 
and, finally an effective Hamiltonian in terms of spin-i operators on the triangular 
lattice is derived: 


H eff ~ X! ib { y«+a Q + (< + \ 

i a =1 f ' 


+I X 2 V < + a Q (c 


ri+a a +a a+ i 


y r,+a a +a a _i 

h “n+a, 


/ f !/ <+a^r 1 +a a+2 


I 3 V < J t i+aa ^r i+aa+3 


<T ri+ 2a a + H.C.] } + .. . 


(19) 


with the set of parameters 

jxy 

U = 0.1417 (-> oo), I z ^ 0.028 U, -j- 


1 I_^_ _ 1 If_ 
3 ’ l 7 ~ 6 ’ l 7 


i i_z 

8’ I z 



Here the vectors a a , a = 1, ..., 6 , connect the nearest neighbors on the triangular 
lattice, see Ref. |42] 

The effective spin-Z Hamiltonians m (one-dimensional case) and (1191) (two- 
dimensional case) are much easier to study compared to the initial models. In the 
one-dimensional case a bosonization technique can be applied whereas in the two- 
dimensional case a mean-field treatment can be elaborated. The results of the inves¬ 
tigation of the fate of the CDW state at lattice filling above t'max are different for the 
sawtooth chain and the kagome lattice. For the sawtooth chain the long-range or¬ 
der of the CDW state will be immediately destroyed by proliferation of free domain 
walls upon increasing the density of bosons, i.e., the exact ground state |CDWi) 
melts due to delocalization of domain walls. The universal low-energy properties 
at filling v = ^ + 5u are described by a commensurate-incommensurate transition. 
For the kagome lattice the doping of bosons into interstitial sites of the CDW is 
energetically favorable versus adding bosons as domain walls. Since the interstitial 
bosons do not destroy the CDW, the long-range order of the CDW state is expected 
to survive and to coexist with a condensate formed by mobile interstitial bosons, 
i.e., one deals with a supersolid phase. In Ref. [42]one can find further details con¬ 
cerning the sawtooth chain and the kagome lattice above t'max obtained on the 
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basis of effective Hamiltonians 03 and m, as well as discussions about possible 
realizations of the kagome system using ultracold atoms and frustrated quantum 
magnets in high magnetic fields. A similar study for the frustrated two-leg ladder 
is reported in Ref. 11131 (see also Ref. I114D . 

6.2. Valence-bond plateau states of the antiferromagnetic kagome 
spin-lattice 

The spin-1 Heisenberg antiferromagnet on the kagome lattice is one of the most 
interesting and challenging problems in highly frustrated magnetism. Although nu¬ 
merous studies of the kagome-lattice antiferromagnet exist, the physics of the model 
in zero field is still not fully understood, see, e.g., Refs. HI 151 1751 1751 The application 
of a magnetic field probes the low-energy states of subspaces with various values 
of S z and leads to new intriguing features, see Fig. [HJ Thus, the model exhibits 
a well-pronounced plateau in the (ground-state) magnetization curve at of the 
saturation magnetization.^ 27 Il°-Il61ll7] N 0 te 7 however, that quite recently the 
existence of this ^-plateau was questioned .® 0 It is noteworthy that recent experi¬ 
ments on two kagome compounding have found indications for a plateau close to 
(but not exactly at) magnetization m = 1 . 

In Ref. Ill7i it was argued that the states of the ^-plateau in the Heisenberg 
model and the corresponding one in the Ising limit belong to the same phase. 
In particular, a quantum valence-bond-crystal state was proposed as a candidate 
for the 1-plateau ground state. On the other hand, it is known that there is a |- 
plateau preceding the jump to saturation due to localized-magnon states.^^The |- 
plateau ground state can be rigorou sly constructed: It is a threefold degenerate close- 
packed magnon-crystal statcPSEi] that is also of quantum valence-bond nature. Its 
structure is sketched in Fig. [8] (left panel): In the background of polarized (“up”) 
spins, on every third hexagon marked by a dashed circle one flipped (“down”) spin 
is distributed over the sites of the hexagon. The width of the 1-plateau is about 
0.1 j}H| 4 3 | 44 | F ur thermore, the numerical data for the magnetization curve shown 
in Fig. [ 8 ] (right panel) give evidence for another plateau at m = |. Having the 
valence-bond picture of the m = | = 5 and m = | plateau states in mind, in 
Ref. [44] a unified description of the magnetization-plateau states at m = |, | and 
1 was proposed, that also provides a clarification of earlier controversial proposals 
concerning kagome plateaus. ^ ” G! 

For the unified description of plateau states the exact wave function for the 
m = | plateau is generalized to corresponding valence-bond-crystal wave functions 
for the candidate plateaus at m = ^ and |. As mentioned above the independent 
localized-magnon state for m = | is the so-called magnon-crystal state (see Fig. [ 8 ]) 

i 4 Bc>=ni^iiit><. (21) 

3 1 

Here the first product runs over an ordered pattern of all non-overlapping trapping 
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Fig. 8. (Color online) Magnon-crystal state (left) and magnetization curve (right) for the quantum 
antiferromagnet on various finite kagome lattices. 


cells (hexagons) denoted by the dashed circles and the second product runs over 
the remaining sites. The localized-magnon state of a hexagon is 


\L,i) = lUttttt))*- 


( 22 ) 


We have introduced here the momentum eigenstates for a hexagon 11cr 0 ... < 75 ))^ = 
]Cr=o e lfcr |oy ... 05 +r}, where a n =t,4-, n + r has to be read modulo 6 , and AT is 
a normalization factor ensuring k{(&o ■ ■ ■ & 5 \\cro ■ ■ ■ a 5 ))k = 1 (A f = 6 unless the state 
repeats under less than 6 translations). The magnon-crystal state is the threefold 

degenerate ground state in the subspace with S z = A- and its energy per site at 
1 . 

h = 0 is ey BC = i- 

Note that the quantum valence-bond crystal state proposed in Ref. 11171 for the 
4-plateau is described by a similar wave function: The global pattern is again as 
sketched in Fig. [HI but the dashed circles were supposed to represent a combination 
of the two Neel states of a hexagon. This picture works well approaching the Ising 
limit of the model, but for the Heisenberg model it is favorable. I to consider the 
true ground state of the Heisenberg model on the hexagon, |L,|).|), instead of 
the combination of the two Neel states. The corresponding threefold-degenerate 
valence-bond-crystal state reads 

I^vbc) = ni^m t>* (23) 

3 l 


with 


\L,Ui) = 

3(3 


+ 


\/l95 
\/l3) 


3||IUttt», 


5l/[3 L 

iitut)), -11 unit)).) 


V3(4+Vl3)||ntm))J- (24) 


Although this state is not an exact eigenstate of the Hamiltonian it is a good 


representation of the true ground state in the subspace with S z = IA energy 


per site at h = 0 is ey BC = — — 


VTsj 

18 


-0.3114195153J. 
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J’/J 


p.724 

p.616 


Fig. 9. (Color online) Exact-diagonalization data for the overlaps (^Sactl^VBc) °f the va l ence_ 
bond-crystal states with the exact ground states for m = ^ and m = •§ as function of ratio of 
bond strengths , where J is bond strength of the trapping cells hosting localized magnons in 
the magnon-crystal state and J' is bond strength of all other bonds. The arrows at the right side 
indicate the overlap of a symmetric linear combination of the three degenerate valence-bond-crystal 
states with the true ground state for the uniform kagome lattice, i.e., at J' = J. 


Inspired by the valence-bond-crystal states for the plateaus at m = ^ and |, it 
is natural to consider a corresponding valence-bond state at to = |: 

i4Bc>=ni^>,nit>. (25) 

j i 


with 


L,U) = [(^5 - 1)|| ttttw »0 - (V5 + 1)11 tmn»o + 2y/2\\ mm»o] -(26) 


Its energy per site at h = 0 is ey BC = — ~ —0.124 225 998 7J. Again, this is 

not the true ground state although it is not far from it. 

In Ref. 02] large-scale exact-diagonalization studies are compared with theoreti¬ 
cal predictions based on the variational wave functions for the plateau states. This 
comparison demonstrates that the introduced model states provide a good descrip¬ 
tion of the physics within the plateau regimes. In particular, the comparison of the 
(connected) spin and dimer correlation function, (sfs|) c = (sfs|) — (s?)(sj) and 
((s i ■ Sj)(s k ■ s i)) c = ((si • Sj)(s k ■ s;)) - ((s i ■ Sj))((s k ■ s;)), calculated using ex¬ 
act diagonalizations and the introduced valence-bond-crystal picture as well as the 
analysis of the low-lying excitations above the plateau ground states, yield clear ev¬ 
idence that the valence-bond-crystal picture holds also for the |- and the |-plateau 
states. To support this conclusion we present in Fig. [9] exact-diagonalization data 
for the overlaps ('I , ™ act |4/y BC ) of the valence-bond-crystal states with the exact 
ground states for m = | and m = § as function of ratio of bond strengths for 
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a finite lattice of N = 36 sites. Here J is the bond strength of the trapping cells 
corresponding to the valence-bond-crystal state and J' is the bond strength of all 
other bonds. Indeed the overlap is large in almost the whole region 0 < -y < 1. 
Only in the close vicinity of the uniform kagome limit (J' = J) it drops down to 

(^Jxactl^vBc) ~ °- 418 and (^fxactl^vBc) ~ °-357. Taking at J' = J instead of only 
one valence-bond-crystal state a symmetric linear combination of the three degener- 

5 5 

ate valence-bond-crystal states gives even a large overlap (H > exactl 4 ' s 9 ym vBc) ~ 0-724 

3 3 

and (^exactl^symVBc) ~ 3 .616. Note that for to = | the overlap is always unity, 

7 

since | v I'y BC ) is an exact eigenstate. We may compare the calculated overlaps with 
corresponding data for the Neel-long-range ordered s = | square-lattice Heisenberg 
antiferromagnet. For a finite lattice with N = 32 sites the overlap of a symmetric 
linear combination of the two Neel states is (^/gxactl'f'symNeei) ~ 0.203, i.e., it is even 
lower than those for the valence-bond-crystal model states. Thus the obtained large 
overlaps may serve as another indication of valence-bond crystalline order in the 
plateaus of the kagome antiferromagnet. 

Finally we note that, complementary to the approach of Ref. 0U S. Nishimoto 
et alPE2 used another route to the same problem. They developed a grand-canonical 
density-matrix-renormalization-group numerical technique to obtain the magneti¬ 
zation curve and to identify the nature of the ground states. The density-matrix- 
renormalization-group resultJSHI agree well with the valence-bond ansatzes and the 
exact-diagonalization data of Ref. 44j Moreover, Ref. [43] predicts an exotic (tiny) 
quantum plateau at to = y 

6.3. Localized-magnon ground states in zero magnetic field: The 

sawtooth chain with ferro- and antiferromagnetic interactions 

The sawtooth-chain (also dubbed as delta-chain or one-dimensional Tasaki-lattice) 
Heisenberg antiferromagnet is a well studied example exhibiting localized-magnon 
physics. If the exchange coupling along the zigzag path J 2 is two times larger than 
the exchange coupling along the basal straight line Ji, see Fig. [Q the lower one- 
magnon band is dispersionless and localized-magnon physics emerges at low tem¬ 
peratures around the saturation field h SSb t = 2 J 2 = 4Ji. Interestingly, very recently 
it has been found that the spin-^ Heisenberg model on this lattice with ferromag¬ 
netic zigzag J 2 < 0 and antiferromagnetic basal Ji = i| J 2 | > 0 interactions has a 
close relation to localized-magnon physics SSI (We notice here that in Ref. 147)1 the fer¬ 
romagnetic zigzag interaction is denoted by J\ < 0 and the antiferromagnetic basal 
interaction is denoted by J 2 = \\Ji\> 0.) Note that as = ( ferro-antiferromagnetic 
sawtooth-chain Heisenberg system may be use d for the description of the compound 
[Cu(bpy) (H 2 O)] [Cu(bpy) (mal)-(H 2 O)] (C10 4 )2 

Since the results of Ref. 06] may open a new window to find exact localized 
many-body ground states in quantum spin systems, we will illustrate here the main 
findings for the ferromagnetic-antiferromagnetic sawtooth chain in some detail. We 
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consider the Hamiltonian of the spin--| sawtooth-chain Heisenberg model 


H 


N_ 



(Following Ref. \4£\ we use here the one-index numbering of sites, see Fig. [I] i.e., 
m, 1 —> 2n — 1, m, 2 —> 2n and n = 1,..., ^.) For dominating ferromagnetic ex¬ 
change J 2 the ground state is the fully polarized ferromagnetic state. At a critical 
point a = -rjh = i a transition to a partly polarized ferrimagnetic ground state 
takes place. Just at the critical point the ground-state manifold of the model is given 
by a macroscopically degenerate set of localized-magnon-like eigenstates, where the 
ground-state degeneracy is even larger than that for the antiferromagnetic saw¬ 
tooth chain, where the latter one was calculated in Ref. [TO] A striking difference 
between the ferromagnetic-antiferromagnetic and the purely antiferromagnetic saw¬ 
tooth chain is the appearance of the localized ground states at zero field for the 
ferromagnetic-antiferromagnetic case. 

To construct the ground states at the critical point a = \ it is convenient to put 
J 2 = — 1 and to represent the Hamiltonian (l27l) as a sum of local Hamiltonians: 


N_ 

2 


3 


I 




(28) 


i=1 


Here i enumerates up-triangles which constitute the sawtooth-chain lattice whereas 
i\ and *3 are the (left and right) basal sites and *2 is the apical site in the *th 
triangle. The eigenvalues of Hi are: 0 (S = |), 0 (S = i), and | (S = |), and, 
therefore, the ground-state energy of H (l28l) is bounded from below by JA 0 = 0. 
Since the ferromagnetic state with maximal total spin S'max = y is an eigenstate 
of H [see Eq. ( 1271 ) ] with zero energy, the ground-state energy of H equals this lower 
bound, i.e., it is zero. 

Let us now illustrate how the complete (and linearly independent) ground-state 
manifold can be constructed.®^ It is convenient to consider the sawtooth chain with 
both open and periodic boundary conditions. First we consider the one-magnon 
subspace with S z = (7 — 1. There are y- linear independent localized states of the 
form ifii | f ■ ■ • t) w ith 



(29) 


V 7 * — s 2i A ^ S 2i+1 A S 2(i+1)’ * — 1, ■ ■ • , A/", Af — 2 (29, 

(periodic boundary conditions). Note here that similar states exist for the anti¬ 


ferromagnetic sawtooth chain, however, with the important difference, that the 
corresponding operator, ipf fm = — 2s ^ +1 + i contains a minus sign. For 

open boundary conditions (and odd N) the first and the last states in this se¬ 
quence are generated by the operators (pi = 2 s]“ + and = sj f _ 1 + 2 s]^, 
respectively, where N = 2AA. It is easy to prove that Hnpi\ f t) = 0 and 
Hi+\tpi\ f • ■ • t) = 0 hold. For all other i we can use the commutation to find 
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Hi(fi\ t ... t) = t • • ■ t) = 0. Hence <pi\ f ■ • ■ t) is an eigenstate of H with 

zero energy. The states created by operators pi are eigenstates of H and S z , but not 
of S 2 . Rather a state tpi | t ■ ■ • t) is a mixture of eigenstates of S 2 with S = ^ and 
S = -y — 1. On the other hand, states created by ipf im are eigenstates of S 2 . Note, 
that the symmetric linear combination JA ipi\ f ■ • ■ t) = 25“ | f ■ • ■ t) belongs to 
the subspace with S = There are ^ — 1 other linear combinations of ipi\ f ... t) 
(but not the states themselves) which belong to the subspace with S = ^ — 1. It 
is convenient to introduce another set of linearly independent operator functions 
(instead of pi defined above), namely 

m 

= to = 1,-N. (30) 

i=1 

The function <F(IM)| f ... t) belongs to the subspace with S = while the remain¬ 
ing -y — 1 functions yield contributions to the subspace with S = y- — 1. 

Next we consider the two-magnon subspace with S z = y — 2 and imply for 
simplicity open boundary conditions. Bearing in mind the two-magnon localized 
states of the antiferromagnetic sawtooth chain, one immediately finds that ipiipj | t 
... f) with j > i + 1 are exact ground states of the Hamiltonian H with zero 
energy, since they are exact eigenstates of each local Hamiltonian Hi with zero 
energy. The number of such hard-dimer states is given by the binomial coefficient 
= i(N — 1)(N — 2) (this is the canonical partition function of two dimers on 
the open chain of N sites). Another choice for the two-magnon ground states with 
zero energy is given by $(mi) ( h(m 2 )| t ■ ■ • t) with 1 < m\ < m 2 < N — 1 (there 
are also states of this kind). There are additional zero-energy ground states in 
the subspace with S z = ^ — 2 which belong to the subspace with S' = y — 1 and 
S = y : The operator function $(N) and mi < N — 1 yield the state $(mi)<l)(N)| t 
... t) = 2S _ 4 , (mi)| t • • ■ t) which belongs to the subspace with S = ^ — 1 and 
S = Y- The number of such states is One more state in the two-magnon 

subspace is given by $(N)<f)(N)| f • • ■ t) = (25' — ) 2 1 f ■ • ■ t) which belongs to the 
subspace with S = -j- In summary, the total number of the exact ground states in 
the two-magnon subspace is Cf^ + C ^ J „ 1 + -1 > c n-i- 

We pass to the fc-magnon subspace. We choose the fc-magnon states in the form 


<f>(mi)<f>(m 2 )... 4>(TOfe_i)$(TO fc )| t ■ • ■ t), 

1 < mi < m 2 < .. • < m/c-2 < mk -1 < mk < N — 1 (31) 

(this set belongs to the subspaces S > ^ — k, i.e., these states carry components of 
S 2 eigenstates with S = ^ — k and also components of S 2 eigenstates with higher 

5 ); 


$(mi)4>(m 2 )... 4>(m fc _i)$(N)| t ■ • ■ t), 
1 < mi < m 2 < ... < mk -2 < mk -1 < N — 1 


( 32 ) 
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(this set belongs to the subspaces S > y — k + 1); 

$(mi)$(m 2 ) ■ • ■ 3>(N)$(N)| f • ■ • t), 

1 < mi < m 2 < ... < mk-2 < N — 1 (33) 

(this set belongs to the subspaces S > A — k + 2); etc. and finally 

4>(N)$(N)...$(N)$(N)| f t) (34) 

(this set belongs to the subspace S = y). All the states are eigenstates of the 
Hamiltonian with the zero-energy eigenvalue. The total number of these states is 

ff° BC (k) = Cfi_ 1 + ... + C^_ 2 + Cft-i + dj-i> (35) 

see TableQ] It is important to note that the degeneracy 5 ° BC (fc) calculated for states 
of the form (1311) (l34)l is larger than that of the hard-dimer states constructed by 

the operators ipi defined above. The total number of degenerated ground states of 
the chain with open boundary conditions is 

N —1 N —1 

w OBC = 2 £ g OBC (fc) = 2 £ (N - k)C^ (36) 

k —0 k —0 

= 2[N2 n " 1 -(N-1)2 n - 2 1 = (N+1)2 n_1 —► e Nln2 —■» e^ ln2 . 

N—too N —^oo 

The resulting residual entropy (per site) at zero magnetic field for open boundary 
conditions is limjv-yoo Sq BC = limjv-voo lnW ^- — = \ In 2 . 

Table 1. Exact-diagonalization data for degeneracies: 


Open chain of N 


17 sites (i.e., 

N = 

9). 




S z — S 

17 

15 

13 

11 

9 

7 

5 

3 

1 


2 

2 

2 

2 

2 

2 

2 

2 

2 

k 

0 

i 

2 

3 

4 

5 

6 

7 

8 

a OBC 

9 17 

i 

8 

28 

56 

70 

56 

28 

8 

1 


Let us also briefly consider the case of periodic boundary conditions. Now instead 
of N = A±i we have W = y. Using operators E51) we can construct the ground 
states in the subspace S z = — k in the form <p ix ... ip ik | f ■ ■ • t) ■ The number of 

these hard-dimer localized-magnon states equals to (i- e -; ec l ua l s to the 

canonical partition function of k hard dimers on a simple periodic chain of N sites). 
As for open boundary conditions, again the hard-dimer states do not exhaust all 
possible ground-state eigenstates in the subspace with S' 2 = y — fc > 2. For 
example, for k = 2 we can check by inspection that ipi(ipi -1 + tpi + pi +i)l t - - - 1 > is 
a zero-energy eigenstate of the local Hamiltonians Hi, H i+ i, H L _\, and of all other 
ones. The authors of Ref. [46] conjecture (and check by exact diagonalizations for 
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finite chains up to N = 28) that the number of the ground states in the subspace 

^,..., -y is given by the formula 


with S z = — k for k = 0,1,2 , 


9N BC (k ) = 



o < k < 4 , 

4,y, 4r < k < AT 


(37) 


(the Kr oneck er delta corresponds to the resonating-valence-bond singlet 
eigenstatcA^I), see Table 0] From Eq. (1371) we get the residual entropy (per site) at 
zero magnetic field so = for periodic boundary conditions 

fE PBC = 2^ 1 5 P BC (fc)+^ BC (AT) -► ln2 , 

k—0 

Jim So BC = \ In 2 = lim s? BC . (38) 

N—>oo Z N—> oo 

Although, on first glance Eqs. (1551) and (1571) look quite different they yield finally 
for N —> oo the same result for the residual entropy. 


Table 2. Exact-diagonalization data for degeneracies: Periodic 
chain of N = 20 sites (i.e., A f = 10). 


S z = S 10 

9 

8 

7 

6 

5 

4 

3 

2 

1 

0 

k 0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

<?20 BC 1 

9 

35 

75 

90 

42 

0 

0 

0 

0 

1 


In Ref. 06] analytical arguments are supported by extensive exact-diagonalization 
studies. In particular, the gaps in fc-magnon subspaces were calculated. It was found 
that the gaps for the fc-magnon states with k > 2 decrease rapidly with increasing 
k. As a result (and by contrast to the antiferromagnetic sawtooth chain) the contri¬ 
bution of the excited states to the partition function cannot be neglected even for 
very low temperatures. 

The thermodynamics of the ferro-antiferromagnetic sawtooth Heisenberg chain 
at the critical point exhibits some interesting low-temperature features. There are 
additional low-temperature maxima in the specific heat, see Fig. 1 101 The susceptibil¬ 
ity diverges for T —> 0 according to % = c x T~ a with c x ~ 0.317 and a ~ 1.09. This 
divergence is naturally related to the ground states with finite magnetic moment. 
However, the exponent a is significantly lower than that for the one-dimensional 
Heisenberg ferromagnet «fm — 2 that is valid on the left side of the critical point, 
i.e., for Ji < J 2 | (see, e.g., Ref. Il2l ). The low-temperature behavior of the en¬ 
tropy per spin is given by \ ln2 + c s T x with c s ~ 0.245 and A ~ 0.12. The 

existence of a residual entropy leads also to an enhanced magnetocaloric effect. 

It is useful to summarize the d iffer ences in localized-magnon physics in the ferro- 
antiferromagnetic sawtooth chairPSI and the antiferromagnetic sawtooth chain. — 
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T 


Fig. 10. (Color online) Exact-diagonalization data for the specific heat per site of the ferro- 
antiferromagnetic sawtooth Heisenberg chain at the critical point (Ji = J 2 = —1)- The 
inset shows the low-temperature part in logarithmic scale. In the main panel the curves for 
N = 14,..., 22 practically coincide indicating small finite-size effects down to T « 0.001. 

Most importantly, the localized-magnon states are ground states in zero magnetic 
field, but not at the (typically large) saturation field. All the effects related to the 
localized-magnon states are therefore better accessible in experiment, in particular, 
the enhanced magnetocaloric effect is observed when switching off the applied mag¬ 
netic field. The excitation gaps above the degenerate ground-state manifold are very 
small. These low-lying excitations together with the ground-state manifold lead to 
characteristic power-law dependences for several thermodynamic quantities at low 
temperatures. 

7. Solid-State Localized-Magnon Systems 

The properties of localized-magnon systems are useful for explaining the low- 
temperature physics of frustrated quantum antiferromagnetic compounds at high 
magnetic fields. Most promising hitherto known candidates to show localized- 
magnon physics are the spin-i (distorted) diamond-chain compound azurite 
Cu3 (C0 3 ) 2 (OH) 2 33 and the spin-i W-like two-dimensional dimer magnet 
Ba2CoSi206Cl2. J ^ The recent experiments on these systems raise hopes for di¬ 
rect comparison with theoretical models elaborated within the localized-magnon 
context. 

7.1. Azurite 

The natural mineral azurite CuACOfi lAOH)? is a frustrated magnet displaying 
interesting magnetic behavior ! ^ l ^ | 49 | 50 | 51 | 52 j n particular, H. Kikuchi et alF-^ 
discovered a magnetization plateau at ^ of the saturation value followed by a 
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steep (jumpwise) increase to the saturation value in the low-temperature magne¬ 
tization curve. For the magnetic properties, the most important feature of azurite 
are diamond-likc chains which are formed by spin-i copper ions. There have been 
several attempts to derive a microscopic model for magnetic properties of azurite. 
In the Letter by H. Jeschke et alP^ the authors performed first-principle density- 
functional-theory calculations as well as model computations using the density- 
matrix-renormalization-group and exact-diagonalization methods and obtained a 
microscopic model for the frustrated quantum magnet azurite. As a result, they 
proposed an effective generalized spin-i diamond-chain model which provides a 
consistent description of a wide range of experiments. The Hamiltonian for the 
generalized diamond chain shown in Fig. [L] is as follows 
At 

H = 'y ' [J2S m ,l • S mj 2 + Jl (s m ,2 + S m +l,l) ■ S m> 3 + J3 (s m> l + S m +1,2) ' S mj 3 

m= 1 

3 

~\~Jm^m,3 ’ s m+1,3 99 ( 5 m,l S m,2 S m,s)] > ( s m,i) ■ (^9) 

Here the sum runs over all A f cells and the total number of spins is N = 3 M. H is 
the external magnetic field an d m b is the Bohr magneton. The exchange constants 
and the g -factor are as follows!^ 

Jr = 15.51 K, J 2 = 33 K, J 3 = 6.93 K, J m = 4.62 K, g = 2.06. (40) 

This model does not take into account a magnetic anisotropy which is clearly present 
in azurite®! anc j interchain couplings which are responsible for the low-temperature 
magnetic ordering of azurite. Nevertheless, the proposed model is a convenient start¬ 
ing point for state-of-the-art numerical many-body calculations and derivations of 
effective low-energy theories. 

The set of parameters for the spin model of azurite does not fit exactly to the 
localized-magnon conditions, since J± ^ J 3 and J m ^ 0. However, we are not far 
from the localized-magnon (i.e., flat-band) point in the parameter space and one may 
expect to see some remnant features of the localized magnons at low temperatures 
and high fields. Starting from Eq. (l39l) it is possible to derive a simpler spin model 
in the low-temperature high-field regime which mimics the properties of azurite!®^ 
By means of perturbation theory, see Sec.0 we arrive at the (pseudo)spin-i XXZ 
Heisenberg chain with the Hamiltonian 
At 

n cS = ^2 [C - ^ T m + J 0CT * +1 + « +1 ) + JT* +1 ] , (41) 

m=1 

cf. Eq. Using the strong-coupling approach)!^! supplemented by an analysis 

of the exact-diagonalization data of the initial full model (l39l) for excitations on 
the -|-plateau and the other experimentally accessible features of the magnetization 
curveone may find 

h = (1.384H — 44.860) K, J = 1.249 K, J Z = 0.565K. 


(42) 
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Here H is the experimentally applied magnetic field measured in teslas. On the other 
hand, within the localized-magnon perturbation theory (and neglecting J m = 0) one 
may finJED 

h = (1.384H — 44.778) K, J = 0.845 K, J z = 0.287 K, (43) 

cf. Eqs. (fl3l) and (fill) . The effective low-energy Hamiltonian (TUI) may be useful for 
making theoretical predictions for low-temperature high-field properties of azurite. 

Among the observable properties of azurite which are related to localized- 
magnon physics there are the plateau at of the saturation magnetization and 
the further steep increase to saturation reached at about 32.5 qt4Z|48| ag we n ag 
the inelastic neutron scattering spectra on the ^-plateau (H « 14 T), see Ref. [49l 
illustrating a band of dimer excitations whose dispersion is strongly suppressed by 
the competition of Ji and J 3 . As it was predicted in Ref. [5T] (but not measured yet) 
azurite should also exhibit a large magnetocaloric effect in the high-field region. 

The successful description of available experimental data for azurite requires 
further theoretical studies. In particular, an accurate treatment of a magnetic 
anisotropy (magnetization curve depends on the direction of the external magnetic 
field) and interchain couplings (azurite orders below about 2 K) is clearly necessary. 
These issues are not accounted within the discussed above models and remain to 
be explored. 

7.2. Spin dimer magnet BazCoSizOeC^ 

Very recently^ E^CoS^OgC^ single crystals have been synthesized, their crys¬ 
tal structure has been determined, and the magnetic properties have been investi¬ 
gated. The magnetism of this compound is determined by the Co 2+ ions. Structural 
studies show that Ba 2 CoSi 20 gCl 2 is a two-dimensional coupled spin dimer system 
with an exchange network corresponding to the frustrated bilayer shown in Fig. [2] 
Furthermore, the Co 2+ ion is in an octahedral environment, and for temperatures 
much lower than 250 K its effective magnetic moment is represented by fictitious 
spin-i operators. When the octahedral environment exhibits tetragonal symmetry 
the effective interaction between fictitious spins is described by the spin-^ XXZ 
Heisenberg model with a strong easy-plane anisotropy, i.e., the effective exchange 
interaction is W-like. In Ref.[55]it was argued that the interaction of z-components 
is only about of 15% of the isotropic XY interaction. Importantly, the interdimer 
exchange interactions are strongly frustrated and satisfy the relations Ju = J 22 , 
J 11 + J 22 = J 12 + J 21 , see Fig. [2] The authors of Ref. l53lestimate that the intradimer 
exchange interaction J 2 is about 4.6 times larger than the interdimer interactions, 
i.e., J 2 > 4 J which is required to realize the localized-magnon scenario. Certainly, 
one has to expect (small) deviations from the ideal frustrated-bilayer geometry, 
however, the main peculiarity is a strong W-like magnetic anisotropy. 

Let us briefly discuss some of the experimental results reported in Ref. [53] The 
performed electron-spin-resonance measurements on Ba 2 CoSi 20 eCl 2 at low tern- 
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peratures allow to obtain the g-factors. The magnetization process was measured 
at T = 1.3 K up to a magnetic field of H = 70 T. The saturation of the Co 2+ 
spin occurs at H sat ~ 57 and 41 T for H || c* (i.e., parallel to the c* axis) and 
H || ab (i.e., parallel to the ab plane), respectively; g c * ~ 2.0 and g a b — 3.86. For 
both field directions, there is a wide zero-magnetization plateau indicating a gapped 
non-magnetic ground state. The authors argue that two Co 2+ spins located on the 
bases of neighboring C0O4CI pyramids are coupled to form an antiferromagnetic 
dimer singlet. The plateau at zero magnetization is followed by a jumpwise change 
of the magnetization to a next plateau at half of the saturation magnetization. Then 
again a jumpwise change to the saturation magnetization is found. The transition 
fields where the jumps take place are He = 46.8 T and Hj] at = 56.7 T for H || c* 
and = 32.0 T and H^ t = 41.0 T for H || ab. The singlet ground state is stabi¬ 
lized in a wide field range below H c . Thus the measured magnetization curves for 
Ba 2 CoSi 2 C> 6 Cl 2 reported in Ref. f53l strongly resemble the theoretically predicted 
magnetization curve for the frustrated bilayer® 7 see the left panel in Fig. [4] 

Summarizing their experimental findings, the authors of Ref. [53] come to the 
conclusion that Ba 2 CoSi 20 gCl 2 is a two-dimensional coupled spin-dimer system 
with AW-likc exchange interactions. A stepwise magnetization process with plateaus 
at zero magnetization and at 1 of the saturation, irrespective of the magnetic field 
directions, together with the absence of a structural phase transition, show that the 
frustration for the interdimer exchange interactions is almost perfect. The ^-plateau 
state is almost exactly given by the alternate product of singlet and triplet dimers, 
which is called in Ref. [53] a “Wigner crystal of triplons’f] or a “Wigner crystal of 
magnons”. Obviously this Wigner crystal is just the magnon-crystal state, i.e., the 
state where the trapping cells are completely occupied by localized magnons in 
accordance with the hard-core rule, cf. Sec. [2TT1 and Sec. 16.21 

The above illustrated experimental results lead to the exciting suggestion that 
Ba2CoSi20gCl2 is a solid-state realization of the quantum Heisenberg antiferromag- 
net on a frustrated bilayer lattice discussed in Sec. [3] (see also Ref. [37). Certainly, 
some differences are also evident: The model discussed in Sec. [3] has isotropic ex¬ 
change interactions and ideal flat-band geometry. An extension of the localized- 
magnon picture for nonideal geometry has been discussed in Sec. [5] and apparently 
can be done for the frustrated bilayer too (at least within the strong-coupling ap¬ 
proach) . Within the framework of the strong-coupling approach it is s traig htforward 
to take into account also the anisotropy of the exchange interactions ! 1 ^ 

The most intriguing feature of the frustrated bilayer is a phase transition which 
occurs if the magnetic field is in the range hi < h < hi (see Sec. [31). In the ter¬ 
minology of the Ref. [531 that is the range H c < H < H sat . Therefore it would be 
very interesting to explore the magnetic field versus temperature phase diagram of 
Ba 2 CoSi 20 eCl 2 performing, for example, precise measurements of the specific heat 
and seeking for traces of the phase transition. However, it might be that the required 


The s z = 1 component of the spin triplet is called the triplon. 
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critical field H| =46.8 T is too high to perform the corresponding thermodynamic 
measurements, whereas for the other field direction, when H^- = 32.0 T, theoreti¬ 
cal predictions should be elaborated yet. [For the isotropic (i.e., XXX) Heisenberg 
interaction the held orientation is irrelevant.] One may also expect that the critical 
fields of Ba 2 CoSi 20 eCl 2 may vary under pressure, this way putting them into a 
better accessible region. 

8. Flat-Band Hubbard Model at Low Temperatures and Low 

Electron Densities 

We consider now electronic flat-band Hubbard models. Initially the Hubbard model 
was introduced to understand the origin of metallic ferromagnetism.^ Indeed, to 
avoid the repulsive on-site interaction U a parallel alignment of electron spins is 
favorable. On the other hand, typically the kinetic energy prevents the emergence 
of ferromagnetism. The suppression of the kinetic energy may open the route to 
ferromagnetism in the Hubbard model. Flat bands may lead to a fascinating example 
of ferromagnetism in electronic systems due to complete quenching of the kinetic 
energy. For the class of such models A. Mielkc^and H. Tasak^®' proved rigorously 
the existence of a fully polarized ferromagnetic ground state at certain electron 
numbers. 

To be specific, we consider the standard repulsive one-orbital Hubbard model 
with the Hamiltonian 

H = ^2 Ho t g. + U 
cr —t A * 

Ho,<7 — ^ ' tjj A A ^ ^ , 

to) * 

= A.cAter, ; Cj,crj } = ^<Tij } = 0 (44) 

with tij > 0. We have chosen in Eq. (14411 nonstandard sign conventions for the 
hopping integrals tij and the chemical potential /i to have closer analogy to the 
Heisenberg antiferromagnets in a magnetic field considered in the previous sections. 
Thus, hopping integrals Uj correspond to the antiferromagnetic exchange bonds Jij 
and chemical potential u corresponds to the applied magnetic field h. There are 
several reviews on Mielke-Tasaki Hubbard model s 0 | J |JJ| and, therefore, we touch 
below only some aspects which were less emphasized previously. Our basic strategy 
is to look for the correspondence between localized-magnon and localized-electron 
systems stressing the similarities and contrasting the differences. 

We begin with lattices, which, although having lowest-energy flat bands, do not 
show ferromagnetism in the ground state. Then we pass to some one-dimensional 
chains which exhibit Mielke-Tasaki ground-state ferromagnetism. In all cases we 
completely characterized the many-electron ground states up to a macroscopic num¬ 
ber of electrons n = n max oc N corresponding to half filling of the flat band. They 
constitute a set of independent ferromagnetic clusters. Furthermore, we calculate 
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their degeneracy. As a result, we obtain the thermodynamics of the standard Hub¬ 
bard model on these lattices at low temperatures and low electron densities. 

8.1. Frustrated diamond chain, frustrated ladder etc. 

In this section we summarize some results of Refs. ED msi mm There is a class of 
lattices which have a lowest-energy flat band, but do not belong to the class of 
Miclke’s or Tasaki’s flat-band ferromagnets. Examples are the frustrated diamond 
chain, the frustrated two-leg ladder, the double-tetrahedra chain, the frustrated 
three-leg ladder (see Fig. [[]), the dimer-plaquette chain etc. in one dimension or 
the square-kagome lattice (see Fig. 0 , the modified checkerboard lattice etc. in 
higher dimensions. The localized states which represent the flat-band states are 
located within traps, which, in contrast to Mielke’s or Tasaki’s lattices, do not have 
common sites, i.e., the traps are isolated from each other. 

We consider as an example the frustrated diamond Hubbard chain, see Fig. [T] 
The one-electron spectrum consists of three bands: 

t [t? 

£i = ~h + Ad £2,3 (k) = Y -j + 4tf (1 + cos k) + fi. (45) 

If 1 2 > 2t± the flat band with energy E\ is the lowest one. The corresponding 
eigenstates can be localized on the vertical bonds. Explicitly they can be written as 

4,Jvac) = -j= (cL,i, ct - <t Xa ) |vac) (46) 

with m running over all Af = y trapping cells, see Fig. [T] The application of n 
distinct operators l' m a to the vacuum state |vac) yields n-electron eigenstates with 
the energy nei which are the ground states in the n-electron subspace. For U > 0 
their degeneracy is easily calculated as gj^(n) = 2 where C%- = (^) = n! (n-tfy. 
is the binomial coefficient. Obviously, the traps (vertical bonds) are disconnected 
and therefore each traps can be occupied by an electron with arbitrary spin f or 
Surely, within the ground-state manifold there is the fully polarized ferromagnetic 
state, but non-magnetic states are predominant and averaging over all n-electron 
ground states yicldsPOl ( s \ n = —> 0. These arguments can be repeated for 

the other lattices mentioned above. We only have to change accordingly the flat- 
band energy and the critical value of i 2 above which the described localized-electron 
picture holds. In the case of the systems having triangular traps we have to take 
into account the additional chiral degrees of freedom to determine the degeneracy, 
see Sec. 0] As a result, each cell may be either empty, or occupied by an electron 
with spin a =t, 4 - and one of two values of the chirality \ = ±, or occupied by 
two electrons with different chiralities which belong to a spin-triplet state, that is, 
< 7 i (0) = 1, < 7 i(1) = 4, <?i(2) = 3. This circumstance will manifest itself in the explicit 
formula for the grand-canonical partition function, see Eq. (14811 below. 

The thermodynamic properties of all considered systems are similar. At low 
temperatures and for the chemical potential around rq = a* — £i they are dominated 
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by the huge number of localized-electron states. Their contribution to the grand- 
canonical partition function is given by 

N 

E(T,»,N) = 'f2g Ar (n)e^ =(l + 2zf, z = e~*. (47) 

71—0 

For the double-tetrahedra chain or the triangular tube due to the chiral degrees of 
freedom we get 

E(T,g,N) = ■■■ E 9i(ni)...giMz n ' + - +n " = (l + 4z + 3z 2 f m 

ni=0,l,2 7i _yv=0,1,2 

Knowing the grand-thermodynamical potential = — TlnH(T, /x, IV) one 

can easily obtain all thermodynamic quantities. For example, the entropy per 
cell is given by s(T, /i) = — , the specific heat per cell is given by 

c(T, g) = T , the average number of electrons is given by n(T, g) = dn ( 1 g'f N ' > ■ 

We can also obtain thermodynamic quantities in the canonical ensemble such as 
s(T, n ) or c(T, n). The thermodynamic quantities for the electron system show many 
similarities to the corresponding ones of the spin system. Thus, the entropy s(T, go) 
remains finite as T — > 0 (residual entropy), the grand-canonical specific heat c(T, g ) 
exhibits an extra low-temperature peak if g is below or above go, and the ground- 
state dependence of the average electron density versus the chemical potential ex¬ 
hibits a jump at g = go- More usual canonical quantities also provide fingerprints 
of localized-electron states. Thus, the entropy as a function of the electron density 
at low temperatures is especially interesting in the case of chiral localized-electron 
states, see Fig. HI] In the left panel of Fig. [FT] we show the residual ground-state 
entropy sl ~ T -^ n ^'> versus electron density for the Hubbard model on the double- 
tetrahedra chain as well as on the frustrated three-leg ladder with t 2 > 2<i (the data 
for both systems are identical). If the hopping integral between neighboring sites 
along the triangular traps acquires a small pure imaginary component ig (e.g., due 
to a magnetic field perpendicular to the triangular trap) the degeneracy due to the 
chirality is lifted and as a result the perturbed model exhibits a different dependence 
of the residual ground-state entropy on the electron density, see the right panel of 
Fig.CLl] 

Let us stress again that the considered lattices, although they have flat bands, do 
not show ferromagnetic ground-state order. However, ferromagnetic ground states 
may appear if the flat band acquires small dispersion provided the Hubbard inter¬ 
action is strong enoughfSS see Sec. [TU] below. 


8.2. Sawtooth and kagome chains 

In this section we summarize some results of Refs. mmm The ground-state man¬ 
ifold discussed in the previous subsection belongs to lattices with only disconnected 
trapping cells. As a result a trap may represent an isolated one-cell ferromagnetic 
cluster where the associated s = \ spin can be flipped irrespective of the other cells. 
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Fig. 11. (Color online) Entropy versus electron density for the Hubbard model on 

the double-tetrahedra chain as well as on the frustrated three-leg ladder with t 2 > 2fi (T = 0, left) 
and in the presence of small imaginary component of the hopping integral between neighboring 
sites of the triangular traps ig, g = 0.1 (slightly perturbed system) (T —> 0, right). Note that the 
curves for both systems are identical. 


Now we consider examples of Miclke’s line graphs and Tasaki’s decorated lattices in 
one dimension, where neighboring trapping cells share a common site, i.e., they are 
connected. We will consider here in some detail the sawtooth chain (one-dimensional 
Tasaki lattice) and a kagome chain (linegraph of an unfrustrated two-leg ladder), 
see Fig. |T] Again we will construct many-electron eigenstates and calculate their 
degeneracy. 

The sawtooth chain consists of J\f = -j cells (see Fig. [T]) and there are two 
branches of the one-electron energies: 

£ 1 , 2 (k) = cosk =F \Jt\ cos 2 k + 2 f|(l + cosk) + r. (49) 

For t 2 = \/2t\ > 0 the lowest one-electron band is flat, Ei(k) = Ei = —2t± + r = 
—2 1 + r. The flat-band eigenstates may be taken as the ones localized in V-shaped 
valleys, i.e., Z^-|vac) = (4j-i ~ ^ c lj + 4j+i)l vac )- 

The kagome chain consists of Af = y cells of diamond shape (see Fig. [T]) and 
there are three branches of the one-electron energies: 

£i(k) = £1 = — 2f + fJL, £ 2 (k) = 2tcOSN + H, £ 3 (k) = 2t(l + cos k) + r. (50) 

The localized states |vac) are created by the operators l] i:) = — c lj+i ~ c 3 j +2 + 

c 3 j+ 3 - Note, however, that for even A f, in contrast to the case of the sawtooth chain, 
a state with k = n from the dispersive band e 2 {k) has the flat-band energy. This 
eigenstate has the form a^ n |vac) = _ c 3 j+ 2 )l vac )- Being in 

this state the electron is trapped on the two legs (two-leg state). Thus we are faced 
with a set of M + 1 one-electron states having the energy —2 1 + rl. Note that a 
similar construction can be elaborated for another version of a kagome chain with 
■y cells of hexagon shape, see Ref. [21] 

We pass to the man y-electron states for Hubbard repulsion U > 0 and begin 

with the sawtooth chain P^EIES] 

In analogy to the last subsection, Sec. 18.11 we first 
consider those n-electron states where the electrons are located in non-contiguous 
(i.e., disconnected) traps. They are ground states in the n-electron subspace with 
energy nei. As long as we consider the occupation of disconnected traps only, the 
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AA 1 A 1 AA 

allowed allowed forbidden 

Fig. 12. (Color online) Illustration of the occupation rules of the trapping cells for the Hubbard 
model on the sawtooth chain. 


spin of the electrons is again irrelevant and we have a 2 n -fold degeneracy of such a 
state, if we fix the location of the trapped electrons (n independent ferromagnetic 
clusters, each consists of one cell). However, these states are not the only ground 
states in the n-electron subspace. An example of another type of states is given 
by 4ji t ■ • ■ 4(fi+n-i) tl vac )> this state corresponds to n electrons, all with spin 
sitting within n contiguous valleys, i.e., it represents a fully polarized ferromag¬ 
netic cluster of n connected cells, cf. Fig. [12] (left) for n = 2. If only ^-electrons 
are considered, the difference to the previously discussed localized multi-magnon 
states becomes particularly evident, since for the magnon states the occupation of 
neighboring trapping cells was not allowed. Although the traps have common sites, 
the Hubbard repulsion is inactive, since the electrons have the same spin, i.e., due 
to the Pauli principle the simultaneous occupation of a common site is ruled out. 
Clearly, because of the SU(2) invariance of the Hubbard model this state is (n-l-1)- 
fold degenerate, i.e., it is one component of a SU(2)-multiplet. Let us emphasize, 
that the occupation of neighboring cells with electrons of different spin is forbidden, 
i.e., such states are not the eigenstates, cf. Fig. [12] (right). To find all ground states 
in the n-electron subspace we must consider all possible splittings of n into a sum 
n = ni + n ,2 + ■ ■ ■, rii > 0 such that clusters of rii contiguous occupied valleys 
are separated. Furthermore, we have to include into the consideration the (n, + 1)- 
fold degeneracy of each connected cluster which gives an [(m + 1)(«2 + 1).. .]-fold 
degeneracy of a state with fixed position of the occupied traps. 

Next we consider the total degeneracy of the ground-state manifold in the n- 
electron subspace gj\f(n). There are two ways to calculate gj\[(n). The first approach 
was presented in Ref. [19} Interestingly, the localized-electron states can also be 
mapped onto a hard-dimer problem on an auxiliary simple chain. This mapping is 
similar (but not identical) and more sophisticated than the mapping for localized- 
magnon states mentioned in Secs. 12.ll andl3l Thus, in contrast to the mapping for 
Heisenberg spins each valley of the sawtooth chain is assigned to two sites of the 
simple chain, i.e., the auxiliary simple chain has 2J\f sites. Again, each state from the 
n-electron ground-state manifold corresponds to a spatial configuration of n hard- 
dimers on the simple chain. Only the easy case of maximum filling of n = n max = M 
electrons (half filling of the flat band) has to be treated separately. Finally we getP^ 
for the degeneracies 

g^{n) = Z hd (n, 2M) + (A f — l)5 n ,jv, (51) 

where Zh d (n, 27V") denotes the canonical partition function of n hard dimers on the 
2A/"-site simple chain. Using the classical transfer-matrix approach we calculate the 
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grand-canonical partition function 

S hd (*. 2AT) = Y * n 2hd(n, 2N) = Xf + X 2 ", 

n= 0 

/ 1 MO — M 

A ± = -±\j- + z, z = e~no = 21. (52) 

The canonical partition function can be determined through the relation 

Z hd (n,2A0= ^" 5h d d i: ,2A0 • 

2—0 

An alternative way to compute the contribution of the localized states to the 
partition function, Z{n,N) = Zh d (n, 2A/”), for the initial A/"-cell Hubbard chain 
was used in Refs. ED Em Let us illustrate this alternative approach in some detail. 
Within this approach we take into account not only the geometrical degeneracy 
related to various possibilities of occupying n cells of the A/”-cell chain but also the 
components of the SU(2)-multiplets. Instead of the true component of the multiplet, 
one may choose, for counting purposes, a representative state, e.g., the one where 
all spins f are in the left of each n contiguously occupied valleys and the spins | are 
on the right side. Then the counting problem becomes as follows: We have to count 
the number of configurations Z(n,Af) of three states in the trapping cells, namely, 
empty and occupied with a =t or a =■!, subject to the constraint that no states | 
is allowed to appear as the right neighbor of a j.. This problem can be solved with 
the help of a 3 x 3 transfer matrix in the grand-canonical-ensemble setup, namely, 

3(z, AT) = Y ~~ n Z{n,M) = TrT* = £+£ + 

n— 0 

= -+ z ± + z, £o = 0. (53) 

We have used a transfer matrix 

/T(0,0)T(0,t)T(0,|)\ /111\ 

T = T(t,0) T(t,t) T(t4) = [z zz , (54) 

\T(|,0) T(|,t) T(l,l)J \z0zj 

where z denotes the activity defined in Eq. (l52l) . Comparing Eqs. (l52ll and (l53ll we 
find that both expressions for E(z,Af) = S hd (z, 2A f) are equivalent, since = A^.. 

Now let us consider the kagome chain for electron numbers n > 1. There is a 
difference to the sawtooth chain, since in addition to the hard-dimer states (relevant 
for both systems) there are additional leg states. First we consider the case n = 2. 
In addition to the Z{2 1 M) hard-dimer states (electrons sit in the diamond-shaped 
trapping cells only) we can put one electron into one of J\f diamond cells and the 
other electron in a two-leg state. Taking into account the SU(2) symmetry we obtain 
3A f extra states in the subspace with n = 2 electrons compared to the sawtooth 
chain. Similarly, there are extra states for n > 2. In particular, localized states exist 
up to maximum filling n = n max = Af + 1, where all the A f diamonds and the 
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Fig. 13. (Color online) The Hubbard model on the sawtooth chain of IV = 12, 16, 20 sites: Specific 
heat per site for n = and U —>■ oo. The results for ideal flat-band geometry ( t 2 = \/2ti, dashed 
curves) are compared to the ones for slightly violated flat-band condition (£2 = 1.3£i, solid curves). 


legs are occupied. Note, however that these extra states become irrelevant for the 
thermodynamics in the limit N 00 , see below. The degeneracy of the ground 
states for the periodic kagome chain is given b}EI 

5 ka gom e (n ) = (1 _ ^ iA , +1 )^wtooth (n) + (1 _ Sn o) [(„ + 1)Q-1 + «J ni2 ] . (55) 


Knowing the ground-state energies and degeneracies for n = 0,..., n max , n max = A/" 
for the sawtooth chain or n max = M + 1 for the kagome chain, we can cal¬ 
culate the contribution of these states to the grand-canonical partition function 
E(T,fj,,N) = ]Cn=(T 9M{ n ) zn - This contribution dominates the thermodynamics 
at low temperatures for a chemical potential around the value /xo = 2 1. Although 
E(T, /x, N) is different for finite sawt ooth and kagome chains, it approaches the same 
function in the thermodynamic limifPH 


Z(T,v,N) = 




(56) 


Thus the sawtooth and the kagome chains exhibit identical behavior for N —> oo 
with a difference only in the relation between N and M. 

A detailed discussion of the thermodynamic properties and a comparison to 
exact-diagonalization data can be found in Refs. HTll and [2T1 H ere we focus only on 
the low-temperature dependence of the specific heat T A^ see Fig. [13] and on the 

ground-state magnetic moment, see Fig. 1141 The specific heat in the gran d-cano nical 
ensemble shows an extra low temperature maximum if /i is around /x 
specific heat of the localized states in the canonical ensemble vanishes, since the 
localized-electron energies depend only on the number of electrons [in fact, they 
constitute the (degenerate) ground-state manifold]. This is illustrated by exact- 
diagonalization results for the Hubbard sawtooth chain for n = ^ in Fig. [I3l 

where obviously C T (dashed curves) remains very small up to about y ~ 0.05. 
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However, a deviation form the flat-band condition, i.e., in case of £2 7 ^ v^ii, pro¬ 
duces a splitting of the ground-state manifold and an additional low-temperature 
maximum of the specific heat emerges as demonstrated in Fig. [13] for t 2 = l-3ii- 

The study of the integrated number of states above the degenerate ground state 
in the subspaces with different n < AT and U > 0 shows a noticeable density of 
states at rather low energies for small and even for rather large values of U. One 
may conclude that these low-lying excitations originate from states which belong to 
the ground-state manifold for U = O.^The question whether the density of low- 
energy excitations increases with increase of N leading to a vanishing excitation gap 
for N (X 0 remains open and requires further studies. For the sawtooth chain the 
separation of the ground-state manifold (for t 2 = V%ti) within the grand-canonical 
setup is controlled by the charge gap A/i = [E(A/"+l) — E(Af)} — [E(AT) — E(AT — 1)]. 
Exact-diagonalization calculations showSl that it opens linearly in U and reaches 
a saturation value as U —t 00 showing no visible finite-size effects. 

Next, we consider the ground-state magnetic properties of the sawtooth chain 
and the kagome chain. Since both chains belong to the class of Mielke’s or Tasaki’s 
flat-band ferromagnets, there is an exact proof of a fully polarized ferromagnetic 
ground state for n = n max , see Refs. [15, 16:, [17| Based on our discussion given 
above there is a simple geometrical explanation: For n = n max all cells are occu¬ 
pied and we are faced with a fully polarized ferromagnetic connected cluster that 
includes all cells. Hence no independent spin flips of disconnected clusters are possi¬ 
ble. Lowering n < n max just the possibility of independent spin flips of disconnected 
clusters leads to a decreasing averaged total magnetic moment. The mapping of the 
localized-electron states onto a classical transfer- matrix pr oblem is the crucial step 
to calculate the averaged total magnetic momentJH32I231 For that we consider the 
operator S 2 and perform the equal-weight average of S 2 over all degenerate ground 
states for the given number of electrons n. Then we consider the quantity 


0 < 


(S 2 ) n ^ S max (S max + 1) 


A f 2 


< 


Af 2 


5Vnax — 


(57) 


It achieves its maximal value, if the ground state in the subspace with n elec¬ 
trons is the saturated ferromagnetic state. If it has a nonzero value less than the 
maximal value (nonsaturated ferromagnetism) the ground-state manifold in the 
subspace with n electrons contains a substantial part of ferromagnetic ground 
states. Again both transfer-matrix techniques (see above) are suitable to cal¬ 
culate (S 2 ) n = 3 (S z2 ) n . We illustrate here briefly the approach based on the 
3x3 transfer-matrix technique, see Eqs. m and (l54l) . to calculate this quan¬ 
tity. The average over all degenerate ground states for a given number of electrons 
is ( S z2 ) n = A/’X]^Q 1 (S'oSj)n, where is the z-conrponent of the spin opera¬ 
tor of the trap j. Within the grand-canonical ensemble we may find = 

J 2 n= 2 ( n Z{n,Ar)(SQSj) n [note that here ( = is the activity which was de- 


d The absence of the gap would imply quantitative corrections to the hard-dimer description l56t - 
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Fig. 14. (Color online) The sawtooth chain (left) and the kagome chain (right): The averaged 
ground-state magnetic moment per cell ' as a function the electron density jjfor Af = 
16, 32, 64, 128, 256. 


noted by z in Eq. (1521) ] which yields the required canonical zz correlation function 


1 d"(SZS*) c 


{SoSj) n via the relation Z(n,Af)(SoSj)„ = m df „ 
3x3 transfer matrix technique is obtained by 


C=o 


Finally, within the 


/ 0 0 0 \ 

{S$S?) C = Tr(ST-'ST^), S= 0± 0 , (58) 

\ 00 -\) 

where T is the transfer matrix given in Eq. (1541) . 

Based on Eq. (1551) in Ref.[7T]the averaged magnetic moment ' was calculated 
up to M = 256, i.e., up to N = 512 (N = 768) for the sawtooth (kagome) chain (cf. 
Fig. 1171) . While for the sawtooth chain (left panel) only hard-dimer states exist, for 
the kagome chain (right panel) also the ferromagnetic states related to the existence 
of leg states have to be taken into account. However, their contribution vanishes 
in the thermodynamic limit, cf. also the discussion of thermodynamics. Obviously, 
the averaged magnetic moment becomes smaller with increasing of M for any fixed 
0 < jj < 1. In fact decreases linearly with jj, cf. Ref. [23]and also Eq. (f70l) 

in Sec. [9] where we discuss the flat-band ferromagnetism in one dimension using 
a percolation representation. Finally, the region of ground-state ferromagnetism 
shrinks to one particular electron density A- = 1 in the thermodynamic limit. In 
the next section we will illustrate that the existence of ground-state ferromagnetism 
in a finite region of electron density jj can be observed when going from dimension 
one to dimension two. 


9. Flat-Band Ferromagnetism as a Pauli-Correlated Percolation 
Problem 


A recently elaborated percolation representation of localized electrons in flat 
bands'^ provides an important step in understanding of flat-band ferromagnetism 
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and the corresponding para-ferromagnetic transition. Some ideas on the per¬ 
colation represent ation were introduced in the early papers by A. Mielke and 
H. TasakJiSll® — without presenting a detailed analysis. Here we report on the 
percolation analogy published in Refs. [22; 23. [24] [25] in some detail. As discussed 
in the previous section the emergence of ground-state ferromagnetism is related to 
ferromagnetic clusters of connected cells occupied by electrons. Due to the Pauli 
principle the occupation of neighboring cells (sharing common sites) by electrons, 
and therefore the occupation of a cluster built by connected cells, is possible only if 
such electrons form a symmetric (i.e., ferromagnetic) spin state. As the density of 
electrons increases, ferromagnetic clusters of increasing size appear. The ferromag¬ 
netic transition corresponds to the emergence of a ferromagnetic cluster of connected 
cells containing a nonzero fraction of the electrons, i.e., we face a percolation prob¬ 
lem. The special features of this percolation problem are related to the degeneracy, 
m + 1, of a ferromagnetic cluster C of |C| cells (|C| sites of the auxiliary lattice) 
containing m = |C| electrons. This degeneracy gives differing weights to different 
clustering of electrons. This unusualpercolation problem is called Pauli-correlated 
percolatiorPS or Pauli percolationPSl 

Within the canonical ensemble (i.e., for fixed electron density jj>) the percolation 
problem deals with all possible geometric configurations q of n particles distributed 
over J\f sites of the auxiliary lattice with a nontrivial weight of each configuration q 

M q 

W{q) = n (|C»I + x ) ■ (59) 

i—1 

Here the weight W ( q ) arises from the above mentioned spin degeneracy of a fer¬ 
romagnetic cluster of size \Ci\ in the configuration q. M q denotes the number of 
clusters in the configuration q, |Ci| = n. The expectation value of a quantity 

A is given by the usual expression (A) = A{q)W(q), where the canonical par¬ 
tition function is Z = W(q). Here the sum over q runs over all configurations for 

a given number of electrons n. As mentioned already above, the appearance of fer¬ 
romagnetism is related to the emergence of a cluster containing a nonzero fraction 
of elec trons. Thus it can be pinned down using standard tools of the percolation 
theoryP- ^ l 125 | 12 ]T] detect ferromagnetism the relevant observable is the square 
of the total spin S 2 , see also the last section. For n particles distributed over the 
lattice with a particular geometric configuration q it can be written as 



In the first form the contribution from each cluster Cj, i = 1,..., M qi of the con¬ 
figuration q is taken into account, while the second form contains the normalized 
number of cluster of size l in the configuration q denoted by n q (l). Such quantities 
[i.e., n q (l)\ are typical ones in percolation studies. The crucial differences is owing 
to the nontrivial weight factor W{q) (which reflects m + 1 possible orientations 
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of the total spin of a ferromagnetic cluster of m electrons) since for the standard 
percolation all geometric configurations have the same weight and W(q) = 1. 

Instead of the canonical ensemble with fixed number of occupied sites n (or 
electron density 4)-), one may consider the grand-canonical ensemble letting each 
site of the lattice to be occupied with an a priori probability p or left empty with 
an a priori probability 1 — p. Eq. (1591) now becomes 

M 1 M i , \ I Ci\ 

W(q) = l[e^(\C i \ + l) = l[ -T- (ICil + 1), (61) 

i— i i —i V P/ 

where is the fugacity which is used to tune the number of electrons. 

Eq. (RTT1) leads to the same formula for the expectation value, (A) = A{q)W(q), 

and the grand-canonical partition function, Z = W ( q ), where, however, now the 

sum over q runs over all configurations of n = 0,1,..., J\f particles (electrons). 

Pauli percolation exhibits a number of characteristic features. First, it corre¬ 
sponds to an effective (entropic) repulsion between particles which, although is 
long-range, saturates, is range free and does not depend on the shape of the cluster, 
and is many-body. Therefore it should not be astonishing that such quantities as the 
percolation threshold or the cluster-size distribution differ from the corresponding 
quantities for the standard percolation. 

Pauli percolation is interesting in its own right as a rare example of a discontin¬ 
uous (first-order) percolation transition with qualitatively more local interactions 
than in all currently considered examples, see Ref. [21] and references therein. It has 
a simple representation as a particular classical two-color, or contagion, percola¬ 
tion problem!^ Indeed Pauli percolation corresponds to the following variation of 
the standard one: (i) sites can come in two colors, say, green (uninfected) or red 
(infected), (ii) each site of a lattice is occupied and colored either green or red or 
left empty, (iii) only configurations where every cluster contains no more than one 
red site are taken into account. Since merging two clusters of size m and n reduces 
their overall weight from (m + l)(n + 1) to (to + n + 1), it closely resembles the 
“product rule” in the famous explosive percolation,but, unlike the latter, it does 
not require the elevated degree of non-locality in the configuration-update process 
as the new site in the system gets occupied, see Ref. EH 

Although the initial studies of Pauli percolation were performed on one and two- 
dimensional lattices, see below, in the recent study concerning infinite dimensions^ 
the analy tical and numerical calculations refer to a regular random graph of N 
sites.H2S1 Such graphs have a vanishing density of short cycles and mostly contain 
loops of size In N: hence they are locally tree-like. A rigorous analytical solution 
of Pauli percolation reported in Ref. (25] gives an explicit example of a rather rare 
discontinuous percolation transition in a random graph. A sketch of this solution is 
as follows. First the Pauli percolation is considered on a Cayley trecEl (setting for 


The distinction between a Bethe lattice and a Cayley tree is that the latter contains a root site 
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concreteness the coordination number z = 3) by solving the corresponding recursive 
relations in the two-color representation described above. Examining the probability 
Poo that the root site of an infinitely large tree is connected to its boundary one 
finds that it is zero (i.e., there is no percolation) until p < For p > | the nonzero 
(i.e., percolating) solutions emerge. Next, individual Cayley trees could be merged 
into the uniform Bethe lattice via the central site or bond addition. The value 
of p at which the percolation solution in the Bethe lattice Voo ^ 0 first emerges 
remains unchanged, p = although the new percolation probability itself changes, 
Voo 7 ^ Poo • The rigorous analysis in Ref. [25] illustrates the discontinuous nature of 
the percolation transition. The value of p at w hich it takes place is determined from 
the behavior of the central-site free energyP ^ | 131 | b een found that Voo ( P ) 

at p c = 0.823506 ... has a discontinuous jump between zero and nonzero values. 

The study of Pauli percolation on a 3-regular random graph leads to an in¬ 
teresting conclusion for understanding ferromagnetism. Namely, it implies that for 
the standard repulsive Hubbard model on the Tasaki-decorated 3-regular random 
graph the ferromagnetic ground state emerges as the electron density exceeds 
0.208176 .. . (see the red line in Fig. 3 of Ref. [25]) and its magnetization increases 
until it reaches a saturation value at ^ |. Interestingly, since the Bethe lattice is 

an exa mple of t he expander graph it intrinsically possesses a gap in the excitation 
spectraP^H34| w fo c h in turn can protect the ferromagnetism in such structures. 

The rest of this section on Pauli percolation contains a brief discussion of the ex¬ 
act solution on a simple chairP^and Monte Carlo simulations on a square lattice!^ 


9.1. Pauli percolation on a simple chain 

Alternatively to the approaches used in Sec. 18.21 for the investigation of the one- 
dimensional sawtooth-chain or kagome-chain Hubbard model we present here the 
percolation approach to consider these models, see Refs. Turn The treatment of 
the strongly correlated electron problem in one dimension as a Pauli percolation 
problem is useful (i) as a preparation for the treatment of the two-dimensional 
problem, see Sec. 19.21 (ii) to provide one-dimensional data to compare with two- 
dimensional ones, and (iii) as comparison with standard percolation. Moreover, it is 
valuable in its own right, since an exact solution of this unconventional percolation 
problem is possible. 

We analyze one-dimensional Pauli percolation using a transfer-matrix method.^ 
To illustrate the percolation approach and to provide results for comparison we first 
consider the standard percolation in one dimension which is exactly solvablel^SJ 


with a coordination number z — 1 while the bulk of a Bethe lattice is fully homogeneouJSD (see 
also Ref. tl29l . The standard percolation transition on the Bethe la ttice is the continuous one with 
the percolation threshold p c = —W, i. e ., p c = ^ for z = 3 [130|125 ] 

f ' In the Ref. l25l all the calculations are performed using p as an a priori probability of site to be 
occupied and colored green or red; it is connected to p in Eq. JED by the relations: p = and 


P 


p 

1 -p‘ 
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We work in the grand-canonical ensemble and seek for the grand-canonical partition 
function Z = Z(z,Af) of a percolating system as a sum of probabilities of all possible 
random realizations. Z(z,Af) can be written in terms of a 2 x 2 transfer matrix T 
as 


Z(z,Af) = TrT^, 



(62) 


The matrix element T(ni,rii+ 1 ) corresponds to the pair of neighboring sites i and 
i + 1 on a simple (periodic) chain and acquires the value 1 if the site i is empty or 
0 if the site i is occupied independently of the occupation number of site i + 1. The 
2 ^ terms of TrT' ,v ” correspond to all possible realizations of the percolating system. 
To determine the fugacity z let us calculate the average occupation number of the 
site i: 


( ru ) 


Tr(T-^N) /0 0\ 

z ’ ~ Vo i ) ' 


(63) 


The matrix elements of the matrix NT acquire the value 0 if the site i is empty or 
z if the site i is occupied this way counting only those realizations (among 2 ^ ones) 
of the percolating system that have an occupied site i. Diagonalizing the transfer 
matrix T one easily finds in the thermodynamic limit AT —> oo that (nf) = On 
the other hand, the average site-occupation number should be equal to p, ( m) = p. 
This yields the relation between z and p: = p or z = . 

We consider now the calculation of the average number of clusters of size l (nor- 
mali zed b y the lattice size Af) to be denoted by n(l). Simple arguments immediately 
yield:^^ n(Z) = (1 —p) 2 p l . To fix the cluster of length l, within the transfer-matrix 
method we start with an empty site, then we have a string (cluster) of l occupied 
sites, and the last site of the string is followed by an empty one. This has to be 
written in terms of matrices. Let us introduce the matrix S with only one nonzero 
element 1 if the site i is empty and the site * + 1 is occupied, the matrix C with 
only one nonzero element 2 if both sites i and i + 1 are occupied, and the matrix F 
with only one nonzero element z if the site i is occupied and the site i +1 is empty. 
To calculate n(l) we have to replace in Eq. (l62l) the product of a sequence of l + 1 
T-matrices by the product SC ,_1 F, i.e., 


n{l) 


Tr ( T .V_,_i sci _i F ) ^ _ ( 0 \\ 

Z ’ V 00 / 




Straightforward calculations in Eq. (1641) in the thermodynamic limit Af —> oo re¬ 
produce the cited above formul a n(l ) = (1 — p) 2 p l ■ Similarly we get the pair (site- 
occupation) correlation function^- 

TrfT-V~ i NT i N') 

9(1 ) = (n ini+l ) - ( ni )(n i+l ) = —^ -p 2 = P ( 1 - p)S t , o, (65) 


and the pair connectivity (the probability that two sites n and n + l are both 
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(66) 


occupied and belong to the same cluster) 

„ Tr^-'NC 1 ) , 

r(n, n + l) = -—- = p- p 


which were obtained by more elementary means in textbooks on percolationM 
Advantages of the transfer-matrix approach become undoubted while applying 
it to the Pauli percolation!^ We have simply to introduce the transfer-matrix for 
the Pauli percolation in one dimension 


/T(0,0)T(0,t)T(0,|)\ / 1 1 1\ 

T = T(t,0) T(t,t) T{ t4) \ = [zzz\, (67) 

\Ta,0)T4,t)T44)/ \z0zj 

as well as the other matrices emerging in the transfer-matrix calculations 

/000\ /011\ /000\ /000\ 

N = | 0 1 0 ) , S = | 0 0 0 ) , C= 02Z , F= zOO . (68) 

\0 0 1 / Vooo/ \00zJ [zoo/ 


These are now 3x3 matrices, cf. Sec. 18.21 because of the two spin values of an 
electron: The site can be empty, occupied by the electron with spin or occupied 
by the electron with spin 4-- The zero in the third row of the transfer matrix T in 
Eq. (1671) appears, since the configuration is forbidden. After straightforward 

although lengthy calculations one gets in the thermodynamic limit 


p(z) = 1 + 4z — 


1 + 4 z ’ 


zip) 


p(2 ~p) 
4(1 — p) 2 ’ 


n(l) 


4(1 ~pf 
(2 -P) 2 


{1 + l)a l , 


a = 


P 

2 -p 


g(l) = -(l-p) 2 a m , 


r(n, n + l) = p • 




a 


(69) 


The Pauli percolation shows a number of striking differences from the standard 
percolation. The most obvious differences can be seen in the average number of 
clusters of size l, n(l ), and in the pair correlation function g(l), see Fig. [15l Although 
the percolation threshold remains unchanged, p c = 1, the quantity n(l) for Pauli 
percolation has a maximum l* > 1 for p > 0.8 moving along l* « —(1 + ^2—) for 
p —> 1. For the standard percolation n{l ) decays monotonically with l. The pair 
correlation function g(l) decays with a correlation length £ = — 2I ^ a that diverges 
as ( p c — p) -1 when p —> p c — 1. Since the pair correlation function g(l) is negative, 
the interaction is repulsive, i.e., particle (electron) positions anticorrelate. (For the 
standard percolation there are no nontrivial pair correlations.) 

From the point of view of magnetism, Eqs. (l60l) and (l69l) allow to calculate the 
macroscopic magnetic moment (in the thermodynamic limit) 


<S 2 } 


3p(2 - p) 

8(1 -p) 


A f, p < 1. 


(70) 
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/ / 

Fig. 15. (Color online) The average number of clusters n(l) (left) and the pair correlation function 
g(l) (right) for p = 0.95 and p = 0.975 for standard (dotted red lines) and Pauli (solid blue lines) 
percolation in one dimension, see Eqs. ED , Ep and Eq. (1691) . 


This formula provides an extension to Af — > oo of the finite-size results reported in 
Fig. HU Clearly Eq. (I70l) implies that —>■ 0 if 0 < p < 1, as already discussed in 

Sec. 18.21 Only at p = 1 we have limj^oo > 0. 


9.2. Flat-band ferromagnetism, and Pauli percolation in two 
dimensions 

To study flat-band ferromagnetism of Hubbard electrons in two dimensions we con¬ 
sider the Tasaki lattice (Fig. [2J which is the two-dimensional counterpart of the 
sawtooth chain. The trapping cells, where the electrons can be localized, contain 
five sites, namely one site of the underlying square lattice [label (m, 1) in Fig. [2] 
and the four sites connected to site (m, 1 ) by the hopping integral t 2 (red zigzag 
paths in Fig. [2]). Since every cell of the initial Tasaki lattice of N sites contains 
precisely one site of the underlying square lattice, the corresponding lattice of the 
related Pauli percolation problem is just the underlying square lattice of M = ^ 
sites. 

The two-dimensional perc olation problem is more complicated and it can be 
studied numerically, onlyU^ However, having in mind that the initial model is a 
strongly correlated quantum many-body model the treatment of the relevant degrees 
of freedom by a classical percolation description provides a significant simplifica¬ 
tion of the initial problem. In contrast to the standard percolation, simple random 
sampling of geometrical configurations is not sufficient for averaging in the Pauli- 
percolation case, since various geometrical config urations q have different weights 
W(q). Going beyond standard numerical schemesJ1MI137] take ; n ^ 0 acc0 unt ef¬ 
ficiently various geometrical configurations q according to their weights W ( q ), one 
has to implement importance sampling choosing samples according to the distri¬ 
bution of W(q). In the canonical ensemble for a fixed number of electrons, a new 
configuration q 2 is generated from the given one qi by a random permutation of 
two sites in order to ensure the fixed number of electrons. In the grand-canonical 
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ensemble, a site is simply chosen and if it is empty (occupied), a particle is proposed 
to be inserted (removed). The new configuration is accepted with the Metropolis 
probability min[l, ]. In addition, for the grand-canonical simulations exchange 

Monte Carlo steps have been employedLl-SSi As a result, the numerical investiga¬ 
tion of the Pauli percolation becomes more challenging, since it suffers from critical 
slowing down (that does not occur in standard-percolation simulations). Labeli ng of 
clusters is done in two different ways: (i) using a modified Newman-Ziff algorith m^' ! 
which locally updates cluster labeling for f ixed number of occupied sites and (ii) 
using the Hoshen-Kopelman algorithnfL® which makes a global update. Central 
results of numerical simulations concern (S 2 ) and the average number n(l) of clus¬ 
ters of size l normalized by the size of the lattice A f. For the numerical simulations 
lattice sizes M = C x C up to C = 270 were used, where periodic or open boundary 
conditions were imposed. 

The average square of the magnetic moment (S 2 ) is the quantity which di¬ 
rectly signalizes the appearance of the ground-state ferromagnetism. In canonical 
simulations, one can obtain the square of the magnetic moment of the largest clus¬ 
ter M 2 = ^ Sm gy luster ^ , where S^ ax = S(§ + 1), or the square of the magnetic 

S max ^ ^ 

/g2\ 

moment A 7 as a function of electron density p = -S- = 3-S-. Furthermore, one 
can perform a finite-size scaling of the square of the magnetic moment to get 
liin/^oo M 2 , see Fig. [16] (left panel). We find that the scaling behavior changes 

from M 2 = a + bC~ 2 H-for low electron densities p to M 2 = a + bC~ x + • • • at 

higher p. While for p = 0.62 the macroscopic moment scales to zero with system 
size, for p = 0.70 it approaches a nonzero value of about 0.40. Further increasing of 
the electron density, i.e., to p = 0.78, yields an almost size-independent magnetic 
moment of about 0.98. Additionally, the cluster-size distribution n{l) indicates the 
emergence of a large component for p > 0.62, see right panel of Fig. [16] The average 
square of the magnetic moment M 2 is presented in Fig. [IT] for finite systems along 
with the extrapolated values. 

Altogether, from these findings one may conclude, that the critical density for 
the Pauli percolation exceeds that of the standard counterpart p c = 0.592746 ..., 
see, e.g., Ref. (1241 For higher densities, there appears a regime of unsaturated fer- 
romagnetism, which is characterized by 0 < ^ < 1. Canonical simulations of the 

'max 

Pauli percolation exhibit some features arising in a phase-coexistence regime. The 
high-density phase appears to form first as a compact nonpercolating object with 
a macroscopic magnetic moment. For higher densities, including p = 0.7, the fer¬ 
romagnetic phase spans across the system. Therefore, the phase-separated regime 
is drastically different from standard percolation. Grand-canonical simulations ex¬ 
hibit a jump at a certain chemical potential between densities (jj)~ around 0.63(1) 
and (tj-)+ ~ 0.75(2). All findings are consistent with the visual analysis of typical 
snapshots of configurations for the standard and Pauli percolation, see Fig. [18] For 
electron densities below the threshold the system tends to form many small clus¬ 
ters which increase the weight. Slightly above the threshold the larger cluster is 
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Fig. 16. (Color online) The Hubbard model on the two-dimensional Tasaki lattice. Left: Finite- 

/ S 2 \ 

size scaling of the square of the magnetic moment ^—-. Right: Normalized number of clusters 

'max 

n(l ) for C = 200. 



Fig. 17. (Color online) The Hubbard model on the two-dimensional Tasaki lattice. Square of the 

/ S 2 v 

magnetic moment of the largest cluster, M 2 = — maxciuster _ versus electron density P = ■ 


still suppressed by the existence of many small clusters which help to increase the 
weight of a configuration. 

In summary, the simulations of the Pauli percolation in two dimensions show 
the following. Effective repulsive interaction leads to a breaking up of the clusters, 
and thus to a first-order grand-canonical transition in two dimensions. The criti¬ 
cal density is higher than that of standard (site) percolation. For the underlying 
two-dimensional Tasaki-Hubbard model these findings imply ground-state ferro¬ 
magnetism in a range of electron density from 0.21(1) to 
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Fig. 18. (Color online) Snapshots of configurations for standard (upper row) and Pauli-correlated 
percolation (middle and bottom rows) on a square lattice of size jCxjC, £ = 200, for small deviations 
from the critical concentration p c . The pink color corresponds to the largest cluster for the given 
concentration. The concentrations (from left to right) are pi = 0.574 and p 2 = 0.6 for standard 
percolation, where p c = 0.592746..., and p 3 = 0.62 (paramagnetic), p 4 = 0.65 (phase-separated), 
= 0.7 (phase-separated), and pe = 0.78 (ferromagnetic) for Pauli percolation. 
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10. Dispersion-Driven Ferromagnetism in a Flat-Band Flubbard 
System 

Introducing dispersion (violation of the flat-band geometry) typically modifies the 
balance of interaction and kinetic energy and tends to destabilize ferromagnetism 
in the Hubbard model. However, it was demonstrated by several studies that the 
flat-band ferromagnetic ground states are robust if the flat band becomes (slightly) 
dispersive and the Hubbard repulsion U is larger than a threshold U r > 0, where U c 
depends on the degree of the deviation from the flat-band geometry!^ By contrast, 
recently it has been found in Ref. [54] that unexpectedly for the class of flat-band 
Hubbard systems without ground-state ferromagnetism, see Sec. 18.11 just the dis¬ 
persion, i.e., the kinetic energy, will open the route to ferromagnetic ground states. 
This dispe rsion-driv en ferromagnetism resembles the famous order-from-disorder 
mechanism} - - —^4 51 i.e., due to distortions an ordered ground state is selected from 
the degenerate flat-band ground-state manifold, where paramagnetic states domi¬ 
nate. 

To this end we investigate the frustrated diamond chain, where for ideal geom¬ 
etry the lowest one-electron band is flat. However, the ground state even for full 
occupation of all trapping cells is nonmagnetic for arbitrary on-site repulsion U, 
since the trapping cells for the frustrated diamond chain are isolated from each 
other and the resulting huge set of paramagnetic states prevails against the ferro¬ 
magnetic eigenstate, cf. Sec. 18.11 Below we will illustrate how a dispersion of the 
former flat band leads to ground-state ferromagnetism. 

We consider the Hubbard model (l44l) on a distorted frustrated diamond-chain 
lattice, see Fig. [l] i.e., we assume ti ^ 1 3 but fix their average ti + £3 = 2 1. The 
electron density is fixed to n = J\f , i.e., for ideal geometry all trapping cells are 
occupied and the flat band is half filled. For the sake of comparison we first give 
the bandwidth Wd of the dispersive band above the flat band for t\ = £ 3 : Wd = 

(*3+*l) 2 


| +8i 2 -f 


(i 


t v t -C 1). The bandwidth Wf of the former flat band due 

to distortion t\ ^ £3 is Wf ns . We introduce the dimensionless parameter 

. Then Q 2 ns is a appropriate measure of the strength of deviation 
from the ideal flat-band geometry. 

To determine the ground-state properties, in Ref. 1541 exact diagonalization and 
fourth-order perturbation theory are used. The main result, namely the ground- 
state phase diagram, is illustrated in Fig. [T9] For not too large deviations from 
ideal geometry controlled by H 2 < 1 and for sufficiently large on-site repulsion 
U > U c , the ground state of the Hubbard diamond chain with half-filled lowest 
band is ferromagnetic (the region denoted as “FM”). 

The Lanczos exact diagonalization for distorted chains were performed up to 
IV = 30 to determine the critical value U c above which the ground state is the 
ferromagnetic SU(2)-multiplet with (S 2 ) = + 1). Note that for the largest 

system N = 30, U = 00 was assumed to find the maximal dispersion fl 2 ns 0.884 
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Fig. 19. (Color online) Phase diagram of the distorted Hubbard diamond chain with electron 
density -ji = (half-filled lowest band). Ferromagnetism (denoted by “FM”) appears for on-site 
repulsions U above a critical value U c . U c is shown as a function of the dimensionless bandwidth 
parameter Q?. The various critical lines U c (£l 2 ) are obtained by fourth-order perturbation theory 
(symbol “PT”), see Eq. cm. and by exact diagonalization for finite systems of N = 5 (open 
boundary conditions) and N = 12,18, 24 (periodic boundary conditions) for hopping parameter 
sets with £3 +£1 =2 and £2 = 3 and £2 = 6 . 


above which ground-state ferromagnetism does not exist. There is an obvious finite- 
size dependence of t/ c (fi 2 ), but the general shape of the curve is the same for all 
TV , cf. Fig. 1191 Interestingly the region of ferromagnetism increases with growing 
system size TV. That provides evidence that the dispersion-driven ferromagnetism 
does not disappear for TV — > oo. The influence of the vertical hopping integral t 2 on 
the phase boundary U C (Q 2 ) is visible from the corresponding curves for TV = 18. For 
small bandwidth f l 2 the region of ferromagnetism slightly shrinks with increasing 
of < 2 , whereas U c decreases with increasing t 2 for H 2 > 0.25. 

Complementary to the exact-diagonalization study the fourth-order perturba¬ 
tion theory was used. All essential features of the perturbation approach are already 
provided by a 5-site cluster that contains two trapping cells (the sites m, 1, m, 2, 
7 Tt,3, m + 1,1, and m + 1,2 connected by six bonds, see Fig. [T]) in the subspace 
of n = 2 electrons. The unperturbed Hamiltonian contains the t 2 - and t/-terms, 
the ti- and the f 3 -terms constitute the perturbation. The ground-state energy of 
the unperturbed Hamiltonian is —2^2, it is fourfold degenerate (one singlet, and 
threefold degenerate triplet). The perturbation theory up to fourth order gives for 
the triplet state 


E± — —2 tn — 


(t3 — ti) 2 (ts + fi) 2 (f3 — fi) 2 (t3 — ti) 4 


h 


2 


t 6 


(71) 


i.e., E t is independent of U. In contrast, the energy of the singlet state, E Sl depends 
on U. The fourth-order result can be found in Ref. [5U In the small-!/ limit the 
dominating term in E S {U ) is 

2(i 3 — H) 4 


E b {U->0) = — 


ut 2 


(72) 
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whereas in the limit U 


E S (U —> oo) = -2t 2 ^ 


oo it is given by 

(^3 — ^l) 2 (*3 + *l) 2 (*3 — *l) 2 


*2 


4*1 


3(*3 — *l) 4 

4*1 


+ ....(73) 


Obviously, we have E s < E t for small U, but E s > E t for large U, i.e., there is a 
critical value U c above which ferromagnetism emerges. Moreover, Eqs. ED E3 
show that the appearance of ground-state ferromagnetism for deviations from ideal 
geometry is a result of fourth-order processes. From E t = E S (U C ) one gets a compact 
analytical formula 


U c V16 + 65fl 2 + ^ 

* 2 1 - n 2 [ ’ 

From Fig. [T9] one can conclude, the perturbation theory confirms the exact- 
diagonalization results for the existence of ferromagnetism. The quantitative agree¬ 
ment between perturbation theory and exact diagonalization is the better the 
smaller 9 2 and the larger * 2 . Higher-order perturbation theory naturally would 
enlarge the region of quantitative agreement with exact-diagonalization data. 

In summary, the exact-diagonalization data and the perturbation theory show a 
dispersion-driven appearance of ground-state ferromagnetism of the Hubbard model 
on a frustrated lattice. For a half-filled lowest one-electron band, the ground state 
for the ideal flat-band case has perfect charge order [n = AT electrons are localized 
in A f trapping cells (vertical dimers) thus avoiding the on-site Hubbard repulsion] 
but no spin oder (spins of electrons in different cells are uncorrelated leading to 
2^-fold degeneracy of the ground state). A small deviation from the ideal flat- 
band geometry removes the localization of the electrons and they (i) either form a 
complicate many-electron ground state with a U -dependent energy which increases 
as U increases (when U is small) (ii) or form a symmetric (i.e., ferromagnetic) 
spin multiplet ground state with a ^-independent energy thus avoiding U (when 
U is large enough). We mention finally that such a scenario is not restricted to the 
specific geometry of the frustrated diamond chain and may be expected for other 
lattices with isolated cells discussed in Sec. l8.1l such as the frustrated ladder or some 
two-dimensional latticesP^U 


11. Closely Related Studies on Electron Systems 

11.1. Construction of exact ground states of the Hubbard model 
using positive semidefinite operator properties 

Alternatively to Mielke’s or Tasaki’s description of the flat-band scenario and to the 
route to ferromagnetism described in the previous section, Z. Gulacsi, A. Kampf, 
D. Yollhardt and coworkers use a technique to construct exact ground states (includ¬ 
ing ferromagnetic o nes) of the Hubbard model that is based on positive semidefinite 
operator propertiesPSESl j n w hat follows we will illustrate some basic features of 
this method. For more details on the method and more applications to strongly 
correlated electron systems we refer the interested reader to Refs. (55j [561 Instead 
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of dealing with an (arbitrary) Hamiltonian H we may always consider the operator 
O p = H — Eq, where Eq is the ground-state energy, which is a positive semidefinite 
operator with only non-negative eigenvalues. Therefore we can start by an exact 
rewriting of the initial Hamiltonian in the form 

H = Y J Op,i + C g , H , (75) 

i 

where the c-number C g ,h is a function of the Hamiltonian parameters and O v y are 
positive semidefinite operators. It is reasonable to introduce block operators A na 
which are composed of operators acting on the sites of a finite block n. Then the 
operators O p ^ contain types of expressions of the form A\ a A n ^ a which preserve 
the positive semidefinite form of O p ,j. To get a valid transformation to the form 
H = )>A 0 Pt i + C g ^H a specific relationship between the parameters of the initial 
Hamiltonian and those contained in the expressions of the O p j in terms of the 
block operators A n cr (the so-called matching equations) must be fulfilled. The rea¬ 
son why Eq. m merits attention is that the deduction of the ground state |GS) can 
be performed by constructing the most general state which satisfies the equation 
O p |GS) = 0. Solving this equation opens a route for deducing exact results for non- 
integrable systems independent of dimensionality. The construction of exact ground 
states is then performed in four steps: (i) the transformation of the Hamiltonian in 
a positive semidefinite form, (ii) the solution of the matching equations, (iii) the 
construction of the ground states, (iv) the proof of the uniqueness. In a final step 
then the physical properties are deduced from the obtained ground state. 

As an illustrative example we consider a diamond Hubbard chain with the Hamil¬ 
tonian 


H = H 0 + Hu, 

Hq = 'y ' y ' | te 3 (^ c ln,l,a c m—l,3,a + c ln,3,a c m,l,a 


rr—t 4 m 


A C m,2,cr C m— 1,3,0- + C m,3,tr C "i,2,er^ 
+t±cln i a c m,2,a + t\\ c m - l,3,rT C m,3,a + H.C. 

( C rn,l,<T C ™,,l,cr + C m,2,tr C m,2,o-^ | ; 

Hu = U y ( {jlm,A T * 


(76) 


For the case t± = ty = 0 and 5 = the one-electron spectrum consists of three flat 
bands (“Aharonov-Bohm cage” limit). The exact ground states for electron densities 
ji = -■ and i represent then simply an example of Mielke’s and Tasaki’s flat- 
band ferromagnetism, see Ref. 1551 This example, to some extent, can be understood 
without introducing the positive-semidefinite-operators technique: The O p operator 
is simply a diagonalized Hq and the ground states |GS) are trivially constructed 
from the localized-Wannier-eigenstate operators. 
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One can obtain exact many-electron ground states also for other special cases 
which do not belong to the scenario of Mielke’s and Tasaki’s flat-band ferromag¬ 
netism. We illustrate as an example the case fy > 0, e = — 1 ± + —, and S = n; 
then the one-electron spectrum consists of a dispersive lowest band and two upper 
flat bands. In order to rewrite Hq as a positive semidefinite form one defines the 
noncanonical fermionic block operators 


Am,(7 — 1,3,(T “b ^2Cm,l,cr “b U3Cm,3,cr “b U4Cm,2,(7: (77) 

which fulfill (A m , CT ) 2 = 0 and ^ 6 mi ,m 2 - Using Eq. ((771) one gets 

the expression for A' m a A m ^ and with the requirement — Tim a ^tn o-Am.o- = Hq 
one arrives at the matching conditions a|ai = <2302 = a|a 3 = a\a± = —te 1 ?, 
0 , 2(14 = —t_ l, ajar = — 1 ||, |ai | 2 + la. 3 1 2 = e + |a, 2 1 2 = e + |a. 4 1 2 . As a result, Eq. (1771) 
becomes 


Am.<7 — VH ^m— l,3,er Cm,3,er 2tj_C 2 (Cm,2,er Cm,l,er) 


(78) 


Consequently, the Hamiltonian has been transformed into the positive semidefinite 
form 


H — Am^A^ a + UP + Eqs, 

m,a 

P = ^ ^ (jlm,f 1 ) (j^m, 1 ) 7 
m 

Eg s = (e + U +1 ±) n — ^3 U + 4tj_ + -—^ A/"; M = —. (79) 

For electron number n = 4A f the ground state is 

igs)= c ^i vac ) - F i =n ( 8 °) 

where C is a normalization constant, m runs over all M cells, and i\ ^ ii are 
two arbitrary sites in each cell. The ground state can be also found for higher 
electron densities ^ > | The physical properties of these ground-state solutions 
correspond to a correlated half-metal. 

Let us mention that the procedure described above is not restricted to the case 
when only flat one-electron bands are present. For instance in Ref. [5B]the authors 
provide a solution for a model, where all one-electron bands are dispersive. 


11.2. Flat bands and randomness 

In flat-band systems with perfect geometry the electrons are localized due to de¬ 
structive quantum interference. On the other hand, the standard Anderson local¬ 
ization of electrons in a dispersive band is caused by disorder interp lay 

between flat-band localization and disorder may lead to unusual features E35BI59] 
A spectacular effect is a “inverse Anderson transition” found in three-dimensional 
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flat-band systems, i.e., a disorder-driven delocalization of electrons!^ To demon¬ 
strate this effect, the authors of Ref. consider a system with only flat bands 
- a diamond lattice with fourfold-degenerated orbitals on each site. The system 
is described by the standard tight-binding Hamiltonian and the disorder is intro¬ 
duced via mutually independent random on-site energies. The random variables 
are uniformly distributed over the range from — -y- to To study the transi¬ 
tion of the one-electron states the authors examine the level statistics. They found 
a disorder-induced localization-delocalization-localization transition which in this 
specific model is characterized by the two critical values of the disorder strength, 
W c i « 13.5 and W c 2 « 37.0. While at W C 2 the well-known Anderson transition 
takes place, the first transition at W c \ < W c 2 is an unknown new transition from 
localized states to extended states driven by disorder. The localization for W < W c 1 
is not due to the strength of disorder but to the flat-band localization which survives 
under a certain strength of disorder. 

A further analysis of tight-binding flat-band systems in the presence of disorder 
was reported in Ref. [ 55 ] The authors of Ref. 1551 focus on the weak-disorder limit and 
consider the two-dimensional checkerboard lattice (also called the planar pyrochlore 
lattice or the square lattice with crossings). For the random on-site potential a 
Gaussian distribution is assumed, periodic systems of size L x L, 34 < L < 80, are 
examined numerically. To study localization they use the probability distribution 
of the spacing between adjacent eigenvalues: Extended and localized phases are 
identified by Wigner-Dyson and Poisson distributions, respectively, while the critical 
point corresponds to a distinct, universal distribution. The numerical results for 
the probability distribution of the spacing between adjacent levels is intermediate 
between those for Wigner-Dyson and Poisson distributions, which indicates that the 
system is at a critical point. Thus, weak disorder in a two-dimensional tight-binding 
model with a flat band gives rise to eigenstates that are critical, i.e., they are neither 
Anderson localized nor spatially extended. 

In a series of recent papers 0 ™ the interplay of flat-band geometry and randomness 
has been studied for several lattices in detail. As already mentioned in Sec. 12. li the 
authors also give a scheme for generating flat-band lattices by adding Fano defects to 
an initial lattice with dispersive bands and performing appropriate transformations 
of the Fano lattice to a corresponding flat-band lattice, see Ref. [59: 

In what follows we will illustrate some results of Ref. ES Let us consider as an 
example a tight-binding diamond chain (with zero “frustrating” vertical bond £ 2 , 
see Fig. |TJ). The corresponding one-electron band energies are {0, ±2\/l + cosk}. 
Note that t\ = fa = 1 and £2 = 0; then the flat band has zero energy. We con¬ 
sider the eigenstate of the Hamiltonian of the form IT) = X^m( a ™ c m,a + b + 
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Cmcln C )|vac)@ The basic equations to be discussed further are as follows: 

iu n €n,aQ j n = V ^n+1? 
itn + €n,bbn ~ V ( a n 4” ^n)i 

iCn T Cn,c^n — V ^n-t-1- ($1) 

Here V 2 / n = /„ + f n -i and e n ,a> ct = a,b,c , is the random potential. The tight- 
binding model described by Eq. (1811) contains intersecting dispersive and flat bands 
at zero energy. The disorder is introduced as usual by uniformly distributed random 
on-site energies e njQ , £ [—^-] • 

Various standard quantities such as the participation ratio, the compactness 
index, the localization length £ (asymptotic decay rate of the eigenmode tails) etc. 
were calculated numerically by diagonalization of Eq. ED- The authors of Ref. (59l 
found different behavior for the low-energy (|E| < and high-energy (|E| > 4T) 
regimes. In particular, for the correlation length at the flat-band energy E = 0 
£(IV) ~ W _1 with the unusual exponent 7 = 1.30 ± 0.01 was found in contrast 
to the exponent 7 = 2 for a dispersive band, i.e., the flat-band states are more 
delocalized. Moving slightly away from the flat-band energy the anomalous exponent 
persists as long as \E\ < W. As E increases further, there is a rapid crossover to 
the conventional 7 = 2 exponent. Comparing the behavior of other measures of 
localization, the authors conclude that the flat-band states display criticality in the 
weak-disorder regime. 

As mentioned above, the authors of Ref. [59] found a relation between flat-band 
lattices and lattices with Fano defects. Using the example of the frustrated two- 
leg ladder (see Fig. [D note that in Ref. [59] the notation “cross-stitch lattice” is 
used instead of “frustrated two-leg ladder”) they discuss the connection between 
the existence of flat-band states and the appearance of Fano resonancei^l f or wave 
propagation. This relationship between the representations allows to relate flat-band 
physics to features of Fano resonances.^! For the frustrated two-leg ladder with ti = 
1 the one-electron eigenstates can be written as IT) = X^m( a ”i c m a + ^m c L t,)l vac )0 
To find the coefficients a m and b m one has to solve the equations 

Ea n — Cn,a^n ^n+l a n— 1 b n — \ h n +\ t 2 b n: 

Eb n — € n ,bbn Un+1 On-1 b n —\ ^n+1 ^2 (82) 

In the absence of a potential e ra> b = e ra> b = 0 the resulting energy bands read 

-Ufb = t 2l E(k ) = —4cosk — t 2 . 


E We use here the notations (indices) of Ref. [591 To make the relationship to the notations (indices) 
used throughout our paper (see Fig. m one has to use the relations n, a —> n, 2, n,b n — 1,3, 
and n,c n, 1. 

h Again we use here the notations (indices) of Ref. 1591 To make the relationship to the notations 
(indices) used throughout our paper (see Fig. [TJ one has to use the relations n,a —»■ n, 1 and 
n,b — > n, 2 . 
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By using the transformation 

1 


Pn — 


fn — ^ (&n b n), 

_ 1 _ 1 ( 

^n,+ — 2 v^ n i a "h ^n,— — 2 \^ n 5 a £ n,b) 


(83) 


one obtains a new representation of the initial problem that correspon ds t o a lattice 
with dispersive degrees of freedom p n and side-coupled Fano states /„•*“ 


Ep n = (e„,+ - t 2 )p n + €n,-fn ~ 2(p„+l + Pra-l), 

Ef n ~ (^ra,+ "b ^2)/ra “b €n,—Pn- (84) 

In Ref. [59] such a lattice with side-coupling defects is called “Fano lattice”. This 
procedure can be inverted: Choose a dispersive chain, add a set of uncoupled Fano 
states f n with energies Ef n , assign locally each /„ to a site with p n , and perform 
local rotations (transformations) in the space {p n , /„}. As a result, one obtains the 
original frustrated two-leg ladder lattice. This construction can be generalized and 
it yields complex lattices having flat-band localized states. 

The above described Fano-lattice representation of the frustrated two-leg ladder 
allows to discuss the impact of perturbations of the initial model from a completely 
different viewpoint. We consider a site-dependent on-site potential e n . For e n - = 
5 (era,a — £n,b) = 0 the Fano states will remain decoupled, but their degeneracy is 
lifted, since e„ j+ yb 0, see Eq. (l84]l . A particular case corresponds to uncorrelated 
random numbers e rajQ = e nj {, = e n with uniformly distributed random variables e n € 
[— f’ T 1 ] • The Fano states stay decoupled, however, acquire an energy spread of the 
order of W around Efb = t 2 - The dispersive lattice becomes Anderson localized 
with a localization length £ ~ W ~ 2 for weak disorder W < 4. If the symmetry 
constraint is relaxed, i.e., e n ,_ yb 0, the Fano states are locally coupled to the 
dispersive chain. Because of the pure local coupling, they can be eliminated re sultin g 
in a equation for the dispersive lattice with a Lorentzian (Cauchy) disordeJ^^Sl jf 
\E — Efb\ < If \E — Efb\ > the Lorentzian disorder does not appear and the 
correlation-length exponent is 7 = 2. These analytical arguments are in excellent 
agreement with numerical findings!^ However, at Efb = 0 the localization-length 
exponent is found to be 7 = 1.30±0.01 for the diamond chain, by contrast to 7 = 1 as 
expected for the Lorentzian distribution. On the other hand, for other models, e.g., 
for the frustrated two-leg ladder (or “cross-stitch lattice”), the expected exponent 
7 = 1 is found. This unexpected finding for the diamond chain remains a puzzle to 
be explained. 


11.3. Conductance through flat-band clusters 

The conductance through various nanosystems such as molecular devices or 
nanowires has been extensively studied both theoretically and experimentally in 
recent years. The main nanotransport phenomena (e.g., Coulomb blockade, con¬ 
ductance quantization, resonant tunneling etc.) are meanwhile well studied and 
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- 0 " 

Fig. 20. (Color online) Periodic A ,r -ceII diamond chain (A f = 4) connected to semi-infinite leads 
via hoping amplitudes tj y and tj{ (dashed bonds). 



understood. For nanoclusters with flat-band geometry the existence of localized 
states may lead to distinct features, i.e., localiz ed states may manifest themselves 
in an intriguing way in the conductance! ^ 

We follow Ref. [60] and consider as an example a finite periodic diamond chain 
(i.e., a diamond ring) without frustration (the vertical bond is zero, cf. Fig. [l]) which 
is connected to left and right semi-infinite leads, see Fig. [20] The Hamiltonian of 
the two-terminal device is as follows: 


H — H leads + H cluster A H‘ ln t , 

0 oo 

Pleads =-t^2 (\j - l)(j| +H.C.) -t ^ (li)O' -*-!!+ H.C.) , 

j=—oo j=N-\-l 

AT 

^cluster = -t y, \e 1 ^ (Im, 3)(m, 1| + |m, l)(m - 1,3|) 
m= 1 

+e _i 2 jv (| m , 2 )(m, 3| + |m — 1, 3)(m, 2|) + H.c. , 

Hint = -t L (|0)(1| + H.c.) - t R (\N)(N + 1| + H.c.). (85) 

The phase factors in H c i us ter are related to fluxes, where the fluxes $ and <fq are 
threading the plaquettes and the inner ring, respectively, and 4 > 0 = <E>i + 0 $. 
Furthermore, | j) = cj|vac) (respectively |m,a) = |vac), a = 1,2,3) is a state, 
where the fermion is located at the lattice site j of the lead (at the lattice site m, a 
of the diamond ring). The pairs of indices labeling the sites of the diamond ring 
enumerate the cells m = 1 ,... ,Af and the position a = 1,2,3 of the site within 
a cell, see Fig. [[] The connecting hopping amplitudes t R and t R are assumed to 
be small compared to the hopping integrals t. The site 1 ( N ) stands for the left 
(right) contact site of the cluster, cf. the left and right panel in Fig. [2D] The spin 
variable a =t, 4 - is omitted, since it is irrelevant in this tight-binding treatment, i.e., 
for U = 0. 

Within the Landauer picturepMSH which treats the electric current as a con¬ 
sequence of the injection of a carrier at a contact and its probability to reach the 
other end, the conductance is determined by the formula 

G = G 0 \t(E)\ 2 , 


( 86 ) 
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where Go = ^- is the conductance quantum and |t(E )| 2 is the transmission prob¬ 


ability of a propagating electron with the energy E. For the two-terminal device 
the propagating electron comes from the left lead with the energy E = — 2fcosK, 
—7r < k < 7 T. Going along standard reasoning of quantum scattering theory, one 
seeks for a solution of the stationary Schrodinger equation, H\ijj) = E\i/j), where 
l^) = V'nl 71 ) (n = J or n = m, a). In the left and right lead the amplitudes have 

the form: 



(87) 


i.e., it is a linear combination of the incoming and the reflected waves in the left 
lead and the transmitted wave in the right lead. There are 3A f + 2 equations for 
3A f + 2 unknown quantities r, t and ipi, where i = m, a runs over the sites 
of the diamond ring. Solving these equations we determine t which is, in fact, the 
required amplitude of the transmission probability, t = t (E). [Note that instead of 
the standard approach to calculate t (E) via Eq. (IHTjl the auth ors of Ref. [60] use a 
Green’s function representation of the scattering problem. 

Let us briefly describe some interesting findings reported in Ref. 1601 We begin 
with the case of zero flux. If the leads are connected to a site to, 3 of the diamond 
ring (left panel in Fig. [50]) one gets an expected dependence of the conductance 
on the gate voltage Igate- It exhibits three Breit-Wigner-shaped peaks at energies 
—2 1, 0, 21 . The latter values correspond to the eigenvalues of the cluster which 
follows from the band energies, { 0 , ±2ty/l + cos k} (the flat band has zero energy), 
after substitution n = —7 r, — 0 , ^. By contrast to this usual conductance behavior, 
the situation is changed, if the left lead is connected to a site to, 1 (right panel in 
Fig. l20l) : Then the E = 0-peak in the conductance profile is absent, i.e., there 
is no conductance, if E corresponds to the flat-band value. The authors of Ref. [60 
explain this as follows. The hopping term to a site to, 1 of the diamond ring describes 
a hopping to an itinerant state of the ring as well as to a localized state. Since the 
localized state is decoupled from all other states of the ring, it leads only to a 
reflected wave back into the left lead. For E = 0 this reflected wave interferes 
destructively with the incident wave resulting in the absence of the zero-energy 
peak in the conductance. 

In the presence of finite fluxes the behavior at E = 0 may be different in case of 
the connection of the left lead to a site to, 1 and the right lead to site to, 2 , to — 1 , 2 
etc. Due to fluxes a gap between the flat-band states with zero energy and itinerant 
states opens. The conductance behavior for E in a small region around E = 0 is 
again unusual: The conductance is still zero precisely at E = 0, but in close vicinity 
on both sides of E = 0 two well-pronounced peaks emerge, i.e., the conductance 
exhibits a double peak symmetric to zero energy. This is unexpected, since the 
energies of the itinerant states of the ring are far from this zero-energy peak region 
and this is related to variable extensions of localized states in the presence of fluxes. 
For further details we refer the interested reader to Ref. [601 
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In summary, localized states in a flat-band tight-binding cluster generate a strik¬ 
ing behavior in the two-terminal electronic conductance. The positions of the con¬ 
tact sites are extremely important. The localized states may act as a filter of the 
zero-energy conductance peak in absence of a flux or may lead to the appearance 
of a conductance double peak around zero energy. Although the tight-binding cal¬ 
culations described above do not take into account the interelectron interactions, 
one may expect that the peculiarities will survive when interactions are present, see 
Refs. [HHl 11481IT491IT501 

12. Experiments 

The existence of ferromagnetism in a basic model of strongly correlated electrons, 
the one-band repulsive Hubbard model, requires a particular pattern of hopping 
integrals, i.e., a special lattice topology. The search for this kind of ferromagnetism 
is an attractive goal for many researches and it is an ambitious task to find its 
traces in specific materials. In particular, flat-band ferromagnetism opens the win¬ 
dow to ferromagnets consisting entirely of nonmagnetic elements such as organic 
molecules. There are various directions of ongoing experimental search for flat-band 
ferromagnetism. 

The realization of flat-band ferromagnetism in quantum dot arrays and in quan¬ 
tum atomic wires formed on solid surfaces was discussed in Ref. IED While in real 
materials a lattice distortion may destabilize the ferromagnetism when the flat band 
is half filled, artificial materials such as dot lattices do not have such disadvan¬ 
tages. H. Tamura et a,b-^ demonstrated that flat-band ferromagnetism should be 
observable at reasonable temperatures in semiconductor dot arrays using existing 
fabrication technology. They performed calculations for GaAs, InAs, and Si dot 
arrays of Lieb- and kagome-lattice structures of typical interdot spacing of 5, 10, 
50, and 100 nm and dot diameter two times smaller than the interdot spacing. The 
designed structures, however, may require an ultra-low-temperature technique to 
detect flat-band ferromagnetism, see the paper by K. Shiraishi et al. in Ref. [6T . 

Other examples for flat-band ferromagnets are related to polymers,®^ where the 
sear ch fo r purely organic ferromagnets is known as a challenging target. R. Arita 
et al.®^ found that a chain of five-membered rings (polyaminotriazole) should be 
ferromagnetic with an appropriate doping. Furthermore, the authors demonstrated 
that the predicted flat-band ferromagnetism of a single chain may occur as a bulk 
property, i.e., in a three-dimensional crystal of the polymer. 

Another experimental realization of flat-band ferromagnetism has been dis¬ 
cussed recently in Ref. EU Specific organic molecules [tetracyano-p-quinodimethane 
(TCNQ) molecules] deposited on graphene epitaxially grown on Ru(0001) acquire 
a charge from the substrate and develop a magnetic moment which survives when 
the molecules form a monolayer. The authors of Ref. [65] claim that the intermolec- 
ular bands are almost flat and half-filled and that the TCNQ/graphene/Ru(0001) 
system might be a realization of a flat-band Hubbard ferromagnet. 
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Just recently, the el ectro nic structure of silicene has been examined and flat 
bands have been found P-S-E The authors claim that if we can shift the chemical 
potential close to the flat band (e.g., by chemical doping) flat-band ferromagnetism 
may be expected. 

Although, so far no unambiguous example of an experimental realization of a 
flat-band ferromagnet was found this field of experimental research is very promising 
not only for fundamental physics but also for applied physics (having, e.g., in mind 
the possible realization of organic ferromagnets). Within this context it seems to 
be important that (small) deviations from the ideal flat-band geometry are not 
necessarily contraindicative for the realization of flat-band ferromagnetism. 

13. Summary and Outlook 

In this review we offer a unifying view on spin and electron systems possessing flat 
bands. The existence of flat bands implies similarities in mathematical description, 
although the statistics of particles and the interactions lead to differences in proper¬ 
ties. A study of strongly correlated flat-band systems leads to some surprises. First 
of all, and contrary to what one may expect, it appears that in spite of compli¬ 
cated (often frustrated) lattices which are usually required to get flat bands, the 
theoretical description may be simpler and it involves specific methods of classi¬ 
cal statistical mechanics. Second, there are a number of intriguing properties of 
strongly correlated flat-band systems with important potential applications. For 
example, flat-band ferromagnetism offers a route to ferromagnets without magnetic 
atoms. Other promising examples are the application for magnetic cooling using the 
magnetocaloric effect in frustrated quantum antiferromagnets or for the control of 
the charge transport through flat-band nanodevices. 

We mention again that flat-band physics is a much broader field than it was 
discussed in this review; examples are topol ogica l flat-band models,^ superconduc¬ 
tivity/superfluidity in flat-band systems 11521153 dispersion 0 f phononsflat 
bands in an extended Falicov-Kimball model fLSH flat-band-like states in f racta l 
networksJISSI “flat-band” (i.e., localize d-stat e) effects in magnetic moleculei^EH 
or nonequilibrium flat-band phenomenaP-^ From a general point of view it might 
be interesting to examine systematically other flat-band systems and to find more 
traces of flat-band physics in experiments. We hope this review will encourage fur¬ 
ther exploration of strongly correlated flat-band systems along the described lines. 

Acknowledgements 

We acknowledge useful discussions and related collaborations with S. Capponi, D. 
V. Dmitriev, S.-L. Drechsler, A. Honecker, V. Ya. Krivnov, T. Krokhmalskii, O. 
Krupnitska, A. M. Lauchli, R. Moessner, S. Nishinroto, and K. Shtengel. We spe¬ 
cially thank D. V. Dmitriev, S. Flach, Z. Gulacsi, V. Ya. Krivnov, A. A. Lopes, and 
H. Tanaka for valuable comments on the manuscript. We thank A. Andrusyk and 
T. Verkholyak for critical reading of the manuscript. 




















April 14, 2015 0:38 WSPC/INSTRUCTION FILE derimarevised 


Strongly correlated flat-band systems: The route from Heisenberg spins to Hubbard electrons 65 

The present study was supported by the DFG (project RI615/21-1). O. D. ac¬ 
knowledges financial support of the Abdus Salam International Centre for Theo¬ 
retical Physics (Trieste) through the Senior Associate award. O. D. and J. R. are 
grateful to the MPIPKS (Dresden) for the kind hospitality in October-December of 
2013. O. D. is grateful to the University of Magdeburg for the kind hospitality in 
April-June and October-December of 2014. O. D. acknowledges kind hospitality of 
the organizers of the 567. WE-Heraeus-Seminar on “Integrablc Lattice Models and 
Quantum Field Theories” (Bad Honnef, Germany, June 28 - July 2, 2014) and of 
the Conference on Field Theory Methods in Low-Dimensional Strongly Correlated 
Quantum Systems (ICTP, Trieste, Italy, August 25-29, 2014). 

References 

1. S. A. Parameswaran, R. Roy and S. L. Sondhi, C. R. Physique 14, 816 (2013). 

2. E. J. Bergholtz and Z. Liu, Int. J. Mod. Phys. B27, 1330017 (2013). 

3. Z. Liu, F. Liu and Y.-S. Wu, Chinese Phys. B23, 077308 (2014). 

4. C. Zeng and V. Elser, Phys. Rev. B42, 8436 (1990). 

5. H. Asakawa and M. Suzuki, Physica A205, 687 (1994). 

6. A. P. Ramirez, Annu. Rev. Mater. Sci. 24, 453 (1994). 

7. J. E. Greedan, J. Mater. Chem. 11, 37 (2001). 

8. J. Schnack, H.-.I. Schmidt, J. Richter and J. Schulenburg, Eur. Phys. J. B24, 475 
(2001); J. Schulenburg, A. Honecker, J. Schnack, J. Richter and H.-J. Schmidt, Phys. 
Rev. Lett. 88, 167207 (2002); J. Richter, J. Schulenburg, A. Honecker, J. Schnack 
and H.-J. Schmidt, J. Phys.: Condens. Matter 16, S779 (2004). 

9. J. Richter, O. Derzhko and J. Schulenburg, Phys. Rev. Lett. 93, 107206 (2004); O. 
Derzhko, J. Richter and J. Schulenburg, physica status solidi (b) 242, 3189 (2005); 
O. Derzhko and J. Richter, Phys. Rev. B72, 094437 (2005). 

10. O. Derzhko and J. Richter, Phys. Rev. B70, 104415 (2004). 

11. M. E. Zhitomirsky and A. Honecker, J. Stat. Mech.: Theor. Exp. P07012 (2004). 

12. M. E. Zhitomirsky and H. Tsunetsugu, Phys. Rev. B70, 100403(R) (2004); M. E. 
Zhitomirsky and H. Tsunetsugu, Prog. Theor. Phys. Suppl. 160, 361 (2005); M. E. 
Zhitomirsky and H. Tsunetsugu, Phys. Rev. B75, 224416 (2007). 

13. O. Derzhko and J. Richter, Eur. Phys. J. B52, 23 (2006); J. Richter, O. Derzhko 
and T. Krokhmalskii, Phys. Rev. B74, 144430 (2006). 

14. A. Honecker and J. Richter, Condensed Matter Physics (L’viv) 8, 813 (2005); J. 
Richter, O. Derzhko and A. Honecker, Int. J. Mod. Phys. B22, 4418 (2008); O. 
Derzhko, J. Richter and A. Honecker, J. Phys.: Conf. Ser. 145, 012059 (2009). 

15. A. Mielke, J. Phys. A24, L73 (1991); A. Mielke, J. Phys. A24, 3311 (1991); A. 
Mielke, J. Phys. A25, 4335 (1992); A. Mielke, Phys. Lett. A174, 443 (1993). 

16. H. Tasaki, Phys. Rev. Lett. 69, 1608 (1992). 

17. A. Mielke and H. Tasaki, Commun. Math. Phys. 158, 341 (1993). 

18. H. Tasaki, Phys. Rev. Lett. 73, 1158 (1994); H. Tasaki, Phys. Rev. Lett. 75, 4678 
(1995); H. Tasaki, J. Stat. Phys. 84, 535 (1996); H. Tasaki, Commun. Math. Phys. 
242, 445 (2003); see also K. Kusakabe and H. Aoki, Phys. Rev. Lett. 72, 144 (1994); 
K. Kusakabe and H. Aoki, Physica B194-196, 215 (1994); H. Aoki, Int. J. Mod. 
Phys. B17, 4953 (2003). 

19. O. Derzhko, A. Honecker and J. Richter, Phys. Rev. B76, 220402(R) (2007); A. 
Honecker, O. Derzhko and J. Richter, Physica B404, 3316 (2009). 

20. O. Derzhko, A. Honecker and J. Richter, Phys. Rev. B79, 054403 (2009). 








April 14, 2015 0:38 WSPC/INSTRUCTION FILE de'rimarevised 


66 O. Derzhko, J. Richter, M. Maksymenko 

21. O. Derzhko, J. Richter, A. Honecker, M. Maksymenko and R. Moessner, Phys. Rev. 
B81, 014421 (2010); O. Derzhko, M. Maksymenko, J. Richter, A. Honecker and R. 
Moessner, Acta Physica Polonica A118, 736 (2010). 

22. M. Maksymenko, A. Honecker, R. Moessner, J. Richter and O. Derzhko, “The flat- 
band ferromagnetic transition as a Pauli-correlated percolation transition”, preprint 
of the Institute for Condensed Matter Physics of the National Academy of Sciences 
of Ukraine (L’viv), ICMP-11-14E (2011). 

23. M. Maksymenko, A. Honecker, R. Moessner, J. Richter and O. Derzhko, Phys. Rev. 
Lett. 109, 096404 (2012). 

24. M. Maksymenko, Ph.D. Thesis, Institute for Condensed Matter Physics, National 
Academy of Sciences of Ukraine, L’viv, 2012. 

25. M. Maksymenko, R. Moessner and K. Shtengel, “Reversible first-order transition in 
Pauli percolation”, arXiv: 1401.6172 

26. C. Lhuillier and G. Misguich, in High Magnetic Fields: Applications in Condensed 
Matter Physics and Spectroscopy, Lecture Notes in Physics Vol. 595, ed. C. Berthier, 
L. P. Levy and G. Martinez (Springer-Verlag, Berlin, 2002), pp. 161-190; G. Mis¬ 
guich and C. Lhuillier, in Frustrated Spin Systems, ed. H. T. Diep (World Scientific, 
Singapore, 2013), pp. 235-319. 

27. .J. Richter, J. Schulenburg and A. Honecker, in Quantum Magnetism, Lecture Notes 
in Physics Vol. 645, ed. U. Schollwock, J. Richter, D. J. J. Farnell and R. F. Bishop 
(Springer-Verlag, Berlin, 2004), pp. 85-153. 

28. Introduction to Frustrated Magnetism. Materials, Experiments, Theory, ed. C. 
Lacroix, P. Mendels and F. Mila (Springer-Verlag, Berlin, Heidelberg, 2011). 

29. M. C. Gutzwiller, Phys. Rev. Lett. 10, 59 (1963); J. Hubbard, Proc. Roy. Soc. London 
A276, 238 (1963); J. Kanamori, Prog. Theor. Phys. 30, 275 (1963). 

30. Editorial, Nature Physics 9, 523 (2013). 

31. H. Tasaki, J. Phys.: Condens. Matter 10, 4353 (1998). 

32. H. Tasaki, Prog. Theor. Phys. 99, 489 (1998). 

33. H. Tasaki, Eur. Phys. J. B64, 365 (2008). 

34. J. Richter, Fizika Nizkikh Temperatur (Kharkiv) 31, 918 (2005) [Low Temperature 
Physics 31, 695 (2005)]. 

35. O. Derzhko, J. Richter, A. Honecker and H.-J. Schmidt, Fizika Nizkikh Temperatur 
(Kharkiv) 33, 982 (2007) [Low Temperature Physics 33, 745 (2007)]. 

36. J. Richter and O. Derzhko, in Condensed Matter Physics in the Prime of the 21st 
Century. Phenomena, Materials, Ideas, Methods, ed. J. Jsgdrzejewski (World Scien¬ 
tific, Singapore, 2008), pp. 237-270. 

37. O. Derzhko, T. Krokhmalskii and J. Richter, Phys. Rev. B82, 214412 (2010); O. 
Derzhko, T. Krokhmalskii and J. Richter, Teoret. Mat. Fiz. 168, 441 (2011) [Theo¬ 
retical and Mathematical Physics 168, 1236 (2011)]. 

38. M. Maksymenko, O. Derzhko and J. Richter, Acta Physica Polonica A119, 860 
(2011); M. Maksymenko, O. Derzhko and .]. Richter, Eur. Phys. J. B84, 397 (2011). 

39. A. Mielke, Eur. Phys. J. B85, 184 (2012). 

40. O. Derzhko, J. Richter and O. Krupnitska, Condensed Matter Physics (L’viv) 15, 
43702 (2012). 

41. O. Derzhko, J. Richter, O. Krupnitska and T. Krokhmalskii, Phys. Rev. B88, 094426 
(2013); O. Derzhko, J. Richter, O. Krupnitska and T. Krokhmalskii, Fizika Nizkikh 
Temperatur (Kharkiv) 40, 662 (2014) [Low Temperature Physics 40, 513 (2014)]. 

42. S. D. Huber and E. Altman, Phys. Rev. B82, 184502 (2010). 

43. S. Nishimoto, N. Shibata and C. Hotta, Nat. Commun. 4, 2287 (2013). 

44. S. Capponi, O. Derzhko, A. Honecker, A. M. Lauchli and J. Richter, Phys. Rev. B88, 








April 14, 2015 0:38 WSPC/INSTRUCTION FILE derimarevised 


Strongly correlated flat-band systems: The route from Heisenberg spins to Hubbard electrons 67 


144416 (2013). 

45. H. Nakano and T. Sakai, J. Phys. Soc. Jpn. 79, 053707 (2010); T. Sakai and H. 
Nakano, Phys. Rev. B83, 100405(R) (2011); T. Sakai and H. Nakano, J. Phys.: 
Conf. Ser. 320, 012016 (2011); T. Sakai and H. Nakano, physica status solidi (b) 
250, 579 (2013). 

46. V. Ya. Krivnov, D. V. Dmitriev, S. Nishimoto, S.-L. Drechsler and J. Richter, Phys. 
Rev. B90, 014441 (2014). 

47. H. Kikuchi, Y. Fujii, M. Chiba, S. Mitsudo, T. Idehara, T. Tonegawa, K. Okamoto, 
T. Sakai, T. Kuwai and H. Ohta, Phys. Rev. Lett. 94, 227201 (2005); H. Kikuchi, 
Y. Fujii, M. Chiba, S. Mitsudo, T. Idehara, T. Tonegawa, K. Okamoto, T. Sakai, 
T. Kuwai, K. Kindo, A. Matsuo, W. Higemoto, K. Nishiyama, M. Horvatic and C. 
Bertheir, Prog. Theor. Phys. Suppl. 159, 1 (2005). 

48. H. Kikuchi, Y. Fujii, M. Chiba, S. Mitsudo, T. Idehara, T. Tonegawa, K. Okamoto, 
T. Sakai, T. Kuwai and H. Ohta, Phys. Rev. Lett. 97, 089702 (2006). 

49. K. C. Rule, A. U. B. Wolter, S. Siillow, D. A. Tennant, A. Briihl, S. Kohler, B. Wolf, 
M. Lang and J. Schreuer, Phys. Rev. Lett. 100, 117202 (2008). 

50. H. Jeschke, I. Opahle, H. Kandpal, R. Valenti, H. Das, T. Saha-Dasgupta, O. Janson, 
H. Rosner, A. Briihl, B. Wolf, M. Lang, J. Richter, S. Hu, X. Wang, R. Peters, T. 
Pruschke and A. Honecker, Phys. Rev. Lett. 106, 217201 (2011). 

51. A. Honecker, S. Hu, R. Peters and J. Richter, J. Phys.: Condens. Matter 23, 164211 
( 2011 ). 

52. P. T. Cong, B. Wolf, R. S. Manna, U. Tutsch, M. de Souza, A. Briihl and M. Lang, 
Phys. Rev. B89, 174427 (2014). 

53. H. Tanaka, N. Kurita, M. Okada, E. Kunihiro, Y. Shirata, K. Fujii, H. Uekusa, A. 
Matsuo, K. Kindo and H. Nojiri, J. Phys. Soc. Jpn. 83, 103701 (2014). 

54. O. Derzhko and J. Richter, Phys. Rev. B90, 045152 (2014). 

55. Z. Gulacsi, A. Kampf and D. Vollhardt, Phys. Rev. Lett. 99, 026404 (2007); Z. 
Gulacsi, A. Kampf and D. Vollhardt, Prog. Theor. Phys. Suppl. 176, 1 (2008). 

56. Z. Gulacsi, A. Kampf and D. Vollhardt, Phys. Rev. Lett. 105, 266403 (2010); Z. 
Gulacsi, Int. J. Mod. Phys. B27, 1330009 (2013); M. Gulacsi, G. Kovacs and Z. 
Gulacsi, Phil. Mag. Lett. 94, 269 (2014); M. Gulacsi, G. Kovacs and Z. Gulacsi, 
Europhys. Lett. 107, 57005 (2014). 

57. M. Goda, S. Nishino and H. Matsuda, Phys. Rev. Lett. 96, 126401 (2006); S. Nishino, 
H. Matsuda and M. Goda, J. Phys. Soc. Jpn. 76, 024709 (2007). 

58. J. T. Chalker, T. S. Pickles and P. Shukla, Phys. Rev. B82, 104209 (2010). 

59. D. Leykam, S. Flach, O. Bahat-Treidel and A. S. Desyatnikov, Phys. Rev. B88, 
224203 (2013); S. Flach, D. Leykam, J. D. Bodyfelt, P. Matthies and A. S. Desy¬ 
atnikov, Europhys. Lett. 105, 30001 (2014) and 106, 19901 (2014); J. D. Bodyfelt, 
D. Leykam, C. Danieli, X. Yu and S. Flach, Phys. Rev. Lett. 113, 236403 (2014); 
C. Danieli, J. D. Bodyfelt and S. Flach, “Flatband engineering of mobility edges”, 
arXiv: 1502.06690 

60. A. A. Lopes, A. A. Z. Antonio and R. G. Dias, Phys. Rev. B89, 235418 (2014); A. A. 
Lopes and R. G. Dias, “Simple approach for the two-terminal conductance through 
interacting clusters”, arXiv: 1407.5922 see also A. A. Lopes and R. G. Dias, Phys. 
Rev. B84, 085124 (2011). 

61. H. Tamura, K. Shiraishi, T. Kimura and H. Takayanagi, Phys. Rev. B65, 085324 
(2002); T. Kimura, H. Tamura, K. Shiraishi and H. Takayanagi, Phys. Rev. B65, 
081307 (2002); K. Shiraishi, H. Tamura and H. Takayanagi, Physica E24, 107 (2004); 
M. Ichimura, K. Kusakabe, S. Watanabe and T. Onogi, Phys. Rev. B58, 9595 (1998); 
H. Ishii, T. Nakayama and J.-i. Inoue, Phys. Rev. B69, 085325 (2004). 








April 14, 2015 0:38 WSPC/INSTRUCTION FILE de'rimarevised 


68 O. Derzhko, J. Richter, M. Maksymenko 

62. S. Nishino, M. Goda and K. Kusakabe, J. Phys. Soc. Jpn. 72, 2015 (2003); S. Nishino 
and M. Goda, J. Phys. Soc. Jpn. 74, 393 (2005). 

63. R. Arita, Y. Suwa, K. Kuroki and H. Aoki, Phys. Rev. Lett. 88 , 127202 (2002); 

R. Arita, Y. Suwa, K. Kuroki and H. Aoki, Phys. Rev. B68, 140403(R) (2003); 

Y. Suwa, R. Arita, K. Kuroki and H. Aoki, Phys. Rev. B68, 174419 (2003); H. Aoki, 
Appl. Surf. Sci. 237, 2 (2004). 

64. H.-H. Lin, T. Hikihara, H.-T. Jeng, B.-L. Huang, C.-Y. Mou and X. Hu, Phys. Rev. 
B79, 035405 (2009). 

65. M. Garnica, D. Stradi, S. Barja, F. Calleja, C. Diaz, M. Alcami, N. Martin, A.L. 
Vazquez de Parga, F. Martin and R. Miranda, Nat. Phys. 9, 368 (2013). 

66. C. Wu, D. Bergman, L. Balents and S. Das Sarma, Phys. Rev. Lett. 99, 070401 
(2007); D. L. Bergman, C. Wu and L. Balents, Phys. Rev. B78, 125104 (2008); M. 
Zhang, H.-h. Huang, C. Zhang and C. Wu, Phys. Rev. A83, 023615 (2011). 

67. V. Apaja, M. Hyrkas and M. Manninen, Phys. Rev. A82, 041402(R) (2010). 

68. M. Hyrkas, V. Apaja and M. Manninen, Phys. Rev. A87, 023614 (2013). 

69. K. Noda, K. Inaba and M. Yamashita, “Flat-band ferromagnetism in the multilayer 
Lieb optical lattice”, arXiv: 1410.2166 

70. S. Mukherjee, A. Spracklen, D. Choudhury, N. Goldman, P. Ohberg, E. Andersson 
and R. R. Thomson, “Observation of a localized flat-band state in a photonic Lieb 
lattice”, arXiv: 1412.6342 

71. M. Powalski, K. Coester, R. Moessner and K. P. Schmidt, Phys. Rev. B87, 054404 
(2013); K. Coester, W. Malitz, S. Fey and K. P. Schmidt, Phys. Rev. B88, 184402 
(2013). 

72. T. Nakamura and K. Kubo, Phys. Rev. B53, 6393 (1996); D. Sen, B. S. Shastry, R. 

E. Walstedt and R. Cava, Phys. Rev. B53, 6401 (1996). 

73. Ch. Waldtmann, H. Kreutzmann, U. Schollwock, K. Maisinger and H.-U. Everts, 
Phys. Rev. B62, 9472 (2000); P. Azaria, C. Hooley, P. Lecheminant and A. M. Tsve- 
lik, Phys. Rev. Lett. 81, 1694 (1998). 

74. K. Takano, K. Kubo and H. Sakamoto, J. Phys.: Condens. Matter 8, 6405 (1996); 
H. Niggemann, G. Uimin and J. Zittartz, J. Phys.: Condens. Matter 9, 9031 (1997); 
H. Niggemann, G. Uimin and J. Zittartz, J. Phys.: Condens. Matter 10 , 5217 (1998); 
K. Okamoto, T. Tonegawa and M. Kaburagi, J. Phys.: Condens. Matter 15, 5979 
(2003); L. Canova, J. Strecka and M. Jascur, J. Phys.: Condens. Matter 18, 4967 
(2006). 

75. M. P. Gelfand, Phys. Rev. B43, 8644 (1991); F. Mila, Eur. Phys. J. B6, 201 (1998); 
A. Honecker, F. Mila and M. Troyer, Eur. Phys. J. B15, 227 (2000). 

76. M. Mambrini, J. Trebosc and F. Mila, Phys. Rev. B59, 13806 (1999); O. Rojas and 

F. C. Alcaraz, Phys. Rev. B67, 174401 (2003); C. D. Batista and B. S. Shastry, Phys. 
Rev. Lett. 91, 116401 (2003); D. Antonosyan, S. Bellucci and V. Ohanyan, Phys. Rev. 
B79, 014432 (2009). 

77. V. Subrahmanyam, Phys. Rev. B50, 16109 (1994); A. Liischer, R. M. Noack, G. 
Misguich, V. N. Kotov and F. Mila, Phys. Rev. B70, 060405(R) (2004); G. Seeber, P. 
Kogerler, B. M. Kariuki and L. Cronin, Chem. Commun., 1580 (2004); J. Schnack, 
H. Nojiri, P. Kogerler, G. J. T. Cooper and L. Cronin, Phys. Rev. B70, 174420 
(2004); J.-B. Fouet, A. Lauchli, S. Pilgram, R. M. Noack and F. Mila, Phys. Rev. 
B73, 014409 (2006); N. B. Ivanov, J. Schnack, R. Schnalle, J. Richter, P. Kogerler, 

G. N. Newton, L. Cronin, Y. Oshima and H. Nojiri, Phys. Rev. Lett. 105, 037206 
(2010); T. Sakai, M. Sato, K. Okamoto, K. Okunishi and C. Itoi, J. Phys.: Condens. 
Matter 22, 403201 (2010); M. Lajko, P. Sindzingre and K. Penc, Phys. Rev. Lett. 
108, 017205 (2012). 








April 14, 2015 0:38 WSPC/INSTRUCTION FILE de'rimarevised 


Strongly correlated flat-band systems: The route from Heisenberg spins to Hubbard electrons 69 

78. P. Lecheminant, B. Bernu, C. Lhuillier, L. Pierre and P. Sindzingre, Phys. Rev. B56, 
2521 (1997); Ch. Waldtmann, H.-U. Everts, B. Bernu, C. Lhuillier, P. Sindzingre, P. 
Lecheminant and L. Pierre, Eur. Phys. J. B2, 501 (1998); F. Mila, Phys. Rev. Lett. 
81, 2356 (1998). 

79. R. R. P. Singh and D. A. Huse, Phys. Rev. B76, 180407(R) (2007); S. Yan, D. A. 
Huse and S. R. White, Science 332, 1173 (2011); S. Depenbrock, I. P. McCulloch 
and U. Schollwock, Phys. Rev. Lett. 109, 067201 (2012); H. Nakano and T. Sakai, J. 
Phys. Soc. Jpn. 80, 053704 (2011); A. M. Lauchli, J. Sudan and E. S. Sprensen, Phys. 
Rev. B83, 212401 (2011); O. Gotze, D. J. J. Farnell, R. F. Bishop, P. H. Y. Li and 
J. Richter, Phys. Rev. B84, 224428 (2011); Z. Y. Xie, J. Chen, J. F. Yu, X. Kong, 
B. Normand and T. Xiang, Phys. Rev. X4, 011025 (2014); A. L. Chernyshev and M. 

E. Zhitomirsky, Phys. Rev. Lett. 113, 237202 (2014); O. Gotze and J. Richter, Phys. 
Rev. B91, 104402 (2015). 

80. R. Siddharthan and A. Georges, Phys. Rev. B65, 014417 (2002); P. Tomczak and 
J. Richter, J. Phys. A36, 5399 (2003); J. Richter, J. Schulenburg, P. Tomczak and 
D. Schmalfufi, Condensed Matter Physics (L’viv) 12, 507 (2009); H. Nakano and T. 
Sakai, J. Phys. Soc. Jpn. 82, 083709 (2013). 

81. I. Rousochatzakis, R. Moessner and J. van den Brink, Phys. Rev. B88, 195109 (2013). 

82. H.-Q. Lin and J. L. Shen, J. Phys. Soc. Jpn. 69, 878 (2000); P. Chen, C.-Y. Lai and 
M.-F. Yang, Phys. Rev. B81, 020409(R) (2010); A. F. Albuquerque, N. Laflorencie, 
J.-D. Picon and F. Mila, Phys. Rev. B83, 174421 (2011); Y. Murakami, T. Oka and 
H. Aoki, Phys. Rev. B88, 224404 (2013). 

83. N. B. Ivanov and J. Richter, Phys. Lett. A232, 308 (1997); J. Richter, N. B. Ivanov 
and J. Schulenburg, J. Phys.: Condens. Matter 10, 3635 (1998); A. Koga, K. Okunishi 
and N. Kawakami, Phys. Rev. B62, 5558 (2000); J. Schulenburg and J. Richter, Phys. 
Rev. B65, 054420 (2002); A. Koga and N. Kawakami, Phys. Rev. B65, 214415 (2002). 

84. J. Richter, J. Schulenburg, A. Honecker and D. Schmalfufi, Phys. Rev. B70, 174454 
(2004), T.-P. Choy and Y. B. Kim, Phys. Rev. B80, 064404 (2009); B.-J. Yang, A. 
Paramekanti and Y. B. Kim, Phys. Rev. B81, 134418 (2010). 

85. S. E. Palmer and J. T. Chalker, Phys. Rev. B64, 094412 (2001); W. Brenig and A. 
Honecker, Phys. Rev. B65, 140407 (2002); J.-B. Fouet, M. Mambrini, P. Sindzingre 
and C. Lhuillier, Phys. Rev. B67, 054411 (2003). 

86. B. Canals and C. Lacroix, Phys. Rev. Lett. 80, 2933 (1998). 

87. J. Strecka, L. Canova, M. Jascur and M. Hagiwara, Phys. Rev. B78, 024427 (2008); 
D.-X. Yao, Y. L. Loh, E. W. Carlson and M. Ma, Phys. Rev. B78, 024428 (2008); J. 
Cisarova and J. Strecka, Phys. Rev. B87, 024421 (2013); J. Cisarova, F. Michaud, 

F. Mila and J. Strecka, Phys. Rev. B87, 054419 (2013). 

88. J. Richter, R. Schmidt and J. Schnack, J. Magn. Magn. Mat. 295, 164 (2005); 
J. Schnack, R. Schmidt and J. Richter, Phys. Rev. B76, 054413 (2007); I. 
Rousochatzakis, A. M. Lauchli and F. Mila Phys. Rev. B77, 094420 (2008); J. 
Schnack and R. Schnalle, Polyhedron 28, 1620 (2009). 

89. J. Schnack, Dalton Trans. 39, 4677 (2010). 

90. J. Vidal, R. Mosseri and B. Dougot, Phys. Rev. Lett. 81, 5888 (1998); J. Vidal, 
B. Dougot, R. Mosseri and P. Butaud, Phys. Rev. Lett. 85, 3906 (2000); J. Vidal, 
P. Butaud, B. Dougot and R. Mosseri, Phys. Rev. B64, 155306 (2001); see also B. 
Dougot and J. Vidal, Phys. Rev. Lett. 88, 227005 (2002); B. Dougot, L. B. Ioffe and 
J. Vidal, Phys. Rev. B69, 214501 (2004). 

91. G. Moller and N. R. Cooper, Phys. Rev. Lett. 108, 045306 (2012). 

92. U. Fano, Phys. Rev. 124, 1866 (1961). 

93. A. E. Miroshnichenko, S. Flach and Y. S. Kivshar, Rev. Mod. Phys. 82, 2257 (2010). 








April 14, 2015 0:38 WSPC/INSTRUCTION FILE de'rimarevised 


70 O. Derzhko, J. Richter, M. Maksymenko 

94. Z. Nussinov and J. van den Brink, Rev. Mod. Phys. 87, 1 (2015). 

95. H.-J. Schmidt, J. Phys. A35, 6545 (2002). 

96. H.-J. Schmidt, J. Richter and R. Moessner, J. Phys. A39, 10673 (2006). 

97. R. J. Baxter, Exactly Solved Models in Statistical Mechanics (Academic Press, Lon¬ 
don, 1982). 

98. E. Miiller-Hartmann and J. Zittartz, Z. Phys. B27, 261 (1977); X. N. Wu and F. Y. 
Wu, Phys. Lett. A144, 123 (1990); X.-Z. Wang and J. S. Kim, Phys. Rev. Lett. 78, 
413 (1997); S. J. Penney, V. K. Cumyn and D. D. Betts, Physica A330, 507 (2003). 

99. J.-B. Fouet, F. Mila, D. Clarke, H. Youk, O. Tchernyshyov, P. Fendley and R. 
M. Noack, Phys. Rev. B73, 214405 (2006); J.-D. Picon, A. F. Albuquerque, K. P. 
Schmidt, N. Laflorencie, M. Troyer and F. Mila, Phys. Rev. B78, 184418 (2008). 

100. A. Honecker and A. Lauchli, Phys. Rev. B63, 174407 (2001). 

101. D. J. Klein, J. Chem. Phys. 61, 786 (1974). 

102. P. Fulde, Electron Correlations in Molecules and Solids (Springer-Verlag, Berlin, Hei¬ 
delberg, 1993). 

103. F. H. L. Essler, H. Frahm, F. Gohmann, A. Kliimper and V. E. Korepin, The One- 
Dimensional Hubbard Model (Cambridge University Press, Cambridge, UK, 2005). 

104. V. L. Berezinskii, Zh. Eksp. Teor. Fiz. 59, 907 (1970) [Sov. Phys. JETP 32, 493 
(1971)]; J. M. Kosterlitz and D. J. Thouless, J. Phys. C6, 1181 (1973). 

105. R. Gupta, J. DeLapp, G. G. Batrouni, G. C. Fox, C. F. Baillie and J. Apostolakis, 
Phys. Rev. Lett. 61, 1996 (1988); R. Gupta and C. F. Baillie, Phys. Rev. B45, 2883 
(1992); M. Hasenbusch and K. Pinn, J. Phys. A30, 63 (1997); M. Hasenbusch, J. 
Phys. A38, 5869 (2005); Y. Komura and Y. Okabe, J. Phys. Soc. Jpn. 81, 113001 
(2012); Y.-D. Hsieh, Y.-J. Kao and A. W. Sandvik, J. Stat. Mech. P09001 (2013). 

106. H.-Q. Ding and M. S. Makivic, Phys. Rev. B42, 6827 (1990); H.-Q. Ding, Phys. Rev. 
B45, 230 (1992); K. Harada and N. Kawashima, Phys. Rev. B55, R11949 (1997); K. 
Harada and N. Kawashima, J. Phys. Soc. Jpn. 67, 2768 (1998); J. Carrasquilla and 
M. Rigol, Phys. Rev. A86, 043629 (2012); G. Ceccarelli, J. Nespolo, A. Pelissetto 
and E. Vicari, Phys. Rev. B88, 024517 (2013). 

107. F. Hebert, G. G. Batrouni, R. T. Scalettar, G. Schmid, M. Troyer and A. Dorneich, 

Phys. Rev. B65, 014513 (2001); G. Schmid, S. Todo, M. Troyer and A. Dorneich, 

Phys. Rev. Lett. 88, 167208 (2002). 

108. U. Tutsch, B. Wolf, S. Wessel, L. Postulka, Y. Tsui, H. O. Jeschke, I. Opahle, T. 
Saha-Dasgupta, R. Valenti, A. Briihl, K. Removic-Langer, T. Kretz, H.-W. Lerner, 
M. Wagner and M. Lang, Nat. Commun. 5, 5169 (2014). 

109. J. Ruostekoski, Phys. Rev. Lett. 103, 080406 (2009); B. Damski, H. Fehrmann, H.-U. 
Everts, M. Baranov, L. Santos and M. Lewenstein, Phys. Rev. A72, 053612 (2005). 

110. K. Awaga, T. Okuno, A. Yamaguchi, M. Hasegawa, T. Inabe, Y. Maruyama and N. 
Wada, Phys. Rev. B49, 3975 (1994). 

111. S. M. Girvin and T. Jach, Phys. Rev. B29, 5617 (1984). 

112. L. G. Phillips, G. De Chiara, P. Ohberg and M. Valiente, Phys. Rev. B91, 054103 

(2015). 

113. M. Tovmasyan, E. P. L. van Nieuwenburg and S. D. Huber, Phys. Rev. B88, 
220510(R) (2013). 

114. S. Takayoshi, H. Katsura, N. Watanabe and H. Aoki, Phys. Rev. A88, 063613 (2013). 

115. K. Hida, J. Phys. Soc. Jpn. 70, 3673 (2001). 

116. A. Honecker, J. Schulenburg and J. Richter, J. Phys.: Condens. Matter 16, S749 
(2004). 

117. D. C. Cabra, M. D. Grynberg, P. C. W. Holdsworth, A. Honecker, P. Pujol, J. Richter, 
D. Schmalfufi and J. Schulenburg, Phys. Rev. B71, 144420 (2005); A. Honecker, D. 








April 14, 2015 0:38 WSPC/INSTRUCTION FILE de'rimarevised 


Strongly correlated flat-band systems: The route from Heisenberg spins to Hubbard electrons 71 

C. Cabra, M. D. Grynberg, P. C. W. Holdsworth, P. Pujol, J. Richter, D. Schmalfufi 
and J. Schulenburg, Physica B359, 1391 (2005). 

118. Y. Okamoto, M. Tokunaga, H. Yoshida, A. Matsuo, K. Kindo and Z. Hiroi, Phys. 
Rev. B83, 180407(R) (2011). 

119. Y. Inagaki, Y. Narumi, K. Kindo, H. Kikuchi, T. Kamikawa, T. Kunimoto, S. Okubo, 
H. Ohta, T. Saito, M. Azuma, M. Takano, H. Nojiri, M. Kaburagi and T. Tonegawa, 
J. Phys. Soc. Jpn. 74, 2831 (2005). 

120. T. Hamada, J.-i. Kane, S.-i. Nakagawa and Y. Natsume, J. Phys. Soc. Jpn. 57, 1891 
(1988). 

121. W. Selke, Z. Phys. B27, 81 (1977); M. Yamada and M. Takahashi, J. Phys. Soc. 
Jpn. 55, 2024 (1986); M. Hartel, J. Richter, D. Ihle and S.-L. Drechsler, Phys. Rev. 
B78, 174412 (2008); J. Sirker, V. Y. Krivnov, D. V. Dmitriev, A. Herzog, O. Janson, 
S. Nishimoto, S.-L. Drechsler and J. Richter, Phys. Rev. B84, 144403 (2011). 

122. J. Richter, O. Krupnitska, T. Krokhmalskii and O. Derzhko, J. Magn. Magn. Mat. 
379, 39 (2015). 

123. V. Derzhko and J. J§drzejewski, “How to recognize a nearly-flat-band ferromagnet 
by means of thermodynamic measurements?”, arXiv:1004.2786 

124. D. Stauffer, Physics Reports 54, 1 (1979). 

125. D. Stauffer and A. Aharony, Introduction to Percolation Theory (Taylor & Francis, 
1994). 

126. M. B. Isichenko, Rev. Mod. Phys. 64, 961 (1992). 

127. D. Achlioptas, R. M. D’Souza and J. Spencer, Science 323, 1453 (2009); Y. S. Clio, 
S. Hwang, H. J. Herrmann and B. Kahng, Science 339, 1185 (2013). 

128. S. N. Dorogovtsev, Lectures on Complex Networks (Oxford University Press, Oxford, 

2010 ). 

129. M. Ostilli, Physica A391, 3417 (2012). 

130. M. E. Fisher and J. W. Essam, J. Math. Phys. 2, 609 (1961). 

131. M. Mezard and G. Parisi, Eur. Phys. J. B20, 217 (2001); C. Laumann, A. Scardicchio 
and S. L. Sondhi, Phys. Rev. B78, 134424 (2008); O. Rivoire, J. Stat. Mech. P07004 
(2005); B. Hsu, C. Laumann, A. Lauchli, R. Moessner and S. Sondhi, Phys. Rev. 
A87, 062334 (2013). 

132. J. Barre and B. Gongalves, Physica A386, 212 (2007). 

133. M. F. Thorpe, NATO Advanced Study Institute Series B78, 85 (1982). 

134. C. R. Laumann, S. A. Parameswaran and S. L. Sondhi, Phys. Rev. B80, 144415 
(2009). 

135. A. Weinrib, Phys. Rev. B29, 387 (1984). 

136. J. Hoshen and R. Kopelman, Phys. Rev. B14, 3438 (1976). 

137. M. E. J. Newman and R. M. Ziff, Phys. Rev. E64, 016706 (2001). 

138. K. Hukushima and K. Nemoto, J. Phys. Soc. Jpn. 65, 1604 (1996). 

139. J. Villain, R. Bidaux, J. P. Carton and R. Conte, J. Phys. 41, 1263 (1980). 

140. E. F. Shender, Zh. Eksp. Teor. Fiz. 83, 326 (1982) [Sot;. Phys. JETP 56, 178 (1982)]. 

141. P. Muller, J. Richter and O. Derzhko, in preparation. 

142. P. W. Anderson, Phys. Rev. 109, 1492 (1958); see also B. Kramer and A. MacKinnon, 
Rep. Prog. Phys. 56 1469 (1993) and references therein. 

143. K. Ishii, Prog. Theor. Phys. Suppl. 53, 77 (1973). 

144. H. Ishii and T. Nakayama, Phys. Rev. B73, 235311 (2006); M. Dey, S. K. Maiti and 
S. N. Karmakar, J. Appl. Phys. 110, 094306 (2011). 

145. Y. Imry and R. Landauer, Rev. Mod. Phys. 71, S306 (1999). 

146. S. Datta, Electronic Transport in Mesoscopic Systems (Cambridge University Press, 
Cambridge UK, 1995). 








April 14, 2015 0:38 WSPC/INSTRUCTION FILE de'rimarevised 


72 O. Derzhko, J. Richter, M. Maksymenko 

147. T. Enss, V. Meden, S. Andergassen, X. Barnabe-Theriault, W. Metzner and K. 
Schonhammer, Phys. Rev. B71, 155401 (2005). 

148. Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 (1992). 

149. H. Ishida and A. Liebsch, Phys. Rev. B86, 205115 (2012). 

150. P.-B. Niu, Y.-Y. Zhang, Q. Wang and Y.-H. Nie, Phys. Lett. A376, 1481 (2012). 

151. Y. Hatsugai, K. Shiraishi and H. Aoki, “Flat bands in Weaire-Thorpe model and 
silicene”, arXiv: 1410.7885 

152. S. Miyahara, S. Kusuta and N. Furukawa, Physica C460-462, 1145 (2007); V. I. 
Iglovikov, F. Hebert, B. Gremaud, G. G. Batrouni and R. T. Scalettar, Phys. Rev. 
B90, 094506 (2014). 

153. G. Zhu, J. Koch and I. Martin, “Nematic superfluidity and Wigner crystallization of 
bosons in frustrated kagome lattices”, arXiv: 1411.0043 

154. M. Goda, R. Tanaka and S. Nishino, physica status solidi (c) 1, 3043 (2004). 

155. M. Udagawa, H. Ishizuka and Y. Motome, Phys. Rev. Lett. 104, 226405 (2010); H. 
Ishizuka, M. Udagawa and Y. Motome, Phys. Rev. B83, 125101 (2011). 

156. A. Nandy, B. Pal and A. Chakrabarti, J. Phys.: Condens. Matter 27, 125501 (2015). 

157. M. Biondi, E. P. L. van Niewenburg, G. Blatter, S. D. Huber and S. Schmidt, “In¬ 
compressible polaritons in a flat band”, arXiv: 1502.07854 








